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HIGH THROUGHPUT DIRECTED EVOLUTION BY RATIONAL 

MUTAGENESIS 

' * 

RELATED APPLICATIONS 

Benefit of priority is claimed, to U.S. application Serial No. 
10/022,249, filed December 17, 2001 , to Manuel Vega and Lila Drittanti, 
entitled "HIGH THROUGHPUT DIRECTED EVOLUTION BY RATIONAL 
5 MUTAGENESIS." Benefit of priority is also claimed to U.S. provisional 
patent application serial No. 60/315,382, filed August 27, 2001, to 
Manuel Vega, Lila Drittanti and Marjorie Flaux, entitled "HIGH THROUGH- 
PUT DIRECTED EVOLUTION BY RATIONAL MUTAGENESIS." Where 
permitted, the subject matter of each of these applications is incorporated 
10 in its entirety by reference thereto. 
FIELD OF INVENTION 

Processes and systems for the high throughput directed evolution 
of peptides and proteins, particularly those that act in complex biological 
settings, are provided. The proteins and peptides include, but are not 
15 limited to, intracellular proteins, messenger/signaling/hormone proteins 

and viral proteins. 
BACKGROUND 

Directed evolution refers to biotechnological processes for 
optimizing the activity of proteins by means of random changes 

20 introduced into selected respective genes. Directed evolution involves the 
creation of a library of mutated genes, and then selection of the mutants 
that encode proteins having desired properties. The process can be an 
iterative one in which gene products that have improvement in a desired 
property are subjected to further cycles of mutation and screening. 

25 Directed evolution provides a way to adapt natural proteins to work in 
new chemical or biological environments, and/or to elicit new functions. 
The .potential plasticity of proteins is such that chances exist that for 
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every new challenge, such as a new environment and desired new or 
altered activity, it should be possible, given a sufficient pool of modified 
proteins (or encoding nucleic acids), that an appropriately 'evolved' 
protein could be found that would have a desired activity. The problem is 
5 in generating and then identifying the appropriate sequence. 

There have been practical approaches to this problem (see, e.g., 
U.S. Patent Nos. 6,096,548; 6,117,679; 6,165,793; 6,180,406; 
6,132,970; 6,171,820; 6,238,884; 6,174,673; 6,057,103; 6,001,574; 
5,763,239; 5,837,500; 5,571,698; 6,156,509; 5,723,323; 5,862,514; 

.10 5,871,974; 5,779,434 and others). Typically theses approaches are of 
two types. One is a purely "rational" approach that is based on the 
assumption that the optimized proteins can be rationally designed. This, 
however, requires sufficient information regarding the laws that govern 
protein folding, molecular interactions, intra-molecular forces and other 

15 dynamics of protein activity. This rational approach is extremely 

dependent on a number of variables and parameters that are not known. 
Consequently, although useful in some specific cases and applications, 
the rational approach intended to 'predict' protein structure remains 
limited in applicability. 

20 In contrast to the rational approach, random approaches have also 

been employed. One random approach requires synthesis of all possible 
protein sequences or a statistically sufficient large number of proteins and 
then screening them to identify proteins having the desired activity or 
property. Since the resources to synthesize all possible theoretical 

25 sequences of a single protein is not possible, this approach is 

impracticable. Other random approaches are based on gene shuffling 
methods, which are PCR-based methods that generate random 
rearrangements between two or more sequence-related genes to 
randomly generate variants of the gene. 
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The development and scope of directed evolution, thus, has been 
limited, and its potential remains to be exploited. In order to exploit the 
potential of directed evolution, alternative approaches for generating and 
identifying evolved proteins are needed. It is an object herein to provide 
5 methods and products to exploit the potential of directed evolution. 
SUMMARY 

Provided herein are methods for performing directed evolution for 
the optimization of proteins that function in complex biological settings. 
Methods of high throughput directed evolution of proteins are provided. 
10 In practicing the methods, each molecule is individually designed, 

produced, processed, screened and tested in a high throughput format. 
Neither random or combinatorial methods nor mixtures of molecules are 
used. 

The methods provided herein include the steps of identifying a 
15 protein target of interest; obtaining nucleic acids that encode the target, 
which may be from any source, such as a natural library, a collection 
generated by known gene shuffling techniques and related methods, and, 
then creating variants of the proteins using methods for rational 
mutagenesis provided herein. Whatever method is used to select or 
20 generate the nucleic acids encoding the protein targets, each molecule is 
processed and screened separately in a high throughput format. 

The nucleic acids encoding each variant are individually screened. 
They can be screened in any suitable assay, including cell-based assays 
and biochemical assays. For cell based assays, each nucleic acid 
25 molecule is introduced into an expression vector for expression in a 

bacterial cell or into a vector for expression in a eukaryotic host cell. In all 
instances, the nucleic acids of interest are introduced into host cells in an 
expression vector, such as by transfection for bacterial hosts and 
transduction with viral vectors into eukaryotic hosts with viral vectors. 
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Each variant is introduced into a host and the resulting cells are 
maintained separately, such as in an addressable array of wells in a 
microtiter plate or other substrate with discrete locations for performing 
reactions or retaining molecules of interest. Typical formats are 96 loci, 
5 and multiples thereof (384, 1536, 3072, . . . 96 x n, where n is 1 to any 
number desired, such as 10, 20, 30, 50 . . .100), although any 
convenient number of loci may be employed. 

Since the process is conducted in a high throughput format, for 
many embodiments, it is often important to assess the relative numbers 

10 of transformed, transduced or transfected cells. Hence the relative (or 
actual) titer of the vector, such as the recombinant viral vector, must be 
known to permit analysis of results. For high throughput formats, it is 
important to assess the relative or actual concentration of the viral vector 
(or plasmid) so that results can be compared among all cells and variants. 

15 Methods for titering (determining the concentration) of the nucleic acid 
encoding the variant and/or the recombinant virus are also provided. 

The processes require accurate titering of the viruses in a collection 
or among collections (libraries) so that the activities of the screened 
mutant proteins can be compared. Provided are general methods for the 

20 quantitative assessment of the parameters of activity corresponding to 

the individual variants in the library, based upon intracellular serial dilution 
generated by precise titering with the gene transfer viral vectors. Any 
method that permits accurate titering may be used, including that 
described in International PCT application No. PCT/FR01/01 366, based on 

25 French application n° 0005852, filed 9 May 2000, and published as 

International PCT application No. WO 01/186291 . A method of titering, 
designated Tagged Replication and Expression Enhancement Technology 
(TREE") is provided herein. 
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Each of the different cells is separately screened by a suitable 
assay, and the results analyzed. Methods for assessing the interactions 
in biological systems, such as a Hill-based analysis (see, published 
International PCT application No. WO 01/44809 based on PCT n° 
5 PCT/FROO/03503, Dec, 2000, and the description herein), or a second 
order polynomial or other algorithm that describes the interaction between 
cells and biological agents to select variants that have a desired property 
are employed in the processes herein. 

A semi-rational method for evolution of proteins that is particularly 

10 designed for use in the methods herein or in any method that uses 

"evolved" proteins is also provided. The method, which is based on an 
amino-acid scanning protocol, is for rationally designing the variants for 
use in the directed evolution and selection method, and can employ 
iterative processing of the steps of the high throughput methods provided 

15 herein. In this method, once the target protein or domain is identified, 
nucleic acid molecules encoding variants are prepared. Each variant 
encoded by the nucleic acid molecules has a single amino acid replaced 
with another selected amino acid, such as alanine (Ala), glycine (Gly), 
serine (Ser) or any other suitable amino acid, typically one selected to 

20 have a neutral effect on secondary and tertiary structure. The resulting 
series of variants are separately screened in the high throughput format 
provided herein, and those that have a change in the target activity are 
selected and the modified amino acids are designated "hits." Nucleic acid 
molecules encoding proteins in which each hit position is replaced by the 

25 eighteen remaining amino acids then are synthesized and the resulting 
collection of molecules is screened, such by introduction into host cells, 
and the proteins that result in an improvement of a targeted activity, are 
identified. Such proteins are designated "leads." Leads may be further 
modified by producing proteins that have combinations of the mutations 
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identified in the leads. This method, which does not require any 
knowledge of the structure of a target protein, permits precise control of 
locations where changes are introduced and also the amount of change 

that is introduced. - 
5 The high throughput directed evolution processes provided herein 

include the use of virus libraries containing mutant versions of a gene; 
viral libraries of such mutant genes are also provided. 

Reporter cells are infected with the titered viruses that encode the 
mutant genes. The mutant genes are expressed and read-out data from 

10 either biochemical or cell-based assays, while isolating each mutant/virus 
physically from the others (i.e. one-by-one analysis), are collected and 
analyzed. Serial dilution assays (i.e. a series of dilutions for each 
individual mutant/virus in the library) are used and the biochemical/cell- 
based assays are performed on each single dilution for each individual 

15 mutant/virus. Analysis of the serial dilution readout-data can be 
performed using any method of analysis that permits one-by-one 
comparisons. Hill-based analysis (see, published International PCT 
application No. WO 01/44809 based on PCT n° PCT/FRO0/035O3, Dec, 
2000, and the description herein) are employed for analysis of the data. 

20 Protein/protein domain variants identified using the methods are 

also provided. Also provided are nucleic acid molecules and proteins and 
polypeptides produced by the methods and viruses and cells that contain 
the nucleic acid molecules and proteins. 

In an exemplary embodiment of methods provided herein, the 

25 process of rational directed evolution provided herein is applied to the 
AAV rep gene. The resulting recombinant rep protein variants and rAAV 
are also provided. Among the rep proteins are those that result in 
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increased rAAV production in rAAV that encode such mutants, thereby, 
among a variety of advantages, offering a solution to the need in the gene 
therapy industry to increase the production therapeutic vectors without 
up-scaling manufacturing. 
5 Thus, for exemplification, some methods provided herein have been 

used to identify amino acid "hit" positions in adeno-associated virus 
(AAV) rep proteins that are relevant for AAV or rAAV production. Those 
amino acid positions are such that a change in the amino acid leads to a 
change in protein activity either to lower activity or to higher activity 
10 compared to native-sequence Rep proteins. The hit positions were then 
used to generate further mutants designated "leads." Provided herein are 
the resulting mutant rep proteins that result in either higher or lower 
levels of AAV or rAAV virus compared to the wild-type (native) Rep 
protein(s). 

15 In addition to enhancing AAV production, among the rep mutants 

are those that inhibit papillomavirus (PV) and PV-associated diseases, 
including certain cancers and human immunodeficiency virus (HIV) and 

HIV-associated diseases. 

Systems and computer controlled systems for performing the high 

20 throughput processes are also provided. 
DESCRIPTION OF THE FIGURES 

FIGURES 1 A - 1E summarize various exemplary embodiments of 
the high throughput processes provided herein. FIGURE 1A depicts an 
embodiment of the process in which an amino acid scan is employed to 
25 generate a library of mutants, which are then introduced into viral 

vectors, such as an adeno-associated viral vector (AAV), a herpes virus, 
such as herpes simplex virus (HSV) and other herpes virus vectors, a 
vaccinia virus vector, retroviral vectors, such as MuMLV, MoMLV, feline 
leukemia virus, and HIV and other lentiviruses, adenovirus vectors and 
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other suitable viral vector, each member of the library is individually 
tested and phenotypically characterized to identify HITS. FIGURE 1 B 
summarizes round 2 in which LEADS are developed by mutagenesis at 
and/or surrounding the positions identified as HITS; FIGURE 1C 
5 summarizes the optional next round in which recombination among 
LEADS is performed to further optimize the LEADS; FIGURE ID depicts 
the process in mammalian cells; and FIGURE 1E depicts the process in 
bacterial cells. 

FIGURE 2A depicts an exemplary titering process (in this instance 
10 the TREE" for titering AAV) in a 96 well format; FIGURE 2B shows the 
results and analysis of a titering process performed using the TREE" 
procedure; and FIGURE 2C shows an exemplary calibration curve for the 
calculation of the titer using the TREE" method. 

FIGURE 3A and 3B depict "HITS" and "LEADS" respectively for 
15 identification of AAV rep mutants "evolved" for increased activity. 

FIGURE 4 shows the genetic map of AAV, including the location of 
promoters, and transcripts; amino acid 1 of the Rep 78 gene is at 
nucleotide 321 in the AAV-2 genome. 

FIGURES 5A and 5B show the alignment of amino acid sequences 
20 of Rep78 among AAV-1; AAV-6; AAV-3; AAV-3B; AAV-4; AAV-2; AAV- 
5 sequences, respectively; the hit positions with 100 percent homology 
among the serotypes are bolded italics, where the position is different 
(compared to AAV-2, no. 6 in the Figure) in a particular serotype, it is in 
bold; a sequence indicating relative conservation of sequences among 
25 the serotypes is labeled "C". 
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Legend: 

1 is AAV-1; 2 is AAV-6, 3 is AAV-3, 4 is AAV-3B, 
5 is AAV-4, 6 is AAV-2, and 7 is AAV-5; 
"." where the amino acid is present > 20%; 
5 ":" where the amino acid is present > 40%; 

" + » where the amino acid is present > 60%; 
"*" where the amino acid is present >l 80%; and 
where the amino acid is the same amongst all 
serotypes depicted it is represented by its single letter 
10 code. 

DETAILED DESCRIPTION 
A. Definitions 

Unless defined otherwise, all technical and scientific terms used 
herein have the same meaning as is commonly understood by one of skill 

15 in the art to which this invention belongs. All patents, patent 

applications, published applications and publications, Genbank sequences, 
websites and other published materials referred to throughout the entire 
disclosure herein are, unless noted otherwise, incorporated by reference 
in their entirety. In the event that there are a plurality of definitions for 

20 terms herein, those in this section prevail. 

As used herein, directed evolution refers to methods that adapt 
natural proteins or protein domains to work in new chemical or biological 
environments and/or to elicit new functions. It is more a more broad- 
based technology than DNA shuffling. 

25 As used herein, high-throughput screening (HTS) refers to 

processes that test a large number of samples, such as samples of test 
proteins or cells containing nucleic acids encoding the proteins of interest 
to identify structures of interest or to identify test compounds that 
interact with the variant proteins or cells containing them. HTS 

30 operations are amenable to automation and are typically computerized to 
handle sample preparation, assay procedures and the subsequent 
processing of large volumes of data. 
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As used herein, DNA shuffling is a PCR-based technology that 
produces random rearrangements between two or more sequence-related 
. genes to generate related, although different, variants of given gene. 

As used herein, "hits" are mutant proteins that have an alteration in 
5 any attribute, chemical, physical or biological property in which such 
alteration is sought. In the methods herein, hits are generally generated 
by systematically replacing each amino acid in a protein or a domain 
thereof with a selected amino acid, typically Alanine, Glycine, Serine or 
any amino acid, as long as each residue is replaced with the same 

10 residue. Hits may be generated by other methods known to those of skill 
in the art and tested by the high throughput methods herein. For 
purposes herein a Hit typically has activity with respect to the function of 
interest that differs by at least 10%, 20%, 30% or more from the wild 
type or native protein. The desired alteration, which is generally a 

15 reduction in activity, will depend upon the function or property of interest. 

As used herein, "leads" are "hits" whose activity has been 
optimized for the particular attribute, chemical, physical or biological 
property. In the methods herein, leads are generally produced by 
systematically replacing the hit loci with all remaining 1 8 amino acids, and 

20 identifying those among the resulting proteins that have a desired activity. 
The leads may be further optimized by replacement of a plurality of "hit" 
residues. Leads may be generated by other methods known to those of 
skill in the art and tested by the high throughput methods herein. For 
purposes herein a lead typically has activity with respect to the function 

25 of interest that differs from the native activity, by a desired amount and is 
at by at least 10%, 20%, 30% or more from the wild type or native 
protein. Generally a Lead will have an activity that is 2 to 10 or more 
times the native protein for the activity of interest. As with hits, the 
change in the activity is dependent upon the activity that is "evolved." 
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The desired alteration will depend upon the function or property of 
interest. 

As used herein, MOI is multiplicity of infection. 
As used herein, ip, with reference to a virus or recombinant vector, 
5 refers to a titer of infectious particles. 

As used herein, pp refers to the total number of vector {or virus) 

physical particles 

As used herein, biological and pharmacological activity includes any 
activity of a biological pharmaceutical agent and includes, but is not 

10 limited to, biological efficiency, transduction efficiency, gene/transgene 
expression, differential gene expression and induction activity, titer, 
progeny productivity, toxicity, cytotoxicity, immunogenicity, cell 
proliferation and/or differentiation activity, anti-viral activity, 
morphogenetic activity, teratogenetic activity, pathogenetic activity, 

15 therapeutic activity, tumor suppressor activity, ontogenetic activity, 
oncogenetic activity, enzymatic activity, pharmacological activity, 
cell/tissue tropism and delivery. 

As used herein, "output signal" refers to parameters that can be 
followed over time and, if desired, quantified. For example, when a virus 

20 infects a cell, the infected cell undergoes a number of changes. Any 
such change that can be monitored and used to assess infection, is an 
"output signal," and the cell is referred to as a "reporter cell." Output 
signals include, but are not limited to, enzyme activity, fluorescence, 
luminescence, amount of product produced and other such signals. 

25 Output signals include expression of a viral gene or viral gene product, 
including heterologous genes (transgenes) inserted into the virus. Such 
expression is a function of time ("t") after infection, which in turn is 
related to the amount of virus used to infect the cell, and, hence, the 
concentration of virus ("s") in the infecting composition. For higher 
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concentrations the output signal is higher. For any particular 
concentration, the output signal increases as a function of time until a 
plateau is reached. Output signals may also measure the interaction 
between cells, expressing heterologous genes, and biological agents 
5 As used herein, adeno-associated virus (AAV) is a defective and 

non-pathogenic parvovirus that requires co-infection with either 
adenovirus or herpes virus for its growth and multiplication, able of 
providing helper functions. A variety of serotypes are known, and 
contemplated herein. Such serotypes include, but are not limited to: 

10 AAV-1 (Genbank accession no. NC002077; accession no. VR-645); AAV- 
2 (Genbank accession no. NC001401; accession no. VR-680); AAV-3 
(Genbank accession no. NC001729; accession no. VR-681); AAV-3b 
(Genbank accession no. NC001863); AAV-4 (Genbank accession no. 
NC001829; ATCC accession no. VR-646 ); AAV-6 (Genbank accession 

15 no.NC001862); and avian associated adeno-virus (ATCC accession no. 
VR-1449). The preparation and use of AAVs as vectors for gene 
expression in vitro and for in vivo use for gene therapy are well known 
(see, e.g., U.S. Patent Nos. 4,797,368, 5,139,941, 5,798,390 and 
6,127,175; Tessier et al. (2001) J. Virol. 75:375-383; Salvetti et al. 

20 (1 998) Hum Gene Ther 20:695-706; Chadeuf et al. (2000) J Gene Med 
2:260-268). 

As used herein, the activity of a Rep protein or of a capsid protein 
refers to any biological activity that can be assessed. In particular, 
herein, the activity assessed for the rep proteins is the amount (i.e., titer) 
25 of AAV produced by a cell. 

As used herein, the Hill equation is a mathematical model that 
relates the concentration of a drug (i.e., test compound or substance) to 
the response being measured 



30 
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y = [D] n + [D 50 ] n 

where y is the variable being 
5 measured, such as a response, 

signal, y max is the maximal response achievable, [D] is the molar 
concentration of a drug, [D 50 ] is the concentration that produces a 50% 
maximal response to the drug, n is the slope parameter, which is 1 if the 
drug binds to a single site and with no cooperativity between or among 
10 sites. A Hill plot is Iog 10 of the ratio of ligand-occupied receptor to free 
receptor vs. log [D] (M). The slope is n, where a slope of greater than 1 
indicates cooperativity among binding sites, and a slope of less than 1 
can indicate heterogeneity of binding. This general equation has been 
employed for assessing interactions in complex biological systems (see, 
15 published International PCT application No. WO 01/44809 based on PCT 
n° PCT/FROO/03503, see, also, EXAMPLES). 

As used herein, in the Hill-based analysis {published International 
PCT application No. WO 01/44809 based on PCT n° PCT/FR00/03503), 
the parameters, h f K,T,€,rj t O t are as follows: 
20 it potency of the biological agent acting on the assay (cell- 

based) system; 

k constant of resistance of the assay system to elicit a 
response to a biological agent; 

e is global efficiency of the process or reaction triggered by the 
25 biological agent on the assay system; 

t is the apparent titer of the biological agent; 
9 is the absolute titer of the biological agent; and 
q is the heterogeneity of the biological process or reaction. 
In particular, as used herein, the parameters n (potency) or k 
30 (constant of resistance) are used to respectively assess the potency of a 
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test agent to produce a response in an assay system and the resistance 
of the assay system to respond to the agent. 

As used herein, ^{efficiency), is the slope at the inflection point of 
the Hill curve (or, in general, of any other sigmoidal or linear 
5 approximation), to assess the efficiency of the global reaction (the 
biological agent and the assay system taken together) to elicit the 
biological or pharmacological response. 

As used herein, r (apparent titer) is used to measure the limiting 
dilution or the apparent titer of the biological agent. 
10 As used herein, 9 (absolute titer), is used to measure the absolute 

limiting dilution or titer of the biological agent. 

As used herein, n (heterogeneity) measures the existence of 
discontinuous phases along the global reaction, which is reflected by an 
abrupt change in the value of the Hill coefficient or in the constant of 

15 resistance. 

As used herein, a library of mutants refers to a collection of 

plasmids or other vehicles that carry (encode) the gene variants, such that 

individual plasmids or other vehicles carry individual gene variants. 

When a library of proteins is contemplated, it will be so-stated. 
20 As used herein, a "reporter cell" is the cell that "reports", i.e., 

undergoes the change, in response to introduction of the nucleic acid 

infection and, therefore, it is named here a reporter cell. 

As used herein, "reporter" or "reporter moiety" refers to any moiety 

that allows for the detection of a molecule of interest, such as a protein 
25 expressed by a cell. Typical reporter moieties include, for example, 

fluorescent proteins, such as red, blue and green fluorescent proteins. 

For expression in cells, nucleic acid encoding the reporter moiety can be 

expressed as a fusion protein with a protein of interest or under the 

control of a promoter of interest. 
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As used herein, a titering virus increases or decreases the output 
signal from a reporter virus, which is a virus that can be detected, such 
as by a detectable label or signal. 

As used herein, phenotype refers to the physical or other 
5 manifestation of a genotype (a sequence of a gene). In the methods 

herein, phenotypes that result from alteration of a genotype are assessed. 

As used herein, activity refers to the function or property to be 
evolved. An active site refers to a site(s) responsible or that participates 
in conferring the activity or function. The activity or active site evolved 
10 (the function or property and the site conferring or participating in 

conferring the activity) may have nothing to do with natural activities of 
a protein. For example, it could be an 'active site' for conferring 
immunogenicity (immunogenic sites or epitopes) on a protein. 

As used herein, the amino acids, which occur in the various amino 
15 acid sequences appearing herein, are identified according to their known, 
three-letter or one-letter abbreviations (see, Table 1). The nucleotides, 
which occur in the various nucleic acid fragments, are designated with 
the standard single-letter designations used routinely in the art. 

As used herein, amino acid residue refers to an amino acid formed 
20 upon chemical digestion (hydrolysis) of a polypeptide at its peptide 

linkages. The amino acid residues described herein are presumed to be in 
the "L" isomeric form. Residues in the "D" isomeric form, which are so- 
designated, can be substituted for any L-amino acid residue, as long as 
the desired functional property is retained by the polypeptide. NH 2 refers 
25 to the free amino group present at the amino terminus of a polypeptide. 
COOH refers to the free carboxy group present at the carboxyl terminus 
of a polypeptide. In keeping with standard polypeptide nomenclature 
described in J. Biol. Chem., 243:3552-59 (1969) and adopted at 37 
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C.F.R. § § 1.821 - 1.822, abbreviations for amino acid residues are 
shown in the following Table: 

Table 1 



10 



15 



20 



25 



SYMBOL 




1 -Letter 


3-Letter 


AMINO ACID 


Y 


Tyr 


tyrosine 


G 


Gly 


glycine 


F 

t 


Phe 


phenylalanine 


M 

1 V 1 


Met 


^ methionine 


A 


Ala 


alanine 


s 


Ser 


serine 


I 
1 


He 


isoleucine 


1 


Leu 


leucine 


T 
i 


Thr 


threonine 


v 

V 


Val 

V 14 ■ 


valine 


p 


Pro 


proline 


K 


Lys 


lysine 


H 


His 


histidine 


Q 


Gin 


glutamine 


E 


Gtu 


glutamic acid 


Z 


Glx 


Glu and/or Gin 


w 


Trp 


tryptophan 


R 


Arg 


arginine 


D ■ 


Asp 


aspartic acid 


N 


Asn 


asparagine 


B 


Asx 


Asn and/or Asp 


C 


Cys 


cysteine 
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SYMBOL 




X 


Xaa 


Unknown or other 



It should be noted that all amino acid residue sequences 
represented herein by formulae have a left to right orientation in the 
5 conventional direction of amino-terminus to carboxyl-terminus. In 

addition, the phrase "amino acid residue" is broadly defined to include the 
amino acids listed in the Table of Correspondence and modified and 
unusual amino acids, such as those referred to in 37 C.F.R. § § 1.821- 
1.822, and incorporated herein by reference. Furthermore, it should be 

10 noted that a dash at the beginning or end of an amino acid residue 

sequence indicates a peptide bond to a further sequence of one or more 
amino acid residues or to an amino-terminal group such as NH 2 or to a 
carboxyl-terminal group such as COOH. 

In a peptide or protein, suitable conservative substitutions of amino 

15 acids are known to those of skill in this art and may be made generally 
without altering the biological activity of the resulting molecule. Those of 
skill in this art recognize that, in general, single amino acid substitutions 
in non-essential regions of a polypeptide do not substantially alter 
biological activity (see, e.g., Watson et al. Molecular Biology of the Gene, 

20 4th Edition, 1987, The Benjamin/Cummings Pub. co., p. 224). 

Such substitutions are preferably made in accordance with those 
set forth in TABLE 2 as follows: 
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TABLE 2 





Original residue 


Conservative substitution 




Ala (A) 


Gly; Ser 




Arg (R) 


Lys 


5 


Asn <N) 


Gin; His 




Cys (C) 


Ser 




Gin (Q) 


Asn 




Glu (E) 


Asp 




Gly (G) 


Ala; Pro 


10 


His (H) 


Asn; Gin 




He (I) 


Leu; Val 




Leu (L) 


lie; Val 




Lys (K) 


Arg; Gin; Glu 




Met (M) 


Leu; Tyr; He 


15 


Phe (F) 


Met; Leu; Tyr 




Pro (P) 


Ala; Gly 




Ser (S) 


Thr 




Thr (T) 


Ser 




Trp (W) 


Tyr 


20 


Tyr (Y) 


Trp; Phe 




Val (V) 


He; Leu 



Other substitutions are also permissible and may be determined 
empirically or in accord with known conservative substitutions. 

As used herein, nucleic acids include DNA, RNA and analogs 

25 thereof, including protein nucleic acids (PNA) and mixture thereof. 

Nucleic acids can be single or double stranded. When referring to probes 
or primers, optionally labeled, with a detectable label, such as a 
fluorescent or radiolabel, single-stranded molecules are contemplated. 
Such molecules are typically of a length such that they are statistically 

30 unique of low copy number (typically less than 5, preferably less than 3) 
for probing or priming a library. Generally a probe or primer contains at 
least 14, 16 or 30 contiguous of sequence complementary to or identical 
to a gene of interest. Probes and primers can be 10, 14, 16, 20, 30, 50, 
100 or more nucleic acid bases long. 

35 As used herein, homologous means about greater than 25% nucleic 

acid sequence identity, preferably 25% 40%, 60%, 80%, 90% or 95%. 
The intended percentage will be specified. The terms "homology" and 
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"identity" are often used interchangeably. In general, sequences are 
aligned so that the highest order match is obtained (see, e.g.: 
Computational Molecular Biology, Lesk, A.M., ed., Oxford University 
Press, New York, 1988; Biocomputing: Informatics and Genome Projects, 
5 Smith, D.W., ed., Academic Press, New York, 1993; Computer Analysis 
of Sequence Data, Part I, Griffin, A.M., and Griffin, H.G., eds., Humana 
Press, New Jersey, 1994; Sequence Analysis in Molecular Biology, von 
Heinje, G., Academic Press, 1987; and Sequence Analysis Primer, 
Gribskov, M. and Devereux, J., eds., M Stockton Press, New York, 1991; 

10 Carillo et al. (1988) S/AM J Applied Math 45:1073). By sequence 

identity, the number of conserved amino acids are determined by standard 
alignment algorithms programs, and are used with default gap penalties 
established by each supplier. Substantially homologous nucleic acid 
molecules would hybridize typically at moderate stringency or at high 

15 stringency all along the length of the nucleic acid of interest. Also 

contemplated are nucleic acid molecules that contain degenerate codons 
in place of codons in the hybridizing nucleic acid molecule. 

As used herein, a nucleic acid homolog refers to a nucleic acid that 
includes a preselected conserved nucleotide sequence, such as a 

20 sequence encoding a therapeutic polypeptide. By the term "substantially 
homologous", it is meant having at least 80%, preferably at least 90%, 
most preferably at least 95% homology therewith or a less percentage of 
homology or identity and conserved biological activity or function. 
The terms "homology" and "identity" are often used 

25 interchangeably. In this regard, percent homology or identity may be 

determined, for example, by comparing sequence information using a GAP 
computer program. The GAP program uses the alignment method of 
Needleman and Wunsch (J. MoL Biol. 48:443 (1970), as revised by Smith 
and Waterman (Adv. Appl. Math. 2:482 (1981). Briefly, the GAP program 
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defines similarity as the number of aligned symbols (i.e., nucleotides or 
amino acids) which are similar, divided by the total number of symbols in 
the shorter of the two sequences. The preferred default parameters for 
the GAP program may include: (1) a unitary comparison matrix 
5 (containing a value of 1 for identities and 0 for non-identities) and the 
weighted comparison matrix of Gribskov and Burgess, Nucl. Acids Res. 
14:6745 (1986), as described by Schwartz and Dayhoff, eds., ATLAS OF 
PROTEIN SEQUENCE AND STRUCTURE, National Biomedical Research 
Foundation, pp. 353-358 (1979); (2) a penalty of 3.0 for each gap and an 
10 additional 0.10 penalty for each symbol in each gap; and (3) no penalty 
for end gaps. 

Whether any two nucleic acid molecules have nucleotide sequences 
that are, for example, at least 80%, 85%, 90%, 95%, 96%, 97%, 98% 
or 99% /identical" can be determined using known computer algorithms 

15 such as the "FAST A" program, using for example, the default parameters 
as in Pearson and Lipman, Proc. Natl. Acad. ScL USA 55:2444 (1988). 
Alternatively the BLAST function of the National Center for Biotechnology 
Information database may be used to determine identity 

In general, sequences are aligned so that the highest order match 

20 is obtained. "Identity" per se has an art-recognized meaning and can be 
calculated using published techniques. (See, e.g. : Computational 
Molecular Biology, Lesk, A.M., ed., Oxford University Press, New York, 
1988; Biocomputing: Informatics and Genome Projects, Smith, D.W., ed., 
Academic Press, New York, 1993; Computer Analysis of Sequence Data, 

25 Part I, Griffin, A.M., and Griffin, H.G., eds., Humana Press, New Jersey, 
1994; Sequence Analysis in Molecular Biology, von Heinje, G., Academic 
Press, 1987; and Sequence Analysis Primer, Gribskov, M. and Devereux, 
J., eds., M Stockton Press, New York, 1991). While there exist a number 
of methods to measure identity between two polynucleotide or 
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polypeptide sequences, the term "identity" is well known to skilled 
artisans (Carillo, H. & Lipton, D., SIAMJ Applied Math 45:1073 (1988)). 
Methods commonly employed to determine identity or similarity between 
two sequences include, but are not limited to, those disclosed in Guide to 
5 Huge Computers, Martin J. Bishop, ed., Academic Press, San Diego, 
1994, and Carillo, H. & Lipton, D., SIAM J Applied Math 45:1073 
(1988). Methods to determine identity and similarity are codified in 
computer programs. Preferred computer program methods to determine 
identity and similarity between two sequences include, but are not limited 

10 to, GCG program package (Devereux, J., et al.. Nucleic Acids Research 
72(l):387 (1984)), BLASTP, BLASTN, FASTA (Atschul, S.F., et aL, J 
Molec Biol 2/5:403 (1990)), and CLUSTALW. For sequences displaying 
a relatively high degree of homology, alignment can be effected manually 
by simply lining up the sequences and manually or visually matching the 

15 conserved portions. 

Therefore, as used herein, the term "identity" represents a 
comparison between a test and a reference polypeptide or polynucleotide. 
For example, a test polypeptide may be defined as any polypeptide that 
is 90% or more identical to a reference polypeptide. 

20 For the alignments presented herein (see, Fig. 5) for the AAV 

serotype, the CLUSTALW program was employed with parameters set as 
follows: scoring matrix BLOSUM, gap open 10, gap extend 0.1, gap 
distance 40% and transitions/transversions 0.5; specific residue penalties 
for hydrophobic amino acids (DEGKNPQRS), distance between gaps for 

25 which the penalties are augmented was 8, and gaps of extremities 
penalized less than internal gaps. 

As used herein, a "corresponding" position on a protein, such as 
the AAV rep protein, refers to an amino acid position based upon 
alignment to maximize sequence identity. For AAV Rep proteins an 
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alignment of the Rep 78 protein from AAV-2 and the corresponding 
protein from other AAV serotypes (AAV-1, AAV-6, AAV-3, AAV-3B, 
AAV-4, AAV-2 and AAV-5) is shown in Figure 5. The "hit" positions are 
shown in italics. 

5 As used herein, the term at least "90% identical to" refers to 

* 

percent identities from 90 to 100% relative to the reference polypeptides. 
Identity at a level of 90% or more is indicative of the fact that/ assuming 
for exemplification purposes a test and reference polynucleotide length of 
100 amino acids are compared. No more than 10% (i.e., 10 out of 100) 

10 amino acids in the test polypeptide differs from that of the reference 
polypeptides. Similar comparisons may be made between a test and 
reference polynucleotides. Such differences may be represented as point 
mutations randomly distributed over the entire length of an amino acid 
sequence or they may be clustered in one or more locations of varying 

15 length up to the maximum allowable, e.g. 10/100 amino acid difference 
(approximately 90% identity). Differences are defined as nucleic acid or 
amino acid substitutions, or deletions. 

As used herein, it is also understood that the terms substantially 
identical or similar varies with the context as understood by those skilled 

20 in the relevant art. 

As used herein, genetic therapy involves the transfer of 
heterologous nucleic acids to the certain cells, target cells, of a mammal, 
particularly a human, with a disorder or conditions for which such therapy 
is sought. The nucleic acid, such as DNA, is introduced into the selected 

25 target cells in a manner such that the heterologous nucleic acid, such as 
DNA, is expressed and a therapeutic product encoded thereby is 
produced. Alternatively, the heterologous nucleic acid, such as DNA, 
may in some manner mediate expression of DNA that encodes the 
therapeutic product, or it may encode a product, such as a peptide or 
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RNA that in some manner mediates, directly or indirectly, expression of a 
therapeutic product. Genetic therapy may also be used to deliver nucleic 
acid encoding a gene product that replaces a defective gene or 
supplements a gene product produced by the mammal or the cell in which 
5 it is introduced. The introduced nucleic acid may encode a therapeutic 
compound, such as a growth factor or inhibitor thereof, or a tumor 
necrosis factor or inhibitor thereof, such as a receptor therefor, that is not 
normally produced in the mammalian host or that is not produced in 
therapeutically effective amounts or at a therapeutically useful time. The 

10 heterologous nucleic acid, such as DNA, encoding the therapeutic product 
may be modified prior to introduction into the cells of the afflicted host in 
order to enhance or otherwise alter the product or expression thereof. 
Genetic therapy may also involve delivery of an inhibitor or repressor or 
other modulator of gene expression. 

15 As used herein, heterologous or foreign nucleic acid, such as DNA 

and RNA, are used interchangeably and refer to DNA or RNA that does 
not occur naturally as part of the genome in which it is present or which 
is found in a location or locations in the genome that differ from that in 

» * 

which it occurs in nature. Heterologous nucleic acid is generally not 
20 endogenous to the cell into which it is introduced, but has been obtained 
from another cell or prepared synthetically. Generally, although not 
necessarily, such nucleic acid encodes RNA and proteins that are not 
normally produced by the cell in which it is expressed. Any DNA or RNA 
that one of skill in the art would recognize or consider as heterologous or 
25 foreign to the cell in which it is expressed is herein encompassed by 

heterologous DNA. Heterologous DNA and RNA may also encode RNA or 
proteins that mediate or alter expression of endogenous DNA by affecting 
transcription, translation, or other regulatable biochemical processes. 
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Examples of heterologous nucleic acid include, but are not limited to, 
nucleic acid that encodes traceable marker proteins, such as a protein 
that confers drug resistance, nucleic acid that encodes therapeutically 
effective substances, such as anti-cancer agents, enzymes and hormones, 
5 and DNA that encodes other types of proteins, such as antibodies. 

Hence, herein heterologous DNA or foreign DNA, includes a DNA 
molecule not present in the exact orientation and position as the 
counterpart DNA molecule found in the genome. It may also refer to a 
DNA molecule from another organism or species [i.e., exogenous). 

10 As used herein, a therapeutically effective product introduced by 

genetic therapy is a product that is encoded by heterologous nucleic acid, 
typically DNA, that, upon introduction of the nucleic acid into a host, a 
product is expressed that ameliorates or eliminates the symptoms, 
manifestations of an inherited or acquired disease or that cures the 

15 disease. 

As used herein, isolated with reference to a nucleic acid molecule 
or polypeptide or other biomolecule means that the nucleic acid or 
polypeptide has separated from the genetic environment from which the 
polypeptide or nucleic acid were obtained. It may also mean altered from 

20 the natural state. For example, a polynucleotide or a polypeptide naturally 
present in a living animal is not "isolated," but the same polynucleotide or 
polypeptide separated from the coexisting materials of its natural state is 
"isolated", as the term is employed herein. Thus, a polypeptide or 
polynucleotide produced and/or contained within a recombinant host cell 

25 is considered isolated. Also intended as an "isolated polypeptide" or an 
"isolated polynucleotide" are polypeptides or polynucleotides that have 
been purified, partially or substantially, from a recombinant host cell or 
from a native source. For example, a recombinantly produced version of 
a compounds can be substantially purified by the one-step method 
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described in Smith and Johnson, Gene 67:31-40 (1988). The terms 
isolated and purified are sometimes used interchangeably. 

Thus, by "isolated" it is meant that the nucleic is free of the coding 
sequences of those genes that, in the naturally-occurring genome of the 
5 organism (if any) immediately flank the gene encoding the nucleic acid of 
interest. Isolated DNA may be single-stranded or double-stranded, and 
may be genomic DNA, cDNA, recombinant hybrid DNA, or synthetic 
DNA. It may be identical to a native DNA sequence, or may differ from 
such sequence by the deletion, addition, or substitution of one or more 

10 nucleotides. 

Isolated or purified as it refers to preparations made from biological 
cells or hosts means any cell extract containing the indicated DNA or 
protein including a crude extract of the DNA or protein of interest- For 
example, in the case of a protein, a purified preparation can be obtained 

15 following an individual technique or a series of preparative or biochemical 
techniques and the DNA or protein of interest can be present at various 
degrees of purity in these preparations. The procedures may include for 
example, but are not limited to, ammonium sulfate fractionation, gel 
filtration, ion exchange change chromatography, affinity chromatography, 

20 density gradient centrifugation and electrophoresis. 

A preparation of DNA or protein that is "substantially pure" or 
"isolated" should be understood to mean a preparation free from naturally 
occurring materials with which such DNA or protein is normally 
associated in nature. "Essentially pure" should be understood to mean a 

25 "highly" purified preparation that contains at least 95% of the DNA or 
protein of interest. 

A cell extract that contains the DNA or protein of interest should be 
understood to mean a homogenate preparation or cell-free preparation 
obtained from cells that express the protein or contain the DNA of 
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interest. The term "cell extract" is intended to include culture media, 
especially spent culture media from which the cells have been removed. 

As used herein, receptor refers to a biologically active molecule 
that specifically binds to (or with) other molecules. The term "receptor 
5 protein" may be used to more specifically indicate the proteinaceous 
nature of a specific receptor. 

As used herein, recombinant refers to any progeny formed as the 
result of genetic engineering. 

As used herein, a promoter region refers to the portion of DNA of a 
10 gene that controls transcription of the DNA to which it is operatively 

linked. The promoter region includes specific sequences of DNA that are 
sufficient for RNA polymerase recognition, binding and transcription 
initiation. This portion of the promoter region is referred to as the 
promoter. In addition, the promoter region includes sequences that 
15 modulate this recognition, binding and transcription initiation activity of 
the RNA polymerase. These sequences may be cis acting or may be 
responsive to trans acting factors. Promoters, depending upon the nature 
of the regulation, may be constitutive or regulated. 

As used herein, the phrase "operatively linked" generally means the 
20 sequences or segments have been covalently joined into one piece of 
DNA, whether in single or double stranded form, whereby control or 
regulatory sequences on one segment control or permit expression or 
replication or other such control of other segments. The two segments 
are not necessarily contiguous. For gene expression a DNA sequence and 
25 a regulatory sequence(s) are connected in such a way to control or permit 
gene expression when the appropriate molecules, e.g., transcriptional 
activator proteins, are bound to the regulatory sequence(s). 

As used herein, production by recombinant means by using 
recombinant DNA methods means the use of the well known methods of 
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molecular biology for expressing proteins encoded by cloned DNA, 
including cloning expression of genes and methods, such as gene 
shuffling and phage display with screening for desired specificities. 

As used herein, a splice variant refers to a variant produced by 
5 differential processing of a primary transcript of genomic DNA that results 
in more than one type of mRNA. 

As used herein, a composition refers to any mixture of two or more 
products or compounds. It may be a solution, a suspension, liquid, 
powder, a paste, aqueous, non-aqueous or any combination thereof.* 
10 As used herein, a combination refers to any association between 

two or more items. 

As used herein, substantially identical to a product means 
sufficiently similar so that the property of interest is sufficiently 

* 

unchanged so that the substantially identical product can be used in place 

15 of the product. 

As used herein, the term "vector" refers to a nucleic acid molecule 
capable of transporting another nucleic acid to which it has been linked. 
One type of preferred vector is an episome, i.e., a nucleic acid capable of 
extra-chromosomal replication. Preferred vectors are those capable of 

20 autonomous replication and/or expression of nucleic acids to which they 
are linked. Vectors capable of directing the expression of genes to which 
they are operatively linked are referred to herein as "expression vectors". 
In general, expression vectors of utility in recombinant DNA techniques 
are often in the form of "plasmids" which refer generally to circular 

25 double stranded DNA loops which, in their vector form, are not bound to 
the chromosome. "Plasmid" and "vector" are used interchangeably as the 
plasmid is the most commonly used form of vector. Other such other 
forms of expression vectors that serve equivalent functions and that 
become known in the art can be used subsequently hereto. 
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As used herein, vector is also used interchangeable with "virus 
vector" or "viral vector". In this case, which will be clear from the 
context, the "vector" is not self-replicating. Viral vectors are engineered 
viruses that are operatively linked to exogenous genes to transfer (as 
5 vehicles or shuttles) the exogenous genes into cells. 

As used herein, transduction refers to the process of gene transfer 
and expression into mammalian and other cells mediated by viruses. 
Transfection refers to the process when mediated by plasmids. 

As used herein, "polymorphism" refers to the coexistence of more 
10 than one form of a gene or portion thereof. A portion of a gene of which 
there are at least two different forms, i.e., two different nucleotide 
sequences, is referred to as a "polymorphic region of a gene". A 
polymorphic region can be a single nucleotide, referred to as a single 
nucleotide polymorphism (SNP), the identity of which differs in different 
15 alleles. A polymorphic region can also be several nucleotides in length. 

As used herein, "polymorphic gene" refers to a gene having at least 
one polymorphic region. 

As used herein, "allele", which is used interchangeably herein with 
"allelic variant" refers to alternative forms of a gene or portions thereof. 
20 Alleles occupy the same locus or position on homologous chromosomes. 
When a subject has two identical alleles of a gene, the subject is said to 
be homozygous for the gene or allele. When a subject has two different 
alleles of a gene, the subject is said to be heterozygous for the gene. 
Alleles of a specific gene can differ from each other in a single nucleotide, 
25 or several nucleotides, and can include substitutions, deletions, and 

insertions of nucleotides. An allele of a gene can also be a form of a gene 
containing a mutation. 

As used herein, the term "gene" or "recombinant gene" refers to a 
nucleic acid molecule comprising an open reading frame and including at 



WO 03/023032 



PCT/IB02/03921 



-29- 

least one exon and (optionally) an intron sequence. A gene can be either 
RNA or DNA. Genes may include regions preceding and following the 
coding region (leader and trailer). 

As used herein, "intron" refers to a DNA sequence present in a 
5 given gene which is spliced out during mRNA maturation. 

As used herein, "nucleotide sequence complementary to the 
nucleotide sequence set forth in SEQ ID NO: x" refers to the nucleotide 
sequence of the complementary strand of a nucleic acid strand having 
SEQ ID NO: x. The term "complementary strand" is used herein 

10 interchangeably with the term "complement". The complement of a 
nucleic acid strand can be the complement of a coding strand or the 
complement of a non-coding strand. When referring to double stranded 
nucleic acids, the complement of a nucleic acid having SEQ ID NO: x 
refers to the complementary strand of the strand having SEQ ID NO: x or 

15 to any nucleic acid having the nucleotide sequence of the complementary 
strand of SEQ ID NO: x. When referring to a single stranded nucleic acid 
having the nucleotide sequence SEQ ID NO: x, the complement of this 
nucleic acid is a nucleic acid having a nucleotide sequence which is 
complementary to that of SEQ ID NO: x. 

20 As used herein, the term "coding sequence" refers to that portion 

of a gene that encodes an amino acid sequence of a protein. 

As used herein, the term "sense strand" refers to that strand of a 
double-stranded nucleic acid molecule that has the sequence of the 
mRNA that encodes the amino acid sequence encoded by the double- 

25 stranded nucleic acid molecule. 

As used herein, the term "antisense strand" refers to that strand of 
a double-stranded nucleic acid molecule that is the complement of the 
sequence of the mRNA that encodes the amino acid sequence encoded 
by the double-stranded nucleic acid molecule. 
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As used herein, an array refers to a collection of elements, such as 
nucleic acid molecules, containing three or more members. An 
addressable array is one in which the members of the array are 
identifiable, typically by position on a solid phase support or by virtue of 
5 an identifiable or detectable label, such as by color, fluorescence, 

electronic signal (i.e. radiofrequency (RF), microwave or other frequency 
that does not substantially alter the interaction of the molecules of 
interest), bar code or other symbology, chemical or other such label. 
Hence, in general the members of the array are immobilized to discrete 

10 identifiable loci on the surface of a solid phase or directly or indirectly 
linked to or otherwise associated with the identifiable label, such as 
affixed to a microsphere or other particulate support (herein referred to as 
beads) and suspended in solution or spread out on a surface. 

As used herein, a support (also referred to as a matrix support, a 

15 matrix, an insoluble support or solid support) refers to any solid or 

semisolid or insoluble support to which a molecule of interest, typically a 
biological molecule, organic molecule or biospecific ligand is linked or 
contacted. Such materials include any materials that are used as affinity 
matrices or supports for chemical and biological molecule syntheses and 

20 analyses, such as, but are not limited to: polystyrene, polycarbonate, 
polypropylene, nylon, glass, dextran, chitin, sand, pumice, agarose, 
polysaccharides, dendrimers, buckyballs, polyacrylamide, silicon, rubber, 
and other materials used as supports for solid phase syntheses, affinity 
separations and purifications, hybridization reactions, immunoassays and 

25 other such applications. The matrix herein can be particulate or can be 
in the form of a continuous surface, such as a microtiter dish or well, a 
glass slide, a silicon chip, a nitrocellulose sheet, nylon mesh, or other 
such materials. When particulate, typically the particles have at least one 
dimension in the 5-10 mm range or smaller. Such particles, referred 
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collectively herein as "beads", are often, but not necessarily, spherical. 
Such reference, however, does not constrain the geometry of the matrix, 
which may be any shape, including random shapes, needles, fibers, and 
elongated. Roughly spherical "beads", particularly microspheres that can 
5 be used in the liquid phase, are also contemplated. The "beads" may 
include additional components, such as magnetic or paramagnetic 
particles (see, e.g.,, Dyna beads (Dynal, Oslo, Norway)) for separation 
using magnets, as long as the additional components do not interfere with 
the methods and analyses herein. 

10 As used herein, matrix or support particles refers to matrix 

materials that are in the form of discrete particles. The particles have any 
shape and dimensions, but typically have at least one dimension that is 
100 mm or less, 50 mm or less, 10 mm or less, 1 mm or less, 100//m or 
less, 50 jjm or less and typically have a size that is 100 mm 3 or less, 50 

15 mm 3 or less, 10 mm 3 or less, and 1 mm 3 or less, 100 //m 3 or less and may 
be order of cubic microns. Such particles are collectively called "beads." 

As used herein, the abbreviations for any protective groups, amino 
acids and other compounds, are, unless indicated otherwise, in accord 
with their common usage, recognized abbreviations, or the IUPAC-IUB 

20 Commission on Biochemical Nomenclature (see, (1972) Biochem. 
//:942-944). 

B. High Throughput Process 

Provided herein are high throughput processes for the generation of 
and identification of proteins that exhibit desired phenotypes. The 
25 processes include methods that are particularly adapted for high 

throughput protocols, which require accurate methods for identifying 
modified proteins. 

A general directed evolution process includes the following steps: 
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1 . Generation of diversity at the nucleic acid level, on the gene 
to be 'evolved' 

2. phenotypic characterization of the gene variants generated; 

and 

5 3. identification of optimized gene variants. 

The processes provided herein effect these steps that can be 
performed in a high throughput format (see, FIGURES 1) that is optionally 
automated. A distinguishing characteristic of the processes provided 
herein, is that each candidate nucleic acid molecule is separately 

10 generated and screened, in an automated process at least some of the 
steps are performed without human intervention and are generally 
controlled by software. Most, if not all steps, are performed in 
addressable formats, such as at discrete locations in or on solid supports, 
such as microtiter plates or in other addressable formats, such as linked 

15 to coded supports. The supports can be electronically, physically, 

chemically or otherwise identifiable, such as by an identifiable symbology, 
including a bar code, or can be color coded. 

1 . Generation of Diversity using a semi-rational approach 

A semi-rational approach to creating diversity or evolving genes is 

20 provided herein. The goal is to create diversity but to decrease the 
number of molecules to be screened. By reducing the numbers, the 
molecules can be screened in high throughput format molecule-by- 
molecule (or groups thereof). 

Generation of diversity at the nucleic acid level, in principle, can be 

25 accomplished by a number of diverse technologies like mutagenesis 
(either site-directed or random), recombination, shuffling and de-no vo 
synthesis. These different technologies differ in the degree of diversity 
they generate as well as in the minimal length of the unitary change they 
can introduce (from single base to large domains). The outcome of step 1 
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is a collection of diverse, although highly related, molecules that 

constitutes what is called a 'library'. 

This step is crucial, since it provides the initial conditions for the 

entire process and is determinative of the outcome. The chances of 

5 finding an optimized gene version in a library is a function of the total 

diversity present in the library. In addition, the type of diversity 

introduced (such as, but not limited to, single point mutations, multiple 

point mutations, scarce small rearrangements, recombination of large 

domains, multiple shuffling) condition the outcome, particularly with 

10 respect to the generation of new variants compared to the original gene, 

and the probability that the new variants, not only exhibit the "evolved" 

function or property, but also work in their natural biological networks 

where they are expected to act by interacting, recognizing, and/or being 

recognized, by a large panoply of other proteins and other molecules. 

15 Rapid discovery of protein variants at the amino acid level by 

rational mutagenesis (aa-scan) 

A method, referred to herein as an amino-acid scan method for 
directed evolution, is provided herein for generating protein variants. This 
method can be performed on an entire protein or selected domains 

20 thereof, or can be used to identify benchmark sequences, such as 

functional domains, and, for example, recombine them as exchangeable 
units or restrict the diversity to limited or specific regions of the protein. 
Not only can this method be used with the processes provided herein, but 
it also has applications for any methods that use such variants or require 

25 generation of such variants, such as, but not limited to, searches for 
consensus sequences and homology regions that are used in functional 
genomics, functional proteomics; comparative modeling in protein 
crystallography and protein modeling; searches for natural diversity, 
(e.g., directed evolution methods in 6,171,820, 6,238,884, 6,174,673, 
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6,057,103, 6,001,574, 5,763,239, ); exon- or family-shuffling-based 
diversity {e.g., directed evolution using gene shuffling (see, e.g., U.S. 
Patent Nos. 6,096,548, 6,117,679, 6,165,793, 6,180,406, 6,132,970); 
the optimization of only the CDRs regions [e.g., directed evolution of 
5 antibodies see., e.g., U.S. Patent Nos. 5,723,323, 6,258,530, 5770,434, 
5,862,514) and other methods (see, e.g., U.S. Patent Nos. 5,837,500, 
5,571,698, 6,156,509). 

The amino-acid scanning-based method provided herein has 
advantages that prior methods do not have. For example, prior methods 

10 are based upon the underlying assumption that there are parts of the 
molecule {gene or protein) that are sufficiently adapted to perform their 
respective function, and further changes are not advantageous. Such 
methods do not look at total potential plasticity of a given molecule, but 
at the plasticity still permitted while keeping some basic functions in 

15 place. By choosing this route, however, additional potential variation is 
missed. The potential in the intrinsic plasticity of those regions that are 
presumed 'preserved' is lost. For instance, methods (e.g., those in U.S. 
Patent Nos. 6,171,820, 6,238,884, 6,174,673, 6,057,103, 6,001,574, 
5,763,239) that use natural diversity can miss the potential plasticity 

20 within those regions that are naturally 'conserved', i.e. there where there 
is no natural diversity. Methods that rely on exon- or family-shuffling- 
based diversity can miss the potential plasticity within regions contained 
in the shuffled fragments, i.e. within the fragments exchanged as a block. 
The method provided herein in contrast is sufficiently flexible to 

25 create mutants at a variety of levels, including at the single amino acid 
level; i.e. the method can generate mutants that differ from each other at 
a single amino and not at a larger block level. The challenge solved by 
the method herein is to generate diversity at the single amino acid level, 
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without moving too close to a pure 'random' approach, which results in 
an intractable number of mutants. 

The method provided herein is based on the premise that there are 
single amino acids or small blocks of sequence of amino acids that are 
5 either (1) directly involved in the activities that the methods 'evolve' 
(these amino acids would be at or close to the 'active sites' of the 
protein); or (2) directly involved in maintaining within the protein the intra- 
molecular environment that allows the active site(s) to stay active. 

Potential plasticity at the amino acid level can be exploited if amino 

10 acids or blocks of amino acids directly involved in the active sites for the 
activity to be evolved are known. Often they are not known. The 
problem that is solved herein, however, is how to exploit the potential 
plasticity at the amino acid level when nothing is known about the 
structure of the protein in question or about the position of its single or 

15 numerous active sites. 

The technology referred to herein as amino acid-scanning has been 
used to precisely identify those amino acids directly involved in the active 
sites of some enzymes and receptors (see, e.g., Becl-Sickinger et at. 
(1994) Eur. J. Biochem. 223:947-958; Gibbs et at. (1991)J. Biol. Chem. 

20 256:8923-8931; Matsushita eta/. (2000) J. Biol. Chem. 275:11044- 
1 1049) but has not been employed for directed evolution or for the 
generation of diversity. The amino acid scan as practiced in the prior art 
is used to produce a set of protein mutants, often within the region 
suspected to contain the active site(s), such that in each individual 

25 mutant a selected residue, such as Ala, replaces a different amino acid, 
Ala or other neutral amino acids generally, although not necessarily, is 
selected as a replacement amino acid since, except in instances in which 
the replaced amino acid is directly involved in an active site, it should 
have a neutral effect on the protein activity and not disturb the native 
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secondary structure of the protein. In instances in which the replaced 
amino acid is directly involved in an active site the activity of that site will 
be lost or altered. Amino-acid scanning, such as Ala-scanning, has been 
successfully applied for the identification of active sites in a number of 
5 proteins, and has been performed in computer-based rational drug design 
methods. Other amino acids, particularly amino acids that have a neutral 
effect, such as Glycine, can also be used. 

The amino acid scanning method is employed herein for the 
generation of the mutant proteins for screening for identification of sites 

10 or loci in a target protein or regions in a protein that alter a selected 

activity. In performing this method, the amino acids are each replaced, 
one-by-one along the full-length of the protein, or one-by-one. in pre- 
selected domains, such as domains that possess a desired activity or 
exhibit a particular function. Once sites of interest are identified other 

15 methods for generating diversity from the resulting molecules can be 
employed or the further steps of the method provided herein can be 
performed. 

The method includes the following steps: 

(1) Identification of the active site(s) on the full 

20 length protein sequence. In one embodiment a full-length amino acid- 
scan, typically, although not necessarily, an Ala-scan, or the identification 
and positioning of the active site(s) on proteins of either known or 
unknown function. For purposes herein, an active site is not necessarily 
the natural active site involved in the natural activity of a target protein, 

25 but those amino acids involved in the activities of the proteins under 

'directed evolution' with the purpose of either gain, improvement or loss 
of function. 

The whole process of the 'identification of the active site(s) on the 
full length protein sequence requires the following sub-steps: 
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a. Generation of a mutant library (on the 
gene to be evolved) in which each individual mutant contains a single 
mutation located at a different amino acid position and that includes a 
systematic replacement of the native amino acid by Ala or any other 

5 amino acid {always the same throughout the entire protein sequence); 

b. phenotypic characterization of the 
individual mutants, one-by-one and assessment of mutant protein activity; 

c. identification of those mutants that 

10 display an alteration, typically a decrease, in the selected protein activity, 
thus, indicating that amino acids directly involved in the active site(s) 
have been hit. The aa positions whose aa-scan mutations display an 
alteration, typically a loss or decrease, in activity are named HITS. 
The identification of the active site(s) (HITS) is thus, by this 

15 method, made in a completely unbiased manner. There are no 

assumptions about the specific structure of the protein in question nor 
any knowledge or assumptions about the active site(s). The results of the 
amino acid scan identify such sites. 

Once the active site(s) (the HITS) has(ve) been identified, those 

20 amino acids either at or surrounding the active sites, such as within 1 , 2, 
3, ... 10, 20 or any selected regions, as the unitary elements of 
exchange and generate diversity either at or around one of those sites or 
as a combined diversity at several sites at a time can be assessed. 
This process includes the following steps: 

25 a. Generation of a new mutant library (on the gene to be 

evolved) in which each individual mutant contains either single or multiple 
mutations located at (or surrounding) a specific active site (a HIT) position 
detected by the precedent aa-scan process. In the example these 
mutations include replacement, in each individual mutant, of the native 
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amino acid located either at (or surrounding) the HIT position by one of all 
other possible amino acids, such that, in the library, and at (or 
surrounding) each HIT position the native amino acid has been replaced 
by all possible amino acids. 
5 b. Identification of those mutants that display an 

increase in protein activity, thus indicating that a new sequence at or 
surrounding an active site has been identified with higher activity 
compared to the native sequence. These optimized sequences are named 
LEADS. 

10 The process can be repeated as many times as desired, in search 

for new combinations of optimal amino acids at (or surrounding) the 
different HIT sites. Each time, the process includes the steps of 
generating of a new mutant library (of the gene to be evolved) in which 
each individual mutant contains either single or multiple mutations located 

15 at (or surrounding) a specific active site (a HIT) position; phenotypic 

characterization of the individual mutants, one-by-one and assessment of 
mutant protein activity; and identification of those mutants that display an 
increase in protein activity, thus indicating that a new sequence at or 
surrounding an active site has been identified with higher activity 

20 compared to the native sequence. These optimized sequences are again 
named LEADS (second generation LEADS). 

2. Phenotypic characterization of the gene variants 
This step requires the expression of the gene variants in order to 
allow them to manifest their respective phenotypes. Gene expression can 

25 be accomplished by different means: in vitro, in reconstituted systems or 
in vivo in cellular systems, including bacterial and eukaryotic cells. For all 
exemplification purposes, reference is made to in vivo systems. Those 
of skill in the art can readily adapt these methods for in vitro systems, 
including those using biochemical assays. 
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This step is a crucial step for several reasons: 

(a) Expression system and protein processing. 
Depending on the system used {either bacteria or eukaryotic cells), 

as well as on the specific gene to be 'evolved', the variant proteins may 
5 or may not be appropriately processed, especially when post-translational 
modifications are involved, and therefore be able or not to elicit their 
potential activity. Consequently, the expression system {bacteria vs. 
eukaryotic cells) has to be carefully chosen. 

(b) Standardization of the expression system. 

■ 

10 The technologies available for gene transfer and expression into 

either bacteria or eukaryotic (let's consider mammalian) cells widely vary 
in their intrinsic efficiencies. While it is very easy to efficiently transfer 
and express genes in bacteria by chemical/physical methods 
(transformation), that is not the case for mammalian cells, where the 

15 transformation (here called transfection) process is inefficient and 

unreliable, especially when reproducibility and robustness are necessary in 
miniature, large number and small scale high throughput settings like 
those necessary to analyze gene variant libraries. Therefore, when 
transfection is used on mammalian cells, the specific activity measured 

20 for the individual variants in the library most probably does not accurately 
reflect the real specific activity of the molecules involved. As provided 
herein, transduction, the process of gene transfer and expression into 
mammalian and other cells mediated by viruses, overcomes the limitations 
of transfection. 

25 (3) Characterization. 

A distinction must be made between the 'expression' of the gene 
variants and their 'phenotypic characterization'. The expression system 
(either bacteria or mammalian cells) is only the vehicle to convert the 
gene variants into protein variants. The phenotypic characterization is 
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performed on the protein variants, and may have nothing to do, 
depending on the specific system under study, with the cellular system 
used to express the variants. The phenotypic characterization requires the 
use of specific assays (either biochemical (cell-free) or cell-based assays) 
5 in which the activity of the different cell mutants is challenged and 

assessed. In addition to the implications discussed below, these assays 
must be designed in such a way that they reflect the final environment in 
which the 'evolved' protein is expected to act. As an example, when 
optimizing an enzyme to be used in an artificial industrial setting, the 

10 assay should reproduce those conditions (temperature, pH, media 

composition...) of the real-life industrial reaction mixture, which may be 
relatively easy to do. When the final destination of the 'evolved' protein 
is a complex biological setting, such as the intracellular environment, the 
extra-cellular milieu (example: circulating proteins) or the structure of a 

15 virus, the necessary assay (s) may be quite difficult to setup. With a few 
exceptions, most of the work done so far on directed evolution has been 
made on simple enzymes for which all the necessary settings are 
relatively easy to implement. 

Methods for accurately titering viruses 

20 Much progress in gene therapy, genomics, biotechnology and, in 

general, biomedical sciences, depends on the ability to generate and 
analyze large numbers and small amounts of specific viruses. High 
throughput technologies are employed in disciplines such as functional 
genomics and gene therapy in which the use of viruses plays a key role 

25 for the efficient transfer and analysis of large collections of genes or 

libraries. Also, virus samples and biomedical samples containing viruses 
are routinely analyzed in thousands of hospitals, health centers, 
academic labs and biotech settings. 
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Furthermore in processes herein, accurate titration can be important 
in at least two steps in the process. After preparation of the viruses with 
the mutated variant, and prior to screening, it is necessary to know the 
concentration of titer of the viruses in the sample, so that results among 
5 the samples can be compared. The methods in this section designated 
Real Time Virus Titering (RTVT") and (TREE) are advantageously used. 

The methods in this section are also used in data analysis when 
measuring the output signal. As described below, output signal can be 
assessed by a Hill analysis or a second order polynomial or other 

10 algorithm that describes the interaction of biological molecules in complex 
system. In addition, where the output signal is actually the number of 
viral particles or ip produced, the methods in this section RTVT and TREE 
are advantageously used. 

Prior art virus titration methods (RCA, dRA...),for determining the 

15 amount of virus present in a biological sample, are based on the 

assessment of some kind of output signal, such as cytopathic effect, lysis 
or plaques and cell fusion focuses, induced in a reporter cell following a 
fixed time after infection with varying concentrations of the sample 
containing the virus. The lowest concentration of the sample at which no 

20 signal can be measured is taken as the titer of the virus in that sample. 
These approaches are known as "serial dilution" or "limiting dilution" 
methods. In limiting dilution methods, one single virus concentration, 
measured at a given time end-point gives rise to a single measurement of 
the output signal. These methods tend to be destructive in that 

25 assessment destroys the reaction so that only a single measurement can 
be taken on a sample. 
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Real Time Virus Titering (RTVT m ) 

When a virus infects a cell, the infected cell undergoes a number 
of changes that can be followed over time and quantified. Such changes 
are designated herein as the "output signal". The cell reports an output 
5 signal in response to the infection and, therefore, it is named here a 
reporter cell. One such output signal, is, for example, the expression of 
the genes carried by the virus (whether they are viral genes or 
exogenous genes (transgenes)). The output signal (for instance the level 
of expression of those genes) develops over time and depends, mainly, 

10 on two factors: i) the time point ("t") at which its level is measured after 
infection and ii) the amount of virus infecting the cell; i.e. the 
concentration of the virus preparation used to infect the cell ("s"). 
The output signal, at a given time point after infection, will be higher for 
higher concentrations of the virus infecting the reporter cells; and for any 

15 given concentration of virus, the output signal increases with time after 
infection until it generally reaches a plateau level. 

Real Time Virus Titering (RTVT™) published as International PCT 
application No. WO 01/186291, which is based on PCT/FR01/01 366 and 
EXAMPLES below) uses non-destructive methods for the assessment of 

20 output signal. Real Time Virus Titering is a viral titration method based on 
the kinetic analysis of the development of the output signal in 
virus-infected cells, tested at a single concentration of virus or biological 
sample. Instead of fixing the time point after infection and varying the 
concentration of the sample as is done in limiting dilution methods, in the 

25 Real Time Virus Titering RTVT W method, a fixed concentration of virus is 
used and the generation of a signal over time is assessed. Hence the 
signal is measured as a function of time at a single virus concentration. In 
this situation, a single virus sample (concentration), whose output signal 
is measured at a number of time points, can give rise to as many 
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measurements of the output signal as needed, and, eventually to a 
continuous, over time, assessment of the signal in real time. 

Real Time Virus Titering RTVf can be advantageously used in high 
throughput methods in which large numbers of biological samples are 
5 analyzed at the same time. This is the case, for instance, when titering 
viruses in a virus library. Limiting dilution methods rely on the output 
signal generated by a number of dilutions of each individual sample. If, for 
example, 10 dilutions (or experimental points) of each virus are used for a 
titration using a limiting dilution method, the analysis of a library 

10 containing 10,000 viruses require analysis of 10 5 (i.e., 10 x 10,000) 

experimental points. The Real Time Virus Titering RTVT" method requires 
only one dilution per sample, thereby requiring 10-fold fewer experimental 
points than a limiting dilution method. For a Real Time Virus Titering 
RTVT" titering system, the time (t£) necessary for the output signal to 

15 reach a reference value 05) is a direct function of the concentration of 

virus. Thus, t/3 can be used to directly determine the concentration of the 
virus. 

A limitation of the Real Time Virus Titering RTV"T limiting dilution 

titering method, however, is that not all the viruses (nor the genes carried 

20 by the viruses) generate a readily measured output signal that can be 

followed over time using non-destructive methods. 

Tagged Replication and expression enhancement 
(TREE) 

A method for titering designated Tagged Replication and Expression 
25 Enhancement Technology (TREE™) is provided herein. This system 
includes: i) a cell, ii) a reporter virus carrying a reporter gene, whose 
activity can be followed over time by a non- destructive method (i.e., 
fluorescence), iii) the virus (or virus library to be titered), herein referred 
to as the "titering virus". The elements are employed such that the virus 
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to be titered interferes with any output signal generated by the reporter 
virus, leading to either decrease or increase in the amount of that signal. 
The higher the amount of virus to be titered, the higher is the interference 
with the reporter virus and output signal. In the absence of virus to be 
5 titered, the kinetics of the output signal generated by the reporter virus 
are followed using the Real Time Virus Titering RTX/T" titering method. In 
the presence of increasing amounts of the virus to be titered the output 
signal takes longer (or shorter) to develop as a function of the amount of 

* 

virus to be titered. 

10 Using the TREE™ titering method, tfi, the time necessary for the 

output signal to reach a reference value Off) is a direct function of the 
concentration of virus and, therefore, t/3 can be used to determine the 
concentration of the virus to be titered. It is demonstrated herein (see the 
EXAMPLES) that when using the TREE system for titering, once an 

15 appropriate reference value Off) is determined for the output signal 
generated from the reporter virus, the time tjff is a function of the 
concentration of the virus being titered (see Example). Therefore, the 
concentration (titer) of the virus to be titered, can be assessed by 
assessing the change induced in Xfi by an aliquot of the virus to be titered. 

20 In a calibrated TREE titering assay, only one aliquot virus to be titered is 
needed to determine its titer, which is determined by measuring the shift 
in the t/3 of the system. The only condition is that the virus to be titered 
must "interfere" [i.e., increase or decrease) the output signal of the 
reporter virus. 

25 A calibration curve representing XJS vs. the amount of virus to be 

titered is obtained using aliquots of a reference batch of virus of known 
titer (previously determined using any titering procedure). The calibration 
curve can then be used to determine the amount of virus in a sample of 
unknown titer, based on the change caused by an aliquot of the sample 
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on the t/3 of the system and the corresponding titer read from the 
calibration curve. 

3. Identification of gene variants 
There are at least two considerations in this step: 
5 (a) Selection vs. screening. 

Depending on the specific protein involved, and under certain and 
very specific assay conditions, those variants that have been 'evolved' 
may elicit a selective advantage over the native version. This situation 
represents the most simple case: the cells (bacteria or mammalian) 

10 expressing the library of protein variants, as a pool or mixture, can be 

simply exposed to the selective conditions which by themselves will allow 
to put in evidence the best optimized variants. This situation is however 
very rare and difficult to achieve. It can be hardly believed that for any 
protein that one may want to optimize, a suitable 'selective' assay could 

15 be set up. For the vast majority of the cases, selection will not be 
possible. Therefore pools of molecules cannot be used, because the 
specific readouts of the assays could not be attributed to individual 
variants. When the simplistic selection approach is not possible, then two 
things are absolutely needed: (a) a 'one-by-one' approach, i.e. each 

20 individual variant must be physically separated from the others and its 
activity tested independently; (b) an accurate and quantitative analysis 
that can distinguish slight differences in activity among the different 
variants along a wide range of performance values. 

(b) Accurate quantitative analytics 

25 When selection is not possible, the optimized variants must be 

distinguished from the native variant otherwise. The different degrees of 
optimization among the different variants in the library should, in addition, 
be distinguished if those variants showing the highest optimization level 
are to be identified. A powerful quantitative analytical protocol is then 
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mandatory. These analytics should be able to attribute quantitative 
features (on the activity tested in the specific assay) to each of the 
variants tested and to rank them according to their individual 
performance. This requires, in addition, that each variant in the library is 
5 assayed individually; the use of pools or mixtures of molecules would 
hamper the ability to identify the right variants. 

For such analysis, the output signal can be assessed by a Hill 
analysis (see Examples and (published International PCT application No. 
WO 01/44809 based on PCT n° PCT/FROO/03503),.) or a second order 

10 polynomial (see, Examples and (Drittanti et al. (2000) Gene Then 7: 924- 
929)) or other algorithm that describes the interaction of biological 
molecules in complex system, such as the interaction between cells and 
biological agents. In addition, where the output signal is actually the 
number of viral particles or ip produced, the methods in this designated 

15 Real Time Virus Titering (RTVT W ) and Tagged Replication and expression 
enhancement (TREE™) are advantageously used (for a discussion of 
RTVT W , see, International PCT application No. PCT/FR01/01366 published 
as International PCT application No. WO 01/186291 and the EXAMPLES 

* 

below) or a refinement of that method provide herein and designated 
20 Tagged Replication and expression enhancement (TREE™) described 
above and in the examples. 
C. Practice of the process 

In one embodiment, the process provided herein includes the 

following steps. 

25 1 . Generation of diversity or source of existing diversity 

Generation of a plasmid library containing the genetic variants. The 
genetic variants are physically separated from each other. Any model 
such as, but not limited to, amino acid scanning, mutagenesis, or 
recombination may be used to generate the plasmid library. 
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2. Expression of the genetic variants 
Any method for expression of variants is contemplated. In 
particular the following alternatives are particularly suitable for high 
throughput performance. 
5 a. Expression in bacterial hosts 

The mutated forms of the nucleic acid are prepared or introduced 
into plasmids for expression in bacterial cells. The genetic variants are 
expressed from suitable bacterial cells, which are prepared by 
transformation of aliquots of the cells with each member of the plasmid 
10 library (each genetic variant continues to be physically separated from 
each other). 

b. Expression in eukaryotic host cells 
A virus library is generated from the plasmid library. The virus 
library, in which each different member is separately maintained, is 
15 prepared by: 

(1) Transfection of the plasmid library into 
appropriate virus-producer cells (viruses produced, each one carrying a 
different genetic variant present in the original plasmid library, are 
physically separated from each other); 

20 (2) Titration of the virus library (of each individual 

virus present in the library, separately). Titration is effected by any 
method, but generally by either a method designated Real Time Virus 
Titering (RTVT™) (see, International PCT application No. 
PCT/FR01/01366 published as International PCT application No. 

25 WO 01/186291 and the EXAMPLES below) or a refinement of that 
method provided herein and designated Tagged Replication and 
expression enhancement (TREE™) described above and in the examples; 
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(3) Standardization of the virus library to equal 
concentrations of all the individual viruses in the library (individual viruses 
continue to be physically separated from each other); 

(4) Expression of the genetic variants from 

5 appropriate mammalian cells by transduction with the virus library (each 
genetic variant continues to be physically separated from each other and 
each individual virus is handled separately from the others). 

3. Phenotypic characterization of the variant proteins. 
The variant proteins are expressed (from either plasmids in bacterial 
10 cells (step 2) or viruses in mammalian cells (step 4)) and their activity is 
assessed in one or more appropriate specific assays. The assays can be 
both types: biochemical (cell-free) assays and/or cell-based assays. The 
variant proteins in the library are physically separated from each other and 
their activity is individually assessed on a one-by-one basis. 
1 5 The assays can be performed in one of a variety of ways, 

including, but are not limited to: 

a. Using a single-point dilution for each individual variant 
protein, followed by a kinetic analysis (multiple time points) of the read- 
out by technologies like Tagged Replication and expression enhancement 

20 (TREE™), or any other appropriate technology 

b. Using serial dilutions of each individual variant protein, 
followed by, for example, the Hill-based analysis of the read-out by 
technologies or any other appropriate technology. Hill based analyses 
assess the interaction between cells and biological agents (see, published 

25 International PCT application No. WO 01/44809 based on PCT n° 
PCT/FROO/03503). 

The goal of these methods is to identify proteins having an evolved 
function or property. 
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Lead identification 

Based on the results obtained from the assays described above, 
each individual protein variant is individually tested for the parameters 
that assess the activity, property, function or structure of interest. 
5 Variants are ranked out according to their activity features. Those variant 
proteins best suited for the specificities of each individual project and 
system under study are then selected. The selected leads can be used for 
the desired purpose or further evolved or mutated to achieve desired 
activities. 

10 Typically, as for most directed evolution methods, the process is an 

iterative one, in which mutated variants are produced, screened, the best 
identified and then selected. The selected variants are then subjected to 
further evolution and the screening process repeated. This is repeated 
until the desired goal is achieved. 

15 This further evolution may employ the methods herein or any 

directed evolution method or combinations thereof. The methods for 
variant production will include the amino-acid scan method herein, which 
provides a rational approach to variant generation. Other rounds can 
include combinations of any other method for directed evolution known 

20 and/or combinations thereof. 

D. Directed evolution of a viral gene 

Recombinant viruses have been developed for use as gene therapy 
vectors. Gene therapy applications are hampered by the need for 
development of vectors with traits optimized for this application. The 

25 high throughput methods provided herein are ideally suited for 

development of such vectors. In addition to use for development of 
recombinant viral vectors for gene therapy, these methods can also be 
used to study and modify the viral vector backbone architecture, trans- 
complementing helper functions, where appropriate, regulatable and 
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tissue specific promoters and transgene and genomic sequence analyses. 
Recombinant AAV (rAAV) is a gene therapy vector that can serve as a 
model for application of the methods herein for these and other purposes. 
Adeno-associated virus (AAV) is a defective and non-pathogenic 
5 parvovirus that requires co-infection with either adenovirus or herpes 
virus, which provide helper functions, for its growth and multiplication. 
There is an extensive body of knowledge regarding AAV biology and 
genetics (see, e.g., Weitzman eta/. (1996) J. Virol. 70: 2240-2248 
(1996); Walker et al. (1 997) J. Virol. 77:2722-2730; Urabeefa/. (1999) 

10 J. Virol. 23:2682-2693; Davis eta/. (2000) J. Virol. 23:74:2936-2942; 
Yoon etal. (2001) J. Virol. 75:3230-3239; Deng et al. (1992) Anal 
Biochem 200:81-85; Drittanti etal. (2000) Gene Therapy 7:924-929; 
Srivastava etal. (1983) J. Virol. 45:555-564; Hermonat etal. (1984) J. 
Virol. 5/:329-339; Chejanovsky et al. (1989) Virology 773:120-128; 

15 Chejanovsky etal. (1990) J. Virol. 6*4:1764-1770; Owens etal. (1991) 
Virology 734:14-22; Owens etal. (1992) J. Virol. 66:1236-1240; 
Qicheng Yang etal. (1992) J. Virol. 66:6058-6069; Qicheng Yang eta/. 
(1993) J. Virol. 67:4442-4447; Owens et al. (1 993) J. Virol. 52:997- 
1005; Sirkka etal. (1994) J. Virol. 63:2947-2957; Ramesh etal. (1995) 

20 Biochem. Biophy. Res. Com. Vol 210 (3), 717-725; Sirkka (1995) J. 
Virol. 63:6787-6796; Sirkka et al. (1996) Biochem. Biophy. Res. Com. 
220:294-299; Ryan etal. (1996) J. Virol. 70:1542-1553; Weitzman etal. 
(1996) J. Virol. 70:2440-2448; Walker et al. (1997) J. Virol. 77:2722- 
2730; Walker etal. (1997) J. Virol. 77:6996-7004; Davis etal. (1999) J. 

25 Virol. 73:2084-2093; Urabe et al. (1999) J. Virol. 73:2682-2693; Gavin 
etal. (1999) J. Virol. 73:9433-9445; Davis etal. (2000) J. Virol. 
74:2936-2942; Pei Wu etal. (2000) J. Virol. 74:8635-8647; Alessandro 
Marcello etal. (2000) J. Virol. 74:9090-9098). AAV are members of the 
family Parvoviridae and are assigned to the genus Dependovlrus. 
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Members of this genus are small, non-enveloped, icosahedral with linear 
and single-stranded DNA genomes, and have been isolated from many 
species ranging from insects to humans. 

AAV can either remain latent after integration into host chromatin 
5 or replicate following infection. Without co-infection, AAV can enter host 
cells and preferentially integrate at a specific site on the q arm of 
chromosome 19 in the human genome. 

The AAV genome contains 4975 nucleotides and the coding 
sequence is flanked by two inverted terminal repeats (ITRs) on either side 

10 that are the only sequences in cis required for viral assembly and 

replication. The ITRs contain palindromic sequences, which form a hairpin 
secondary structure, containing the viral origins of replication. The ITRs 
are organized in three segments: the Rep binding site (RBS), the terminal 
resolution site (TRS), and a spacer region separating the RBS from the 

15 TRS. 

Regulation of AAV genes is complex and involves positive and 
negative regulation of viral transcription. For example, the regulatory 
- proteins Rep 78 and Rep 68 interact with viral promoters to establish a 
feedback loop (Beaton et al. (1989) J. Virol 63:4450-4454; Hermonat 

20 (1994) Cancer Lett 37:129-136). Expression from the p5 and p19 

promoters is negatively regulated in trans by these proteins. Rep 78 and 
68, which are required for this regulation, have bind to inverted terminal 
repeats (ITRs; Ashktorab et al. (1989) J. Virol. 63:3034-3039) in a site- 
and strand-specific manner, in vivo and in vitro. This binding to ITRs 

25 induces a cleavage at the TRS and permits the replication of the hairpin 
structure, thus, illustrating the Rep helicase and endonuclease activities 
(Im etal. (1990) Cell 67:447-457; and Walker et al. (1997) J. Virol. 
77:6996-7004), and the role of these non-structural proteins in the initial 
steps of DNA replication (Hermonat etal. (1984) J. Virol. 52:329-339). 
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Rep 52 and 40, the two minor forms of the Rep proteins, do not bind to 
ITRs and are dispensable for viral DNA replication and site-specific 
integration (Im et al. (1992) J. Virol. 66:1119-112834; Ni et al. (1994) J. 
Virol. 65:1128-1138. 
5 The genome (see, FIG. 4) is organized into two open reading 

frames (ORFs, designated left and right) that encode structural capsid 
proteins (Cap) and non-structural proteins (Rep). There are three 
promoters: p5 (from nucleotides 255 to 261: TATTTAA), p19 (from 
nucleotide 843 to 849: TATTTAA) and p40 (from nucleotides 1822 to 

10 1827: ATATAA). The right-side ORF (see FIG. 4) encodes three capsid 
structural proteins (Vp 1-3). These three proteins, which are encoded by 
overlapping DNA, result from differential splicing and the use of an 
unusual initiator codon (Cassinoti et al. (1988) Virology 767:176-184). 
Expression of the capsid genes is regulated by the p40 promoter. Capsid 

15 proteins VP1, VP2 and VP3 initiate from the p40 promoter. VP1 uses an 
alternate splice acceptor at nucleotide 2201; whereas VP2 and VP3 are 
derived from the same transcription unit, but VP2 use an ACG triplet as 
an initiation codon upstream from the start of VP3. On the left side of 
the genome, two promoters p5 and p19 direct expression of four 

20 regulatory proteins. The left flanking sequence also uses a differential 
splicing mechanism (Mendelson et al. (1986) J. Virol 60:823-832) to 
encode the Rep proteins, designated Rep 78, 68, 52 and 40 on the basis 
molecular weight. Rep 78 and 68 are translated from a transcript 
produced from the p5 promoter and are produced from the unspliced and 

25 spliced form, respectively, of the transcript. Rep 52 and 40 are the 
translation products of unspliced and spliced transcripts from the p19 
promoter. 

The rep protein is a adeno-associated virus protein involved in a 
number of biological processes necessary to AAV replication. The 
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production of the rRep proteins enables viral DNA to replicate, 
encapsulate and integrate (McCarty eta/. (1992) J. Virol 65:4050-4057; 
Horer eta/. (1995) J. Virol 65:5485-5496, Berns et at. (1996) Biology of 
Adeno-associated virus, in Adeno-associated virus (AAV) Vectors in Gene 
5 Therapy, K.I. Berns and C. Giraud, Springer (1996); and Chiorini eta/. 
(1996) The Roles of AAV Rep Proteins in gene Expression and Targeted 
Integration, from Adeno-associated virus (AAV) Vectors in Gene Therapy, 
K.I. Berns and C. Giraud, Springer (1996)). A rep protein with improved 
activity could lead to increased amounts of virus progeny thus allowing 

10 higher productivity of rAAV vectors. 

AAV and rAAV have many applications, including use as a gene 
transfer vector, for introducing heterologous nucleic acid into cells and for 
genetic therapy. Advances in the production of high-titer rAAV stocks to 
the transition to human clinical trials have been made, but improvement of 

15 rAAV production will be complemented with special attention to clinical 
applications of rAAV vectors as a successful gene therapy approach. 
Productivity of rAAV (i.e. the amount of vector particles that can be 
obtained per unitary manufacturing operation) is one of the rate limiting 
steps in the further development of rAAV as gene therapy vector. 

20 Methods for high throughput production and screening of rAAV have 
been developed (see, e.g., Drittanti et a/. (2000) Gene Therapy 7:924- 
929). Briefly, as with the other steps in methods provided herein, the 
plasmid preparation, transfection, virus productivity and titer and 
biological activity assessment are intended to be performed in an 

25 automatable high throughput format, such as in a 96 well (or other 
number or multiples thereof, such as 384, 1536 . . . 9600, 9984 . . .) 
formats. 

Since the Rep protein is involved in replication it can serve as a 
target for increasing viral production. Since it has a variety of functions 



WO 03/023032 



PCTYIB02/03921 



-54- 

and its role in replication is complex, it has heretofore been difficult to 
identify mutations that result in increase viral production. The methods 
herein, which rely on in vivo screening methods, permit optimization of its 
activities as assessed by increases in viral production. Provided herein 
5 are Rep proteins and viruses and viral vectors containing the mutated Rep 
proteins that provide such increase. The amino acid positions on the rep 
proteins that are relevant for rep proteins activities in terms of AAV or 
rAAV virus production are provided. Those amino acid positions are such 
that a change in the amino acid leads to a change in protein activity either 

•10 to lower activity or increase activity. As shown herein, the alanine or 
amino acid scan revealed the amino acid positions important for such 
activity (i.e. hits). Subsequent mutations produced by systematically 
replacing the amino acids at the hit positions with the remaining 18 amino 
acids produced so-called "leads" that have amino acid changes and result 

15 in higher virus production. In this particular example, the method used 
included the following specific steps. 
Amino acid scan 

In order to first identify those amino acid (aa) positions on the rep 
protein that are involved in rep protein activity, an Ala-scan was 

20 performed on the rep sequence. For this, each aa in the rep protein 

sequence was individually changed to Alanine. Each resulting mutant rep 
protein was then expressed and the amount of virus it could produced 
measured as indicated below. The relative activity of each individual 
mutant compared to the native protein is indicated in FIG 3A. HITS are 

25 those mutants that produce a decrease in the activity of the protein (in 
the example: all the mutants with activities below about 20 % of the 

native activity). 

In a second experimental round, which included a new set of 
mutations and phenotypic analysis, each amino acid position hit by the 



WO 03/023032 



PCT/IB02/03921 



-55- 

Ala-scan step, was mutated by amino acid replacement of the native 
amino acid by the remaining 18 amino acids, using site-directed 
mutagenesis. 

In both rounds, each mutant was individually designed, generated 
5 and processed separately, and optionally in parallel with the other 
mutants. Neither combinatorial generation of mutants nor mixtures 
thereof were used in any step of the method. 

A plasmid library was thus generated in which each plasmid contained 
a different mutant bearing a different amino acid at a different hit 
10 position. Again, each resulting mutant rep protein was then expressed 
and the amount of virus it could produce measured as indicated below. 
The relative activity of each individual mutant compared to the native 
protein is indicated in FIGURE 3B. LEADS are those mutants that lead to 
an increase in the activity of the protein (in the example: the ten mutants 
15 with activities higher, typically between 6 to 10 times, than the native 
activity). 

Expression of the genetic variants and phenotypic characterization. 
The rep protein acts as an intracellular protein through complex 
interaction with a molecular network composed by cellular proteins, DNA, 

20 AAV proteins and adenoviral proteins (note: some adenovirus proteins 
have to be present for the rep protein to work). The final outcome of the 
rep protein activity is the virus offspring composed by infectious rAAV 
particles. It can be expected that the activity of rep mutants would affect 
the titer of the rAAV virus coming out of the cells. 

25 As the phenotypic characterization of the rep variants can only be 

accomplished by assaying its activity from inside mammalian cells, a 
mammalian cell-based expression system as well as a mammalian cell- 
based assay was used. The individual rep protein variants were expressed 
in human 293 HEK cells, by transfection of the individual plasmids 
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constituting the diverse plasmid library. All necessary functions were 
provided as follows: 

(a) the cellular proteins present in the permissive specific 293 HEK 

cells; 

5 (b) the AAV necessary proteins and DNA were provided by co- 

transfection of the AAV cap gene as well as a rAAV plasmid vector 
providing the necessary signaling and substrate ITRs sequences; 

(c) the adenovirus (AV) proteins were provided by co-transfection 
with a plasmid expressing all the AV helper functions. 

10 A library of recombinant viruses with mutant rep encoding genes 

was generated. Each recombinant, upon introduction into a mammalian 
cell and expression resulted in production of rAAV infectious particles. 
The number of infectious particles produced by each recombinant was 
determined in order to assess the activity of the rep variant that had 

15 generated that amount of infectious particles. 

The number of infectious particles produced was determined in a 
cell-based assay in which the activity of a reporter gene, in the 
exemplified embodiment, the bacterial lacZ gene, or virus replication (Real 
time PCR) was performed to quantitatively assess the number of viruses. 

20 The limiting dilution (titer) for each virus preparation (each coming from a 
different rep variant) was determined by serial dilution of the viruses 
produced, followed by infection of appropriate cells (293 HEK or HeLa 
rep/cap 32 cells) with each dilution for each virus and then by 
measurement of the activity of the reporter gene for each dilution of each 

25 virus. Hill plots (NAUTSCAIST) as described herein (published as 
International PCT application No. WO 01/44809 based on PCT n° 
PCT/FROO/03503, Dec, 2000; see EXAMPLES) or a second order 
polynomial function (Drittanti et al. (2000) Gene Then 7: 924-929) was 
used to analyze the readout data and to calculate the virus titers. Briefly, 
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the titer was calculated from the second order polynomial function by 
non-linear regression fitting of the experimental data. The point where 
the polynomial curve reaches its minimum is considered to be the titer of 
the rAAV preparation. A computer program for calculation of titers has 
5 been developed (see Drittanti et a/. (2000) Gene Then 7: 924-929) to 
assess the minima. 

The TREE method described herein can be used to analyze the 
readout data and to calculate the virus titers. The results are shown in the 
EXAMPLE below. 

10 Comparison between results of full-length Hit position analysis 

reporter here and the literature 

The experiments identified a number of heretofore unknown 
mutation loci, which include the hits at positions: 4, 20, 22, 28, 32, 38, 
39, 54, 59, 124, 125, 127, 132, 140, 161, 163, 193, 196, 197, 221, 

15 228, 231, 234, 258, 260, 263, 264, 334, 335, 341, 342, 347, 350, 
354, 363, 364, 367, 370, 376, 381, 389, 407, 411, 414, 420, 421, 
422, 428, 429, 438, 440, 451, 460, 462, 484, 488, 495, 497, 498, 
499, 503, 511, 512, 516, 517 and 518 with reference to the amino 
acids in Rep78 and Rep 68. Rep 78 is encoded by nucleotides 321- 

20 2,186; Rep 68 is encoded by nucleotides 321-1906 and 2228-2252; Rep 
52 is encoded by nucleotides 993-2186, and Rep 40 is encoded by amino 
acids 993-1 906 and 2228-2252 of wildtype AAV. 

Also among these are mutations that may have multiple effects. 
Since the Rep coding region is quite complex, some of the mutations may 

25 have several effects. Amino acids 542, 598, 600 and 601, which are in 
the Rep 68 and 40 intron region, are also in the coding region of Rep 78 
and 52. Codon 630 is in the coding region of Rep 68 and 40 and non 
coding region of Rep 78 and 52. 
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Mutations at 10, 86, 101, 334 and 519 have been previously 
identified, and mutations, at loci 64, 74, 88, 175, 237, 250 and 429, but 
with different amino acid substitutions, have been previously reported. In 
all instances, however, the known mutations reportedly decrease the 
5 activity of Rep proteins. Among mutations described herein, are 
mutations that result in increases in the activity the Rep function as 
assessed by detecting increased AAV production. 
Lead identification. 

Based on the results obtained from the assays described herein (i.e. 

10 titer of virus produced by each rep variant), each individual rep variant 
was assigned a specific activity. Those variant proteins displaying the 
highest titers were selected as leads and are used to produce rAAV. 

In further steps, rAAV and Rep proteins that contain a plurality of 
mutations based on the hits (see Table in the EXAMPLES, listing the hits 

15 and lead sites), are produced to produce rAAV and Rep proteins that 
have activity that is further optimized. Examples of such proteins and 
AAV containing such proteins are described in the EXAMPLES. 

The rAAV rep mutants are used as expression vectors, which, for 
example, can be used transiently for the production of recombinant AAV 

20 stocks. Alternatively, the recombinant plasmids may be used to generate 
stable packaging cell lines. To create a stable producer cell line, the 
recombinant vectors expressing the AAV with mutant rep genes, for 
example, are cotransfected into host cells with a plasmid expressing the 
neomycin phosphotransferase gene (neor) by transfection methods well 

25 known to those skilled in the art, followed by selection for G418 
resistance. 

Also among the uses of rAAV, particularly the high titer stocks 
produced herein, is gene therapy for the purpose of transferring genetic 
information into appropriate host cells for the management and correction 
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of human diseases including inherited and acquired disorders such as 
cancer and AIDS. The rAAV can be administered to a patient at 
therapeutically effective doses. A therapeutically effective dose refers to 
that amount of the compound sufficient to result in amelioration of 
5 symptoms of disease. 

Gene therapy 

Toxicity and therapeutic efficacy of the rAAV can be determined by 
standard pharmaceutical procedures in cell cultures or experimental 
animals, e.g., for determining the LDS 50 (the dose lethal to 50% of the 

10 population) and the ED50 (the dose therapeutically effective in 50% of 
the population). The dose ratio between toxic and therapeutic effects is 
the therapeutic index and it can be expressed as the ratio LD 50 /ED 5O . 
Doses that exhibit large therapeutic indices are preferred. Doses that 
exhibit toxic side effects may be used; care should be taken to design a 

1 5 delivery system that targets rAAV to the site of treatment in order to 
minimize damage to untreated cells and reduce side effects. 

The data obtained from cell culture assays and animal studies can 
be used in formulating a range of dosage for use in humans. The dosage 
of such rAAV lies preferably within a range of circulating concentrations 

20 that include the ED 50 with little or no toxicity. The dosage may vary 
within this range depending upon the dosage form employed and the 
route of administration utilized. A therapeutically effective dose can be 
estimated initially from cell culture assays. A dose may be formulated in 
animal models to achieve a circulating plasma concentration range that 

25 includes the IC 50 (ie., the concentration of the test compound which 
achieves a half-maximal infection or a half- maximal inhibition) as 
determined in cell culture. Such information can be used to more 
accurately determine useful doses in humans. Levels in plasma may be 
measured, for example, by high performance liquid chromatography. 
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The following examples are included for illustrative purposes only 

and are not intended to limit the scope of the invention. The specific 

methods exemplified can be practiced with other species. The examples 

are intended to exemplify generic processes. 

5 EXAMPLE 1 

Titering or assessment of concentration by a method designated Real 
Time Vector Titering (RTVT™) 

This Example is based on the method described in International PCT 
application No. PCT/FR01/01 366, based on French application 

10 n° 0005852, filed 9 May 2000, and published as International PCT 
application No. WO 01/186291. This method assesses the 
titer or concentration of a biological agent (virus, gene transfer vector) in 
a sample, by measurement of the kinetics of change of a reporter 
parameter following the exposure of cells to the biological agent. 

15 As noted above, reporter parameters may include, but are not 

limited to,: gene / transgene expression related to the gene/transgene 
products, such as enzymatic activity, fluorescence, luminescence, antigen 
activity, binding to receptors or antibodies, and regulation of gene 
expression), differential gene expression, viral/vector progeny 

20 productivity, toxicity, cytotoxicity, cell proliferation and/or differentiation 
activity, anti-viral activity, morphogenetic activity, pathogenetic activity, 
therapeutic activity, tumor suppressor activity, oncogenetic activity, 
pharmacological activity. 

Serial dilution methods 

25 The assessment of the concentration or titer of biological agents 

■ 

using current approaches needs serial dilutions of the agent. 
Serial dilutions of the agent are applied to a cell-based reporter system, 
that elicits an output signal in response, to the exposure to the agent. The 
intensity of the signal is a function of the concentration of the agent. The 
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titer or concentration of the agent is determined as the highest dilution 
that still elicits a measurable response in the output. The higher the 

♦ 

number of dilutions tested, the higher the accuracy of the value obtained 
for the titer. 

5 This approach requires a set of serial dilutions for every biological 

agent whose titer needs to be determined. Thus, the application of this 

approach to the simultaneous titration of a large number of different 

biological agents is limited by the number of experimental points needed 

(example: for 30 biological samples: 20 serial dilutions x 30 biological 

10 agents: 600 experimental points). 

The approach in International PCT application No. 
PCT/FR0 1/0 1366 published as International PCT application 
No. WO 01/186291 

The intensity of the output signal (after exposure of reporter cells 
15 to the biological agent) is not only dependent on the concentration of the 
agent but also on the time after exposure. As time increases, the intensity 
of the signal increases. The kinetics of change of the intensity over time 
depends upon the concentration of the agent. Thus, lower concentrations 
of the agent will require longer times for the intensity to reach a given 
20 value that would be reached in shorter times after exposure to higher 
concentrations of the same agent. 

This approach (designated Real Time Virus Titering (RTVT") uses 
the following: a reference plot representing the relationship between the 
concentration of the agent and the time necessary for the intensity to 
25 reach a given threshold value is obtained using a reference preparation of 
biological agent, whose concentration or titer is known. This plot is then 
used to obtain the concentration of the biological agent under study by 
entering the time that a dilution of that agent needed for the intensity to 
reach the threshold value. 
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Using this approach, there is no need for serial dilutions of the 
biological agent(s) under study. Once the reference plot {t/3 vs c) is 
obtained, it can be used for the determination of the concentration or titer 
of any number of biological agents. Only one dilution of the biological 
5 agent under study is necessary, to obtain the corresponding value of tfi 
that is then used to obtain the concentration or titer using the reference 
plot. 

Thus, the application of this approach to the simultaneous titration 

* 

of large numbers of different biological agents is facilitated by the fact 
10 that only one dilution of each sample is needed (example: for 30 

biological samples: 1 dilution x 30 biological agents: 30 experimental 

points (compared to 600 needed with the current approach). 

This approach is specially suited for the high throughput assessment of 

concentration or titer of large numbers of biological agents. 
1 5 The system 

The system includes the following elements: 

a preparation of the biological agent (virus, gene transfer vector, 

protein,...) whose concentration or titer is unknown and has to be 

determined. 

20 a reporter system including culture of a cell line (or a mixture 

of cell lines) that reacts to the exposure to the biological 
agent by displaying a specific output signal, 
a master preparation of a reference biological agent, of 
known concentration or titer, that is able to generate the 

* 

25 output signal when the reporter cells are exposed to it. 

Practice of the method 

When the reporter cells are exposed to either the biological agent 
under study or the reference biological agent, an output signal is 
generated, that can be measured. 
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The intensity of the output signal is called /; the concentration of the 
biological agent used is called c, the time of exposure of the cells to the 
biological agent is called t. The intensity of the output signal (/) is a 

function of c and t : 
5 / increases as the concentration (c ) of the biological agent applied 

to the cells increases; 

/ increases as the time of exposure of the cells to the biological 
agent (t) increases. 

If the time t after exposure of the cells to the biological agent is 
10 kept constant, then, / will change as a direct function of c . 

If the concentration c of the biological agent is kept constant, 
then, / will be a direct function of t. 

fi is defined here as a threshold value of the intensity of the output 
signal, arbitrarily defined for every system under study. 
15 Xfi is defined here as the time necessary to reach the threshold fi. 

Use of P and tfi to determine the concentration or titer 
of a biological agent. 

The reporter cells are exposed to serial dilutions of a reference 
biological agent, whose concentration (or titer) is previously known. The 
20 intensity of the output signal (/) is measured at several time points (t) for 
every concentration (serial dilutions) c of the reference biological agent. 

/ is plotted vs t , and that, for every concentration c used of the 
reference biological agent. 

Using the plots obtained above, and for every concentration c of 
25 the reference biological agent, the time (xfi) necessary for the intensity of 
the output signal to reach a threshold valued is obtained. 

With the data obtained above, Xfi is plotted vs c. 

This plot represents the time necessary for the intensity of the 
output signal to- reach the threshold value jSasa function of the 
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concentration of the biological agent used. This is a standard plot and 
will be used to determine the unknown concentration of the biological 
agent under study by measuring the time that a given dilution of it needs 
to give an output signal whose intensity equals the threshold /?. 
5 The reporter cells are exposed to a dilution of the biological agent 

under study (whose concentration or titer is to be determined). The 
intensity of the output signal (/) is measured over time until it reaches the 
threshold value The time necessary for / to reach the value /? is 
recorded as tJ3. 

10 The t/3 value recorded above is entered into the standard plot 

obtained above and the corresponding c value is obtained. 

This c value represents the concentration or titer of the biological 

agent under study. 

Example of the Real Time Virus Titering RTVT W 
15 titering method 

Rat-2 cells were infected with serial dilutions of a reference 

preparation of a retroviral vector carrying the green fluorescent protein 

(GFP) gene (vector pSI-EGFPI see, Ropp eta/. (1995) Cytometry 27:309- 

317). At increasing times after infection, the level of expression of the 

20 transgene was determined (as the level of fluorescence due to the GFP 

gene) as the output signal. 

Table 3 represents the values obtained: 





Concentration 


Time after infection 


Output signal 




(1 = 10 6 particles/ml) 


(hrs) 


fluorescence 


25 


0.1 


16 


20.4 






24 


30.1 






40 


95.1 






48 


138.7 






64 


157.3 


30 


0.25 


16 


26.8 






24 


48.5 






40 


173.3 
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48 


228.2 




64 


191.7 


0.5 


16 


38.1 




24 


72 




40 


198.7 




48 


296.2 




64 


203.7 



The threshold value of £ = 100 was arbitrarily selected for this example. 
10 The time (t0) necessary for the output signal to reach the threshold P , for 
every concentration is shown in table 4. 



Table 4 



15 



Concentration 
(1 = 10 s particles/ml) 


(hrs) 


0.1 


42 


0.25 


31 



A plot of t/ff versus concentration for the reference virus shows that 
the concentration and t0 exhibit a clearly defined relationship, that allows 
20 for the calculation of the concentration (c) of a sample, if the 
corresponding tft of that sample is known. 

EXAMPLE 2 

Tagged Replication and Expression Enhancement (TREE) for titering 
As discussed above, TREE is a method for titering and 

25 standardization of preparations of viruses, vectors, antibodies, libraries, 
proteins, genes and any other moiety that is detectable based upon a 
output signal, such as fluorescence. The TREE method is an improvement 
of the Real Time Virus Titering (RTVT) method (see, International PCT 
application No. PCT/FR01/01366 published as International PCT 

30 application No. WO 01/186291). It is performed with a reporter moiety, 
such as a reporter virus (with a known titer) and the test sample (with 
unknown titer). The reporter, such as a reporter virus, has a readily 
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detectable output signal that can be measured as a function of time. The 
effect of the moiety, such as a virus, of unknown titer is assessed. The 
moiety whose titer is assessed either increases or decreases the output 
signal as a function of time. This change in signal is used to assess the 
5 amount or concentration of the moiety of unknown concentration, and 
hence its titer. 

The method is exemplified herein using an AAV system for the 
determination of the titer of an AAV vector and an AAV-reporter vector 
as a competitor and wild type Adenovirus as helper virus. One of skill in 
10 the art readily can adapt the method to other systems, including other 

viruses, and other moieties for which a reporter system can be developed. 

u 

Other such moieties include, but are not limited to, viral vectors, 
plasmids, libraries, proteins, antibodies, vaccines, genes, and nucleic acid 
molecules. 

1 5 Materials and Methods 

1 . Cells and Viruses 
HeLa rep-cap32 cells, a HeLa derived cell line (kindly provided by P. 
Moullier, Laboratoire de Therapie Genique, CHU, Nantes; see, Salvetti et 
a/. (1998) Hum Gene Ther 20:695-706; Chadeuf eta/. (2000) J Gene 

20 Med 2:260-268) was grown in DMEM with 10% fetal calf serum. These 
cells were plated 24 h before infection at a density of 1 x 10 4 cells in 
single well of 96-well plates. rAAV-LacZ (10 10 ip/ml), rAAV-GFP (10 9 
ip/ml) vectors and Human Adenovirus type 5 (Ad5) (10 11 pfu/ml) were 
from CHU, Nantes. 

25 Hela rep-cap32 cells had been produced by cotransfecting plasmid 

pspRC, which harbors the AAV rep-cap genome with the ITRs deleted (bp 
190 to 4484 of wild-type AAV), with plasmid PGK-Neo, conferring 
resistance to G418 on Hela cells (see, Chadeuf eta/. (2000) J. Gene 
Med. 2:260-268 and Salvetti eta/. (1998( Hum Gene Ther. S:695-706). 
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Hela rep-cap 32 cells are a packaging line that harbor one copy of the 
genome with the ITRs deleted (see, also Tessier et al. (2001) J. Virol. 
75:375-383). 

■ 

Plasmid pspRC contains the AAV genome (positions 190-4,484 bp) 
5 with the ITRs deleted and was obtained by excising the rep-cap fragment 
(Xba\ fragment) from the well-known vector psub201 (Samulski et al. 
(1987) J Wra/ 67:3096-3101; also called pSSV9) by Xba\ digestion and 
inserting it into the Xba\ site of plasmid pSP72 (Promega). Plasmid 
psub201 (see, e.g., U.S. Patent No. 5,753,500) is a modified full-length 

10 AAV type 2 genomic clone contains all of the AAV type 2 wild-type 
coding regions and cis acting terminal repeats. 

2. Infection and measurement 
Four serial dilutions of a rAAV-LacZ (0.01, 0.0075, 0.005 and 
0.0025 pi, see Table 2 below, designated samples 1-4, respectively) were 

15 made and used for co-infection of HeLa rep-cap32 cells together with 8 
different Ad5 multiplicity of infection (MOI; from 0.1 to 100/cell) and with 
10* ml (10 6 infectious particles (ip)) or 10 4 ml (10 s ip) rAAV-GFP viral 
vector. All the samples were done in triplicate. After infection, the plates 
were read at different times, from 34.5 h to 80 h (every 30 minutes). 

20 rAAV-GFP is an SSV9-derived vector; SSV9 is a clone containing 

the entire adeno-associated virus (AAV) genome inserted into the Pvull 
site of plasmid pEMBL (see, Du et al. (1996) Gene Ther 3:254-261). The 
rAAV-GFP and rAAV-LacZ plasmids are SSV9 with a GFP or LacZ gene 
under control of the cytomegalovirus (CMV) immediate-early promoter. 

25 All the samples were done in triplicate. 
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3. Process 

Figure 2 shows the overall procedure in 96 well format. Cells were 
plated 24 h before infection. Co-infection of rAAV-GFP with serial 
dilutions of rAAV-LacZ together with Ad5 (different MOI), were done. 
5 Then the plates were read at different times using the Analyst AD&HT 
micro plate-reader (LJL BioSystems). 

4. Analysis 

For this kinetic technique, Fluorescence Intensity (Fl) of the infected 
cells is measured as a function of the time. Serial dilutions of the AAV- 

10 competitor vector AAV-lacZ vector, which decreases the fluorescence 
signal, are performed. For this example, fluorescence was measured for 
AAV-GFP with 10 6 ip and 10 s ip and then 10 6 ip of the AAV-GFP reporter 
virus in the serial dilutions of the competitor virus, AAV-lacZ vector in a 
96-well format (samples 1-4, see Table below). Measurements were 

15 taken of each well and curves of Fl (of the GFP) versus time (hrs) were 
obtained (see Figure 2B). 

An arbitrary one value for Fl (see Fig. 2B, 6 x 10 6 Fl units), 
typically, though not necessarily, near the greatest separation among the 
curves so that the numbers are readily discernable, was selected. The 

20 point at which each of the curves intersect this value is beta time {t/3) for 
that combination of amounts of reporter plus dilution of the virus of 
unknown titer. t&, taken from the Fl vs. Time (hrs) curves, for each 
sample containing a dilution of the unknown plus 10 6 ip of the reporter 
virus is set forth in column 3 of Table 2 below. 

25 To determine the titer of the test virus, the X0 for the AAV GFP 

(reporter virus) is plotted versus quantity of ip (i.e a straight line between 
the tfi for the 10 6 ip and the 10 5 ip) (Fig. 2C). For any tp of the unknown 
virus, the quantity of ip can be determined from this curve. The beta 
time (tfi) of each sample (in this case for the different dilutions of rAAV 
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10 



15 



20 



25 



LacZ mixed with 10 6 infectious particles of rAAV-GFP) is determined, and 
then the residual number of infectious particles of rAAV-GFP for each 
sample. The difference between 10 6 ip of rAAV-GFP put in each sample 
and the number of ip detected by fluorescence in the same well is the 
actual quantity of rAAV-GFP competed (consumed) by the unknown 
rAAV (in this case rAAV-LacZ). This number is determined for each 
dilution. The quantity rAAV-GFP consumed is the same quantity of 
unknown rAAV in the sample. This quantity is present in one volume of 
unknown rAAV, which in this example is 1 ml. Based upon this, the 
infectious titer of the unknown rAAV is determined. The results are 
shown in Table 2. 

TABLE 2 

AAV LacZ titration by TREE titration 



Sample 


volume U/l) 


tfi (hrs) 


Residual 
AAV-GFP 


Consumed 
AAV-GFP 


AAV LacZ 
Concentration 
<i.p./10- 2 //I) 


1 


10.0 x 10' 3 


66.5 


5.56 x 10 5 


4.44 x 10 5 


4.44 x 10 s 


2 


7.5 x 10* 3 


66.5 


5.56 x 10 s 


4.44 x 10 3 


5.93 x 10 s 


3 


5.0 x 10- 3 


65 


7.06 x 10 s 


2.94 x 10 s 


5.88 x 10 s 


4 


2.5 x 10' 3 


64 


.8.06 x 10 5 


1.94 x 10 s 


7.76 x 10 s 



30 



\p/jj\ = 6.01 x 10 5 ip/0.01 //I). The standard deviation was 1.37 x 10 10 
ip/ml with an error of ±23%. 

EXAMPLE 3 

Hill Analysis of the screening assay output 

It is important to have reliable methods for screening and/or 
evaluating the performance of a set of biological agents, such as a library 
of viral or non-viral recombinant vectors, vaccines, recombinant proteins 
and antibodies, in a complex biological system, such as living target cells 
When developing such agents, for example gene therapy vectors and 
other agents for therapeutic use, it is necessary to be able to evaluate 
and compare performance among candidates. 
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The progress of gene transfer into gene therapy depends upon the 
capacity to develop gene transfer vectors into therapeutic drugs. Clinically 
relevant vectors need to be efficient and safe, in reaching and infecting 
target cells and in ensuring a persistent level of expression of the 
5 therapeutic gene with a minimum of adverse effects. The availability of 
standardized quantitative methods, suitable for an accurate and objective 
assessment of titer, performance and safety, is necessary for the 
pharmaceutical development of gene vectors as drugs. 

Any method for assessment is contemplated herein as long as it is 

10 adapted for use in a high throughput format. Of particular interest is the 
Hill equation based method of international PCT application No. 
WO 01/44809 (International PCT application No. PCT/FR00/03503, based 
on French application FR 9915884). 

Two widely used parameters that provide quantitative information 

15 about the potential performance of a gene transfer vector preparation are 
the titer of physical particles and the titer of infectious particles. Vector 
preparations. with high titer of infectious particles and low physical 
particles/infectious particles ratio are considered to be of higher quality. 

The titer in physical particles(pp) (see, e.g., Mittereder et at. (1996) 

20 J. Virol. 70:7498-7509; Atkinson eta/. (1998) Nucl. Acids Res. 26:2821- 
2823; Kechli eta/. (1998) Hum. Gene Ther. 5:587-590; and Nelson eta/. 
(1998) Hum. Gene Ther. S:2401-2405), which represents the total 
number of vector particles, is usually evaluated from the vector content 
by detecting the nucleic acid contents (nucleic acids hybridization and 

25 OD 260 respectively for AAV and AdV), detecting viral protein content (for 
example, reverse transcriptase (RT) activity and p24 content for MLV and 

HIV, respectively). 

Among the physical particles (pp), there are particles potentially 
active in performing transduction (ip, infectious particles), as well as 



WO 03/023032 



PCT/IB02/03921 



-71- 

particles that are inactive (nip, non-infectious particles) (Ruffing et al. 
(1994) J. Gen. Virol. 75:3385-3392; Kechli eta/. (1998) Hum. Gene Ther. 
5:587-590.). The pp and the ip/nip ratio, are features of the packaging 
system, the manufacturing process and the vector itself. 
5 The infectious particles (ip) (infectious units, transducing units, 

etc.) are evidenced by the changes observed in the infected cells (vector 
DNA replication, provirus integration, cell lysis, transgene expression and 
other observable parameters). Infectious particles (ip) measures the 
number of particles effective in performing a process whose output is 

10 being measured; not all particles participate or are capable of participating 
in all processes. 

The precise assessment of ip is not straightforward. Existing 
methods are mainly based on serial dilution experiments followed by 
either linear extrapolation or asymptote approximates. The titer of 

15 infectious particles {ip: infectious unity, transduction unity) (see, e.g., 
Mittereder et al. (1 996) J. Virol. 70:7498-7509; Weitzman et al. (1996) 
J. Virol. 70:1845-1854; Salvetti et al. (1998) Hum. Gene Ther. 5:695- 
706) is evaluated by the studying observed changes in infected cells, 
such as viral replication, provirus integration, cellular lysis and transgene 

20 expression, using methods based on serial dilutions, followed either by a 
linear extrapolation or an asymptotic approximation. Thus, ip measures 
the number of active particles in the measured process; it includes 
physical particle (pp) and inactive particles {nip or non-infectious 
particles) . 

25 In order to resolve the problem of the titer determination and the 

comparison of different recombinant viruses used in gene therapy, the 
variation of the particles/infectious power ratio has been used (see, e.g., 
Atkinson et al. (1998) Nucl. Acids Res. 25:2821-2823; and International 
PCT application No. WO 99/1 1764, which describe a method that uses 
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step of amplification viral genetic material in a host cellular line, 
preparation of vectors of unknown titer obtained by serial dilution and an 
internal check of known titer). In particular, the method uses cells 
infected with a viral preparation in the different wells of a microtiter plate, 
5 viral genome replication in the host cells, nucleic acid hybridization and 
determination of the relative amount of replicated viral nucleic acid in 
each well. 

All of these methods measure the physical particles (pp) titer and/or 
measure infectious particles Up) titer in order to evaluate a gene transfer 

10 vector. A high quality vector preparation is one with an high titer of 
infectious particles and a low pplip ratio. These parameters provide 
quantitative information on the performance of a gene transfer vector. 
Because of the inaccuracy of the procedures used for assessing pp and 
especially ip, these parameters are not informative enough to precisely 

15 define the features of a gene therapy vector nor those of a particular 

preparation thereof. The actual procedures used for pp and./p eyaluation 
change with the vector type, are not very reproducible nor exact, so 
these parameters do not contain enough information to allow a very fine 
definition of vector characteristics and performances. 

20 Hill equation-based analyses 

In this method complex biological processes, including those 
involving the response of cells (in vitro and in vivo) to biological agents, 
such as, for example, cells, viruses, vaccines, viral and non-viral gene 
vectors, antibodies, antigens, proteins in general and plasmids, are 

25 characterized using the formal analysis first introduced by Hill (see, Hill 
(1910) J. Physiol. 40:4P; Hill (1913) /. Biochem.J. 7:471-480). 
International PCT application No. WO 01/44809 (based on- 
PCT/FROO/03503, priority claimed to French application FR 9915884) 
describes the use of the Hill equation (see, Hill (1910) J. Physiol. 40:4P; 
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Hill (1913) A Biochem.J. 7:471-480; see, International PCT application 
No. PCT/FR00/03503) for analysis and characterization of the biological 
and/or pharmacological activity of biological agents (viruses, vectors or 
cells) on biological assay systems in vitro (cell-based) or in vivo. 
5 A number of useful parameters, derived from the Hill equation, are 

scored and used to quantify relevant features of the biological agent, of 
the cells, as well as of the biological process or reaction involved. 

In particular, methods for calculation and analysis of the 
parameters of biological and pharmacological activity of native, 

10 attenuated, modified or recombinant viruses, vaccines, recombinant viral 
and non-viral gene transfer vectors, cells, antibodies and protein factors 
in in vitro (cell-based) or in vivo assays are described. This method is 
adapted for high throughput processes and is sufficiently accurate to 
allow a very fine definition of vector characteristics and performances. 

15 International PCT application No. WO 01/44809 provides, a 

standard process for evaluating the interaction between any biological 
agent, such as a gene therapy vector, with a complex biological system 
(living target cells). It provides a screening process for a pool of complex 
biological agents, in order to select test agents that have a desired 

20 property, activity, structure or whatever is being sought. 

Different biological agents and assay systems (cells) are compared 
and ranked out on the basis of their performance, assessed through the 
Hill parameters. Thus, the accurate analysis and comparison of the 
biological response of complex assay systems (in vitro and in vivo) to 

25 complex biological agents is achieved experimentally. The Hill-based 

analysis (rf,/c,r,er,/7,e) is used for a variety of purposes, including, but not 
limited to: 
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i) validation and optimization of the manufacturing processes used 

to obtain the biological agents; 

ii) development and optimization of the components of the 
biological agents (proteins, genomes, genetic units); 

5 Hi) development and optimization of assays and analytical tests for 

the characterization of the biological agents. 
The method includes the steps of: 

(a) preparation, for each biological agent, of a sample scale, 
obtained by a serial dilution of the biological agent at a R1 concentration, 
10 (b) incubation of each sample of the dilution scale obtained in (a), 

with the target cells at a constant concentration R2, 

(c) determination of the P product from the reaction R1 + R2, at a 
t moment, in each the sample; and 

(d) realization of a theoretical curve H from the experimental points 
15 R1 and P, for each biological agent by iterative approximation of 

parameters of the reaction R1 + R2 -* P, at the t moment, in accordance 
with this equation: 

P - Pmax (/rR1) r /(/c + </rR1) r ) r=1,...,n (2) 

in which: 

20 R1 represents the biological agent concentration in a sample 

from the scale; 

R2 is concentration of target cells (in vitro or in vivo) 

P (output) represents the product from the reaction R1 + 

R2 at a t moment; 
25 P max represents the reaction maximal capacity; 

k represents, at a constant R2 concentration, the resistance 

« -» 

of the biological system for responding to the biological agent (resistance 
constant R2); 
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r represents a coefficient that depends on R1 and 
corresponds to the Hill coefficient; and 

77 represents the intrinsic power of the R1 biological agent to 
induce a response in the biological system (P production at the t 

5 moment), and 

(e) sorting the k and rr values obtained in (d) for each 
biological agent and the biological agent, and then ranking according to 
the values thereof. 

Using the parameters (n t K,T,6 t €,n) the activity of a biological agent 
10 on a complex biological system, as well as its intrinsic features can be 
fully characterized and compared. In addition, different biological 
systems* either in vitro cell-based) or in vivo can be compared. 

Hill Equation 
The Hill equation: 
15 r = n 

P - Z P max . R1 r / (K + (R1) r ) R2 constant (1) 

r =1 

where R1, P, P max and K represent, respectively, the concentration 
of the reagent R1, the concentration of the product, the maximal capacity 

20 of the reaction and the 'affinity' constant between R1 and R2. The Hill's 
coefficient (r) is a function of R1 . The coefficient r is equal to 1 when 
independent non-interactive binding sites are involved between R1 and 
R2, such as in reactions that follow kinetics described by Michaelis- 
Menten; and r varies from 1 to n for systems where the sites involved in 

25 the interaction between the Rl and R2 are not independent from each 

other, and the affinity for R1 at any R2 binding site varies as a function of 
either i) the degree of occupancy of other R2 sites; ii) the concentration 
of R1 itself or iii) the concentration of other (positive or negative) 
regulators. 
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The Hill equation, thus, is a general formalization that describes the 
interaction reaction between molecules. It expresses the amount of 
product formed as a function of the concentration of the reagents and of 
the affinity constant of the system. Originally developed for the study of 
5 the dissociation between haemoglobin and oxygen, the Hill equation 
covers the formal Michaelis-Menten analysis of enzyme kinetics, the 
analysis of ligand-receptor binding and of the allosteric protein systems. 

According to Hill, for a simple reaction like 

R1 + R2 -* P 
10 K 

where the affinity between R1 and R2 changes with concentration 

of each, the Hill equation describes the accumulation of the product P as 

a function of the concentration of one of the reagents (R1) and of the 

intrinsic properties (K) of the system. 

15 This equation can be applied to complex biological systems. For 

example, the response of the cells to infection (P), can be analyzed by 
applying an Hill-type equation. The amount of cells growing in vitro (R2) 
are infected with increasing concentrations of recombinant viruses (R1), 
and (P) is monitored. A Hill equation is iteratively fitted to the 

20 experimental data. 

For analyses of viral output as exemplified herein, 

virus + cell -* transduced cell output (viral genome replication), 
Equation (1) is specifically reformulated as: 

25 P = P max (/7R1) r / (k + {/rR1) r ) r=1,...,n (2) 

where P, P max , R1, rr, r and k , as described above, represent, 
respectively, the output signal (P) (the level of viral gene expression, or 
the level of virus replication), the maximal output signal (Pmax), the initial 
concentration (R1) of infectious viral particles (those susceptible to trigger 

30 the process leading to P), the potency of the vector (tt; a factor that 
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affects the concentration of the vector (R1) by its specific strength or 
activity, the Hill's coefficient (r) and the constant of resistance of the 
reaction or process {k). 

The constant of resistance h 
5 The concept of k is analogous to those of dissociation, kinetics, 

equilibrium or affinity constants concepts for simple chemical and 
biological reactions, k is a feature of the process (reaction) and of the 
biological system tested (cell type), k is a key parameter for the 
characterization of the assay system and the assessment of its 

10 performance as a test for the reaction under study. 

k measures the internal resistance offered by the process or 
reaction triggered by the biological agent, to proceed to P. k is specific to 
a particular process or reaction tested. In addition k is specific to the 
particular biological system tested. Different cell lines and types will 

15 display different k for the same reaction. Moreover, factors affecting the 
performance of a cell to accomplish the reaction (like contaminants, toxic 
agents, etc.) affect k in that cell. 

Variations in k affect equation (2) by shifting the curve to the right 
or to the left, according to whether the value of k increases or decreases, 

20 respectively. All curves differing only in k are parallel each other. 

k finds its direct and practical application in i) assay development and 
validation and ii) assessment of the susceptibility or sensitivity of different 
cell types or tissues to undertake the reaction under study and to be 
affected by it. 

25 The potency n 

it measures the intrinsic potency of the biological agent to 
accomplish P against the resistance (k) offered by the reaction process. 
For every infectious virus particle (R1) added to the assay, the actual 
activity of the virus added is given by /7*R1 . In order to report an output P, 



WO 03/023032 



PCT/IB02/03921 



-78- 

the potency rr has to push forward the reaction inside the cell against k. 
rr is specific to the particular biological agent for the reaction under study. 
rr is a feature of the biological agent. 

Different versions or variants of the biological agent will display 
5 different rr for the same reaction. Thus, mutations, conformational 
changes or other modifications on the biological agent are expected to 
change its rr for a given reaction process. 

The concept of rr is analogous to that of chemical activity by 
opposite to concentration for simple compounds, rr is a correction factor 
10 that affects the concentration (R1) of the biological agent to indicate its 
actual strength or activity for a given reaction process. 

Variations in rr affect equation (2) by shifting the curve to the right 
or to the left, according to whether the value of rr decreases or increases, 
respectively. Curves differing only in rr are not parallel each other. The 
15 slope of the curve given by equation (2) increases as rr increases. 

rr is a key parameter for the characterization of the biological agent 
and the assessment of its performance to accomplish the reaction under 
study, rr finds its direct and practical application in i) biological agent 
optimization and development as it allows to compare the relative potency 
20 of variants of the agent. 

rr is a valuable tool in the field of vaccine, gene transfer vector and 
antibody development, for the comparison between two or more different 
agents or different versions of the same agent, for performance. Two 
agents, for instance, may elicit equivalent potencies for gene transfer, 
25 while their potencies for immunogenicity be different. The use of rr, a 

quantitative and accurate parameter for assessing potency, will allow for 
ranking the candidates according to their potency {i.e., for gene transfer, 
gene expression, immunogenicity and other such properties and activities) 
and to make rational decisions about the relative value of the agent leads. 
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The efficiency e 

e measures the maximal global efficiency of the reaction process 

■ • 

when a biological agent characterized by a given n value interacts with a 
biological system characterized by a given k. e is specific to the particular 
5 couple biological agent {rr) / biological system M for the reaction under 
study. € is a feature of the global reaction process and intervening 
reagents. Changes in either rr , k, or both, will lead to changes in e. 

The efficiency of the reaction process described by equation (2) is 
given by the increase in the output P that can be obtained by increasing 

10 the input R1 . Thus, the first derivative of P with respect to R1, or the 

slope of the curve described by equation (2), gives the global efficiency of 
the reaction at every R1 input. The maximal global efficiency, or e, is 
given directly by either the slope at the inflection point of the curve 
described by equation (2) or by the maximum of its derivative JP/JR1 . 

15 The slope of the curve given by equation (2) and the maximum of <JP/c5R1 
increase as € increases. 

6 is a key parameter for the characterization of the efficiency global 
process, considering the assay conditions and reagents all together. It is 
therefore useful for assay optimization once rr and k have been fixed and 

20 to detect changes in rr when k is kept constant or, inversely, changes in k 
while rr is kept unchanged. 

The heterogeneity index n 
n measures the internal heterogeneity of the reaction process under 
study. Complex processes include a huge chain of individual and causal 

25 events inside a multidimensional network of interrelated and interregulated 
biological reactions. Thus, the constant of resistance (#c) for the particular 
reaction process under study is a macroscopic indicator of the global 
resistance of that process (k = a1/c1 +a2/c2 + ...an/cn/n). If the 
contribution of the individual microscopic constants of resistance {a1/c1, 
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a2/c2, ...an/m) for the individual steps involved in the process were 
homogeneous and no thresholds were present from one step to the next, 
then , no discontinuities in the increase of the Hill coefficient (i.e. in the 
change of k) with R1 should be observed. The existence of a major 
5 heterogeneity among the k\ values corresponding to the microscopic 
individual steps (i.e. the existence of thresholds for the intermediate 
steps) might lead to a macroscopic discontinuity in the system. 
Heterogeneity would cause a change in the rate of variation of the Hill 
coefficient and, which would require a jump in the macroscopic value of k 

10 in order for equation (2) to fit the data. 

The presence of internal heterogeneity in the reaction process can 
be detected by the appearance of steps in the rate of change of the Hill 
coefficient, corresponding to the Hill curve that fits the experimental data. 
rjjs defined as an index of heterogeneity and its value corresponds to the 

15 number of steps in the rate of variation of the Hill coefficient (one step, 
#7=1; two steps, q = 2; n steps, /7= n). 

#7 is a key parameter for the dissection and detailed analysis of the 
reaction process. It is useful for the independent optimization and 
development of every one of the steps identified by rj- 

20 As mentioned, the presence of steps in the rate of change of n 

translates in an abrupt discontinuity in k. Therefore, every step is 
determined by a different macroscopic constant of resistance k. Systems 
with #7 = 2, can thus be described by a Hill equation in which k takes two 
different values (*1 and k2), according to the R1 interval considered. One 

25 part of the curve is described by /c1 and the other by /c2. 

Hill curves describing reaction processes characterized by /7 = 2, are 
hybrids generated from two parallel Hill curves differing only in k. The 
transition from one curve to the other may alter the smooth change in the 
slope of the resulting Hill curve. 
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The apparent titer r 

In the Hill equation (2), when R1 increases, r increases from 1 to 
2,3,4... and P approaches its P max value. On the other direction, on the 
contrary, R1 can only decrease up to a minimal point (R1 m i n h at which r 
5 and P reach their minimal values. The Hill sigmoidal curve is not 

symmetric, only the right arm is asymptotic (towards P max ). On the left 
arm, the curve has an origin at R1 min ; the empirical curve does not fit the 
data for values below R1 m } n . 

From a biological point of view, the fact that P does not exist for 
10 R1 below R1 min/ means that there is no 'product' when the concentration 
of 'substrate' is lower than R1 mln ; e.g. that the system is not responsive 
to concentrations below R1 mln . The minimal concentration of R1 that the 
system can detect and report is R1 min . 

In terms of biological agents, R1 mln represents the minimal amount 
15 that can elicit a response in a given reporter system, and it is represented 
by r. The titer defined this way, is neither an asymptote value nor a value 
approached by extrapolation, but a precise parameter of the Hill equation, 
at the very mathematical origin of the curve. 

t measures the limiting dilution or apparent titer of the biological 
20 reagent. The value of r is determined by the limit of sensitivity of the 
biological assay system and of method used for the measurement of the 
product P; that is why it is said to be apparent titer. 

r is specific to the batch or stock of the biological reagent tested, r 
represents the apparent concentration of the biological agent and is 
25 expressed in units per volume, e.g. the maximal dilution of the biological 
agent that leads to the production of P. r is given by the maximal R1 for 
which the Hill coefficient reaches its minimal value (the Hill coefficient 
becomes constant at a value equal or close to 1). The concept of r 
corresponds to that of titer, of general use for viruses, antibodies and 
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vectors. Variations in r affect equation (2) by shifting the curve to the 
right or to the left, according to whether the value of r decreases or 

* 

increases, respectively. 

t is a key parameter that measures the 'apparent' concentration of 
5 a stock of the biological agent, which is necessary for whatever use it will 

be given. 

The absolute titer 0 

0 is a the parameter that measures the absolute concentration 
(titer) of a stock or batch of the biological agent. The value of 9 is not 
10 determined by nor dependent on the limit of sensitivity of the biological 
assay system or of the method used for the measurement of the product 
P; that is why it is said to be absolute titer. Q is specific to the batch or 
stock of the biological reagent tested. It represents the real physical 
concentration of the biological agent and is expressed in units per volume, 
15 e.g. the maximal dilution of the biological agent that leads to the 
production of P. 

6 is given by the following equation 
Snarls (3),. 
where s is the sensitivity of the detection method. Therefore, for agents 
20 detected using the same method, the following expression is valid: 

Q\ n \ / r i = 02/72 It2 = enrmlTn = constant (4) 
Using equation (4), the ratio of the absolute titer 0, corresponding 
to two biological agent preparations, can be obtained from their 
respective n and r. Variations in 6 affect the equation (2) by shifting the 
25 curve to the right or to the left and/or by changing its slope. 

Compensation between rr and k 

/7 and k may appear to compensate to generate two different Hill 
curves (one differing in rr and the other one differing in k) that would 
apparently fit with the same experimental data. As rr and k have opposite 
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effects, two Hill curves; in which the increase in rr is compensated by the 
decrease in /c, and vice versa, may seem to represent the same curve, 
which could make it difficult to determine whether two Hill curves are 
different because a change in rr or in /c. 
5 Detailed analysis of the Hill curves indicates that rr and k do not 

compensate very well. Although curves differing in compensatory values 
of either rr or k may vary close each other, they do not fit exactly in any 
of the two regions of highest curvature (before and after the inflection 
point). This dispersion is caused by the fact that rr, but not /c, changes 
10 the slope at the inflection point of the Hill curve. Therefore, e, which is 
the slope of the Hill curve at the inflection point, can be used to easily 
differentiate between two Hill curves that apparently compensate for rr 

♦ 

and k. 

Conclusions 

15 The application of the Hill analysis to resolve complex biological 

processes is effective for the precise and objective understanding of 
processes like virus or vaccine action, entry, genome replication, 
transgene expression, vector/transgene immunogenicity, cytotoxicity and 
other such parameters. The analysis is independent of the virus vaccine, 

20 vector and protein type involved and from the output parameter and 

variable measured, such as the internalized vector DNA, transgene mRNA 
level and transgene product activity. 

As in the field of chemical pharmaceuticals, the structure of the 
potential drug (in this case the biological agent) must be optimized to a 

25 maximal possible intrinsic potency. In analytical development, the goal is 
to search for better performing reporter systems (the lowest possible k), 
as analytical tool. Two different systems characterized by constants kA 
and kB, respectively, can be compared (using the same biological agent) 
for their relative resistance or performance. 
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Complex systems involving the interaction of biological agents, 
such as viruses, vaccines, gene transfer vectors, antibodies proteins and 
living cells (either in vitro or in vivo) can be analyzed using the Hill 
equation. A complex succession of unitary processes, each of them 
5 susceptible to be individually analyzed by the Hill equation, as a global 
process, can be also described by the same equation as its constitutive 
steps. 

EXAMPLE 4 

Materials and Methods 
10 Cells: 

293 human embryo kidney (HEK) cells, obtained from ATCC, were 
cultured in Dulbecco's modified Eagle's medium containing 4.5 g/l 
glucose (DMEM; GIBGO-BRL) 10 % fetal bovine serum (FBS, Hyclone). 
Hela rep-cap 32 cells, described above, were obtained from Anna Salvetti 
15 (CHU, Nantes) and cultured in the medium described above. 

Plasmids: 

pNB-Adeno, which encodes the entire E2A and E4 regions and VA 
RNA I and II genes of Adenovirus type 5, was constructed by ligating into 
the polylinker of multiple cloning site of pBSII KS ( + /-) (Stratagene, San 

20 Diego, USA) the Sall-Hindlll fragment (9842-1 1555 nt) of Adenovirus 
type 5 and the BamHI-Clal fragment (21563- 35950) of pBR325. All 
fragments of adenovirus gene were obtained from the plasmid pBHG-10 
(Microbix, Ontario, Canada). pNB-AAV encodes the genes rep and cap of 
AAV-2 and was constructed by ligation of Xbal-Xbal PCR fragment 

25 containing the genome of AAV-2 from nucleotide 200 to 4480 into Xbal 
site of polylinker MCS of pBSIIKS( + /-). The PCR fragment was obtained 
from pAV1 (ATCC, USA). Plasmid pNB-AAV was derived from plasmid 
pVA1l, which contains the AAV genomic region, rep and cap. pNB-AAV 
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does not contain the AAV ITR's present in pAV1. pAAV-CMV(nls)LacZ 
was provided by Dr Anna S^lvetti (CHU, Nantes). 

pCMV{nls)LacZ (rAAV vector plasmid) and pNB-Adeno were 
prepared in DH5a E. coli and purified by Nucleobond AX PC500 Kit 
5 (Macherey-Nagel), according to standard procedures. Plasmid pAAV- 
CMV(nls)LacZ is derived from plasmid psub201 by deleting the rep-cap 
region with SnaB I and replacing it with an expression cassette harboring 
the cytomegalovirus (CMV) immediate early promoter (407 bp), the 
nuclear localized yff-galactosidase gene and the bovine growth hormone 
10 polyA signal (324 bp) (see, Chadeuf et al. (2000) J. Gene Med. 2:260- 
268. pAAV-CMV(nls)LacZ was provided by Dr Anna Salvetti. 

Virus: 

Wild type adenovirus (AV) type 5 stock, originally provided by Dr 
Philippe Moullier (CHU, Nantes), was produced accordingly to standard 
15 procedures. 

Construction of Rep mutant libraries 

25 pmol of each mutagenic primer was placed into a 96 PCR well 
plate. 15 //I of reaction mix (0.25 pmol of pNB-AAV), 25 pmol of the 
selection primer (changing one non-essential unique restriction site to a 

20 new restriction site), 2 //I of 1 0 X mutagenesis buffer (100 mM Tris- 
acetate pH7.5, 100 mM MgOAc and 500 mM KOAc pH7.5) was added 
into each well. The samples were incubated at 98°C for 5 minutes and 
then immediately incubated for 5 minutes on ice. Finally, the plate was 
placed at room temperature for 30 minutes. 

25 The primer extension and ligation reactions of the new strands 

were completed by adding to each sample: 7 jjI of nucleotide mix (2.86 
mM each nucleotide and 1 .43 X mutagenesis buffer) and 3jj\ of a fresh 
1:10 enzyme dilution mix (0.025 U///I of native T7 DNA polymerase and 
1 U///I of T4 DNA ligase were diluted in 20 mM Tris HCI pH7.5, 10 mM 
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KCI, 10 mM 0- mercaptoethanol, 1 mM DTT, 0.1 mM EDTA and 50% 
glycerol). Samples were incubated at 37 °C for 1 hour. The T4 DNA ligase 
was inactivated by incubating the reactions at 72°C for 15 minutes to 
prevent re-ligation of the digested strands during the digestion of the 
5 parental plasmid (pNB-AAV). 

Each mutagenesis reaction was digested with restriction enzyme to 
eliminate parental plasmids: 30 //I solution containing 3 jjI of 10X enzyme 
digestion buffer and 10 units of restriction enzyme were added to each 

i 

mutagenesis reaction and incubated at 37 °C for at least 3 hours. 
10 90 //I of the £. coli XLmutS competent cells (Stratagene, San Diego 

CA; supplemented with 1.5 //I of ^mercaptoethanol to a final concen- 
tration of 25 mM) were aliquoted into prechilled deep-well plates. The 
plates were incubated on ice for 10 minutes and swirling gently every 2 
minutes. 

15 A fraction of the reactions that had been digested with restriction 

enzyme (1/10 of the total volume) was added to the deep well plates. The 
plates were swirled gently prior to incubation on ice for 30 minutes. A 
heat pulse was performed in a 42°C water bath for 45 seconds, the 
transformation mixture was incubated on ice for 2 minutes and 0.45 ml of 

20 preheated SOC medium (2% (w/v) tryptone, 0.5% (w/v) yeast extract, 
8.5 mM NaCI, 2.5 mM KCI, 10 mM MgCI 2 and 20 mM glucose at pH 7) 
was added. The plates were incubated at 37°C for 1 hour with shaking. 

To enrich for mutant plasmids, 1 ml of 2X YT broth medium (YT 
medium is 0.5% yeast extract, 0.5% NaCI, 0.8% bacto-tryptone), 

25 supplemented with 1 00 /yg/ml of ampicillin, was added to each 

transformation mixture and the cultures were grown overnight at 37 °C 
with shaking. Plasmid DNA isolation was performed from each mutant 
culture using standard procedure described in Nucleospin Multi-96 Plus 
Plasmid Kit (Macherey-Nagel). Five hundred /jg of the resulting isolated 
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DNA was digested with 10 units of the selection restriction enzyme in a 
total volume of 30/yl containing 3 fj\ of 10X enzyme digestion buffer for 

overnight at 37°C. 

A fraction of the digested reactions (1/10 of the total volume) were 
5 transformed into 40 jjl of Epicurian coli XL1 -Blue competent cells 

supplemented with 0.68 /j\ of /ff-mercaptoethanol to a final concentration 
of 25 mM. After heat pulse, 0.45 ml of SOC was added and the 
transformation mixtures were incubated for 1 hour at 37°C with shaking 
before to be plate on LB-ampicillin agar plates. The agar plates were 
10 incubated overnight at 37 °C and the colonies obtained were picked up 
and grown overnight at 37 °C into deep-well plates. 

Four clones per reaction were screened for the presence of the 
mutation using restriction enzyme specific to the new restriction site 
introduced into the mutated plasmid with the selection primer. The cDNA 
15 from selected clones was also sequenced to confirm the presence of the 
expected mutation. 

Monitoring rAAV Production 

rAAV from each of the above wells, were produced by triple 
transfection on 293 HEK cells. 3 x 10 4 cells were seeded into each well 

20 of 96 micro-well plate and cultured for 24 hours before transfection. 
Transfection was made on cells at about 70% confluency. 25 kDa PEI 
(poly-ethylene-imine, Sigma-Aldrich) was used for the triple transfection 
step. Equimolar amounts of the three plasmids AV helper plasmid (pNB- 
Adeno), AAV helper plasmid (pNB-AAV or a mutant clone rep plasmid) 

25 and vector plasmid (pAAV-CMV(nls)LacZ) were mixed with 10 mM PEI by 
gently shaking. The mixture was the added to the medium culture on the 
cells. 60 hours after transfection, the culture medium was replaced with 
100 //I of lysis buffer (50 mM Hepes, pH 7.4; 150 mM NaCI; 1 m MgCI 2 ; 
1 mM CaCI 2 ; 0.01 % CHAPS). After one cycle of freeze-thawing the 
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cellular lysate was filtered through a millipore filter 96 well plate and 
stored at -80° C. 

rAAV infection particles (ip) 

Titers of rAAV vector particles were determined on HeLa rep/cap 
5 32 cells using standard dRA {serial dilution replication assay) test. Cells 
were plated 24 hours before infection at a density of 1 x 1 0 4 cells in 96- 
well plates. Serial dilutions of the rAAV preparation were made between 1 
and 1 x 10 6 jj\ and used for co-infection of the HeLa rep/cap 32 cells 
together with wt-AV type 5 (MOI 25). 48 hours after infection the ip 
10 were measured by real time PCR or by the quantification of biological 
activity of the transgene. 
Real Time PCR 

Infected HeLa rep/cap 32 cells were lysed with 50 jj\ of solution 
(50 mM Hepes, pH 7.4; 150 mM NaCI). After one cycle of freeze-thawing 

15 50 jj\ of Proteinase K (10 mg/ml) and the lysate were incubated one hour 
at 55 °C. The enzyme was inactivated by incubation 10 min at 96°C. 

For real time PCR, 0.2 //I of lysate was taken. Final volume of the 
reaction was 10 //I in 384 well plate using an Applied Biosystem Prism 
7900. The primers and fluorescence probe set corresponding to the CMV 

20 promoter were as follows: CMV 1 primer 5'- 

TGCCAAGTACGCCCCCTAT-3' (SEQ ID No. 733) (0.2 /yM) and CMV 2 
primer 5 ' - AG GTC ATGTACTG G G C AT A ATG C -3' (SEQ ID No. 734) (0.2 
/yM) ; probe VIC-Tamra 5'-TCAATGACGGTAAATGGCCCGCCT-3' (SEQ 
ID No. 735) (0.1 jjM). dRA plots were obtained by plotting the DNA copy 

25 number (obtained by real time PCR) vs. the dilution of the rAAV 
preparation. 
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/?-Galactosidase activity 

After 48 hours of infection, cells were treated with trypsin, and 
100 //I of reaction solution (GalScreen Kit, Tropix) was added and 
incubated for one hour at 26 °C. Luminescence was measured in 
5 NorthStar (Tropix) HTS station. dRA plots were obtained plotting the 
intensity of £-Galatosidase activity vs. the dilution of the rAAV 
preparation. 

Mathematical Model for results analysis: 

Results were analyzed using the Hill equation-based analysis 

10 (designated NautScan™; see, Patent n° 9915884, 1999, France; 

published as International PCT application No. WO 01/44809 (PCT n° 
PCT/FR00/03503, Dec, 2000, see EXAMPLE above). Briefly, data were 
processed using a Hill equation-based model that allows extraction of key 
feature indicators of performance for each individual mutant. Mutants 

15 were ranked based on the values of their individual performance and 
those at the top of the ranking list were selected as Leads. 
Results 

Generation of diversity. 

To identify candidate amino acid (aa) positions on the rep protein 
20 involved in rep protein activity an Ala-scan was performed on the rep 
sequence. For this, each amino acid in the rep protein sequence was 
replaced with Alanine. To do this sets of rAAV that encode mutant rep 
proteins in which each differs from wild type by replacement of one 
amino acid with Ala, were generated. Each set of rAAV was individually 
25 introduced into cells in a well of a microtiter plate, under conditions for 
expression of the rep protein. The amount of virus that could be 
produced from each variant was measured as described below. Briefly, 
activity of Rep was assessed by determining the amount of AAV or rAAV 
produced using infection assays on HeLa Rep-cap 32 cells and by 
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10 



15 



20 



measurement of AAV DNA replication using Real Time PCR, or by 

9 

assessing transgene (£-galactosidase) expression. The relative activity of 
each individual mutant compared to the native protein was assessed and 
"hits" identified. Hit positions are the positions in the mutant proteins 
that resulted in an alteration (selected to be at least about 20%), in this 
instance all resulted in a decrease, in the amount of virus produced 
compared to the activity of the native (wildtype) gene (see Fig, 3A). 

The hits were then used for identification of leads (see, Fig. 3B). 
Assays for Rep activity were performed as described for identification of 
the hit positions. Hit positions on Rep proteins and the effect of specific 
amino acids on the productivity of AAV-2 summarized in the following 
table: 



Hit position 


replacing amino acid (effect) 


4 (ttt) F 


(get) A (decrease) 




10 (aag) K 


(gcg) A (decrease) 




20 (ccc) P 


(gec) A (decrease) 




22 (att) I 


(get) A (decrease) 




28 (tgg) W 


(gcg) A (decrease) 




32 (gag) E 


(gcg) A (decrease) 




38 (ccg) P 


(gcg) A (decrease) 




39 (cca) P 


(gca) A (decrease) 




54 (ctg) L 


(get) A (decrease) 




59 (ctg) L 


(gcg) A (decrease) 




64 (ctg) L 


(gcg) A (decrease) 




74 (ccg) P 


(gcg) A (decrease) 




86 (gag) E 


(gcg) A (decrease) 




88 (tac) Y 


(gec) A (decease) 




101 (aaa) K 


(gca) A (decrease) 
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Hit position 


replacing amino acid (effect) 


1 24 (ate) I 


(gee) A (decrease) 




125 (gag) E 


(gcg) A (decrease) 




1 27 (act) T 


(get) A (decrease) 




132 (ttc) F 


(gee) A (decrease) 




1 40 (ggc) G 


(gee) A (decrease) 




161 (acc) T 


(gee) A (decrease) 




1 63 (cct) P 


(get) A (decrease) 




1 75 (tat) Y 


(get) A (decrease) 




193 (ctg) L 


(gcg) A (decrease) 




196 (gtg) V 


(gcg) A (decrease) 




1 97 (teg) S 


(gee) A (decrease) 




221 (tea) S 


(gca) A (decrease) 




228 (gtc) V 


(gcg) A (decrease) 




231 (etc) L 


(gee) A (decrease) 


• 


234 (aag) K 


(gcg) A (decrease) 




237 (acc) T 


(gee) A (decrease) 




250 (tac) Y 


(gee) A (decrease) 




258 (aac) N 


(gee) A (decrease) 




260 (egg) R 


(gcg) A (decrease) 




263 (ate) I 


(gee) A (decrease) 




264 (aag) K 


(gcg) A (decrease) 




334 (ggg) G 


(gcg) A (decrease) 




335 (cct) V 


(get) A (decrease) 




337 (act) T 


(get) A (decrease) 




341 (acc) T 


(gec) A (decrease) 




342 (aac) N 


(gee) A (decrease) 
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Hit position 


replacing amino acid (effect) 


347 (ata) 1 


(gca) A (decrease) 




350 (act) T 


(get) A (decrease 


(aat) N (increase) 


354 (tac) Y 


(gec) A (decrease) 




363 (aac) N 


(gec) A (decrease) 




364 (ttt) F 


(get) A (decrease) 




367 (aac) N 


(gee) A (decrease) 




370 (gtc) V 


(gec) A (decrease) 




376 (tgg) W 


(gcg) A (decrease) 




381 (aag) K 


(gcg) A (decrease) 




382 {atg) M 


(gcg) A (decrease) 




389 (teg) S 


(gcg) A (decrease) 




407 (tec) S 


(gec) A (decrease) 




41 1 (ata) I 


(gca) A (decrease) 




414 (act) T 


(get) A (decrease) 




420 (tec) S 


(get) A (decrease) 




421 (aac) N 


(gec) A (decrease) 




422 (acc) T 


(gee) A (decrease) 




424 (atg) M 


(gcg) A (decrease) 




428 (att) I 


(get) A (decrease) 




429 (gac) D 


(gec) A (decrease) 




438 (cag) Q 


(gcg) A (decrease) 




440 (ccg) P 


(gcg) A (decrease) 




451 (acc) T 


(gec) A (decrease) 




460 (aag) K 


(gcg) A (decrease) 




462 (acc) T 


(gee) A (decrease) 


(ata) I (increase) 


484 (ttc) F 


(gee) A (decrease) 
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Hit position 


replacing amino acid (effect) 


488 (aag) K 


(gcg) A (decrease) 




495 (ccc) P 


(gec) A (decrease) 




497 (ccc) P 


* 

(gec) A (decrease) 


(cga) R (increase) 


497 (ccc) P 


(gec) A (decrease) 


(etc) L (increase) 


497 (ccc) P 


(gec) A (decrease) 


(tac) Y (increase) 


498 (agt) S 


(get) A (decrease) 




499 (gac) D 


(gec) A (decrease) 


- 


503 (agt) S 


(gcg) A (decrease) 




511 (tea) S 


(gca) A (decrease) 




512 (gtt) V 


(get) A (decrease) 




516 (teg) S 


(gcg) A (decrease) 




517 (acg) T 


(get). A (decrease) 


(aac) N (increase) 


518 (tea) S 


(gca) A (decrease) 




519 (aac) D 


(gcg) A (decrease) 




542 (ctg) L 


(gcg) A (decrease) 


(teg) S (increase) 


548 (aga) R 


(gca) A (decrease) 


(age) S (increase) 


598 (gga) G 


(gca) A (decrease) 


(gac) D (increase) 


598 (gga) G 


(gca) A (decrease) 


(age) S (increase) 


600 (gtg) V 


(gcg) A (decrease) 


(ccg) P (increase) 


601 (cca) P 


(gca) A (decrease) 




Hit position 
(within intron) 


replacing sequence (effect) 


630 (tgc) 


gcg (decrease) 


cgc or tea or cct 
(increase) 



25 The hits in other AAV serotypes (see, also FIGS. 5A and 5B) are as 

follows: 



HIT POSITION 
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5 



10 



15 



25 



30 



AAV- 2 


AAV-1 


AAV-3 


AAV-3 B 


AAV-4 


AAV-6 


AAV-5 


4 


4 


4 


4 


4 


4 


4 


10 


10 


10 


10 


10 


10 


10 


20 


20 


20 


20 


20 


20 


20 


22 


22 


22 


22 


22 


22 


22 


29 


29 


29 


29 


29 


29 


29 


■t 

32 


32 


32 


32 


32 


32 


32 


38 


38 


38 


38 


38 


38 


38 


39 


39 


39 


39 


39 


39 


39 


54 


54 


54 


54 


54 


54 


54 


59 


59 


59 


59 


59 


59 


59 


64 


64 


64 


64 


64 


64 


64 


74 


74 


74 


74 


74 


74 




86 


86 


86 


86 


86 


86 


85 


88 


88 


88 


88 


88 


88 


87 


101 


101 


101 


101 


101 


101 


100 


124 


124 


124 


124 


124 


124 


123 


125 


125 


125 


125 


125 


125 


124 


127 


127 


127 


127 


127 


127 


126 


132 


132 


132 


132 


132 


132 


131 


140 


140 


140 


140 


140 


140 




161 


161 


161 


161 


161 


161 


1 58 


163 


163 


163 


163 


163 


163 


160 


175 


175 


175 


175 


175 


175 


172 


193 


193 


193 


193 


193 


193 


190 


196 


196 


196 


196 


196 


196 


193 


197 


197 


197 


197 


197 


197 


194 


221 


221 


221 


221 


221 


221 


217 


228 


228 


228 


228 


228 


228 


224 


231 


231 


231 


231 


231 


231 


227 
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10 



15 



234 


234 


234 


234 


234 


234 


230 


237 


237 


237 


237 


237 


237 


233 


250 


250 


250 


250 


250 


250 


246 


258 


258 


258 


258 


258 


258 


254 


260 


260 


260 


260 


260 


260 


256 


263 


263 


263 


263 


263 


263 


259 


264 


264 


264 


264 


264 


264 


260 


334 


334 


334 


334 


334 


334 


330 


335 


335 


335 


335 


335 


335 


331 


337 


337 


337 


337 


337 


337 


333 


341 


341 


341 


341 


341 


341 


337 


342 


342 


342 


342 


342 


342 


338 


347 


347 


347 


347 


347 
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342 


350 


350 


350 


350 


350 


350 


346 


354 


354 
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354 


354 


350 


363 
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363 


363 


363 


363 


359 


364 


364 


364 


364 


364 


364 


360 


367 


367 


367 


367 


367 


367 


363 


370 


370 


370 


370 


370 


370 


366 


376 


376 


376 


376 


376 


376 


372 


381 


381 


381 


381 


381 


381 


377 


382 


382 


382 


382 


382 


382 


378 


389 


389 


389 


389 


389 


389 


385 


407 


407 


407 


407 


407 


407 


403 


411 


411 


411 


411 


411 


411 


407 


414 


414 


414 


414 


414 


414 


410 


420 


420 


420 


420 


420 


420 


416 


421 


421 


421 


421 


421 


421 


417 


422 


422 


422 


422 


422 


422 


418 


424 


424 


424 


424 


424 


424 


420 . 



20 



25 



30 
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428 


428 


428 


428 


428 


428 


424 


429 


429 


429 


429 


429 


429 


425 


438 


438 


438 


438 


438 


438 


434 


440 


440 


440 


440 


440 


440 


436 


451 


451 


451 


451 


451 


451 


447 


460 


460 


460 


460 


460 


460 


456 


462 


462 


462 


462 


462 


462 


458 


484 


484 


484 


484 


484 


484 


480 


488 


488 


488 


488 


488 


488 


484 


495 


495 


495 


495 


495 


495 


491 


497 


497 


497 


497 


497 


497 


493 


498 


498 


498 


498 


498 


498 


494 


499 


499 


499 


499 


499 


499 


495 


503 


. 503 


503 


503 


503 


503 


499 


511 


511 


511 


511 


511 


511 


529 


512 


512 


512 


512 


512 


512 


530 


516 


516 


516 


516 


516 


51 6 


534 


517 


517 


517 


517 


517 


517 


535 


518 


518 


518 


518 


518 


518 


536 


519 


519 


519 


519 


519 


519 


537 


542 


543 


542 


542 


542 


543 


561 


548 


549 


548 


548 


548 


549 


567 


598 


599 


600 


600 


599 


599 




600 


602 


603 


603 


602 


602 


589 


601 


603 


604 


604 


603 


603 


590 



30 



Sets of nucleic acids encoding the rep protein were generated. The 
rep proteins encoded by these, sets of nucleic acid molecules were those 
in which each amino acid position identified as a "hit" in the ala-scan 
step, were each sequentially replaced by all remaining 18 amino acids 
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using site directed mutagenesis. Each mutant was designed, generated, 
processed and analyzed physically separated from the others in 
addressable arrays. No mixtures, pools, nor combinatorial processing 
were used, 

5 As in the first round (alanine scan), a library of mutant rAAV was 

generated in which each individual mutant was independently and 
individually generated in a independent reaction and such that each 
mutant contains only a single amino acid change and this for each amino 
acid residue. Again, each resulting mutant rep protein was then 

10 expressed and the amount of virus produced in cells assessed and 
compared to the native protein. 
Lead identification 

Since rep proteins that result in increased virus production are of 
interest, those mutants that lead to an increase in the amount of virus 

15 produced (2 to 10 times the native activity), were selected as "leads." 
Ten such mutants were identified. 

Based on the results obtained from the assays described above (i.e. 
titer of virus produced by each rep variant), each individual rep variant 
was assigned a specific activity. Those variant proteins displaying the 

20 highest titers were selected as leads (see Table above). Leads include: 
amino acid replacement of T by N at Hit position 350; T by I at Hit 
position 462; P by R at Hit position 497; P by L at Hit position 497; P by 
Y at Hit position 497; T by N at Hit position 517; G by S at Hit position 
598; G.by D at Hit position 598; V by P at Hit position 600. 

25 Also provided are combinations of the above mutant Rep 78, 68, 

52 and 40 proteins, nucleic acids encoding the proteins, and recombinant 
AAV (any serotype) containing the mutation at the indicated position or 
corresponding position for serotypes other than AAV-2, including any set 
forth in the following table and corresponding SEQ ID Nos. Each amino 
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acid sequence is set forth in a separate sequence ID listing; for each 
mutation or combination thereof there is a single SEQ ID setting forth the 
unspliced nucleic acid sequence for Rep78/68, which for all mutations 
from amino acid 228 on, includes the corresponding Rep 52 and Rep 40 
5 encoding sequence as well. 

Amino acid sequences of exemplary mutant Rep proteins 





Qan nn 


aene 


Dosition(s) 


codon(s) 




seq.1 


rep78 


4 


GCT 




seq.2 


rep68 


4 


GCT 


10 


seq.3 


rep78 


10 


GCG 




seq.4 


rep68 


10 


GCG 




seq.5 


rep78 


20 


GCC 




seq.6 


rep68 


20 


GCC 




seq.7 


rep78 


22 


GCT 


15 


seq.8 


rep68 


22 


GCT 




seq.9 


rep78 


29 


GCG 




seq.10 


rep 68 


29 


GCG 




seq.1 1 


rep78 


38 


GCG 




seq.1 2 


rep68 


38 


GCG 


20 


seq.1 3 


rep78 


39 


GCA 




seq.1 4 


rep68 


39 


GCA 




seq.1 5 


rep78 


53 


GCT 




seq.1 6 


rep68 


53 


GCT 




seq.1 7 


rep78 


59 


GCG 


25 


seq.1 8 


rep68 


59 


GCG 




seq.1 9 


rep78 


64 


GCT 




seq.20 


rep68 


64 ■ 


GCT 




seq.21 


rep78 


74 


GCG 




seq.22 


rep68 


74 


GCG 


30 


seq.23 


rep 7 8 


86 


GCG 




seq.24 


rep68 


86 


GCG 




seq.25 


rep 7 8 


88 


GCC 




seq.26- 


rep68 


88 


GCC 




seq.27 


rep78 


101 


GCA 


35 


seq.28 


rep68 


101 


GCA 




seq.29 


rep 7 8 


124 


GCC 




seq.30 


rep68 


124 


GCC 




seq.31 


rep78 


125 


GCG 




seq.32 


rep68 


125 


GCG 


40 


seq.33 


rep78 


127 


GCT 




seq.34 


rep68 


127 


GCT 




seq.35 


rep 7 8 


132 


GCC 




seq.36 


rep68 


132 


GCC 




seq.37 


rep78 


140 


GCC 


45 


seq.38 


rep68 


140 


GCC 
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con 3Q 
acq .03 


ran7R 


161 


GCC 

VJ v* v* 




pan Afl 


lopuo 


161 

1 VJ 1 


GCC 




0 a n "A 1 

seq.'f 1 


rep / O 


1 63 


GCT 

VJ V* 1 




0 a/-i ^1_9 

seq.f z 


repoo 


1 63 


GCT 




seq .*f 0 


rep / 0 


1 75 


GCT 

VJ V— ' 1 




seq.*t-*r 


repoo 


I/O 


GCT 




seq .h-o 


rep / 0 


1 93 

I Ow 


GCG 




0 An /LA 

seq .*+o 




1 93 
■ w w 


GCG 

\J 




seq .h- / 


rep / 0 


1 96 


GCC 


10 


seq.*t 0 


repoo 


1 96 

I w VJ 


GCC 




can /LQ 

seq .H-y 


ran7ft 

rep / 0 


1 97 


GCC 

VJ V* V* 




con K C\ 

seq . ovj 


repoo 


1 97 


GCC 




0 an R 1 

seq.o 1 


ron7fi 

rep / 0 


991 


GCA 

vj v*r^ 




seq.oz 


ran ft ft 

repoo 


991 


GCA 

VJ v*/*^ 




seq. 00 


rep / 0 


99ft 


GCG 

VJ VVJ 




seq. oh- 


ran K O 

repoz 


99fl 

ZZO 


GCG 

VJ W VJ 




seq. do 


repoo 


99S 


GCG 

VJ V» VJ 




seq . 00 


ran 

rep*rvj 


99ft 


GCG 

VJ vw 




*"* <""V /~l M 

seq.o / 


rep / 0 


931 
zo 1 


GCC 

VJ w w 




ft Aft C O 

seq. 00 


ran K O 

repoz 


931 

ZO I 


GCC 

VJ WW 




can CT O 

seq. oa 


repoo 


931 

ZO I 


GCC 

VJ W Vi» 




seq.ou 


repfu 


931 


GCC 

VJ WW 




seq.o i 


rep / 0 


93/1 


GCG 

VJ WVJ 




seq .oz 


ran K 9 

repoz 


9 3 J. 


GCG 

vj wVJ 


Z.O 


seq .00 


repoo 


9 3A 

ZOH* 


GCG 

VJ W VJ 




seq. oh- 


rep*+u 


934. 
ZO*r 


GCG 

vj w vj 




seq. 00 


rep / 0 


937 

^0 / 


GCC 

VJ w w 




seq. 00 


repoz 


937 


GCC 

VJ ww 




inn C ~7 

seq.o / 


ran CP 

repoo 


937 
zo / 


GCC 

vj w w 


on 


seq. 00 


rep*nj 


937 
ZO / 


GCC 

VJWW 




seq.oy 


rep to 


zow 


GPP 

VJ WW 




seq. /u 


repoz 


ZOVJ 


GCC 

www 




seq. / 1 


MnfiQ 

repoo 


ZOVJ 


GPP 

www 




seq. /z 


rep*fu 


ZOVJ 


Vj ww 


oc 


seq. /o 


rep / a 


ZOO 


GPP 




seq. 74 


repoz 


ZOO 


fiPP 




seq. 75 


repoo 


ZOO 


GPP 

VJ WW 




seq. 76 


rep40 


258 


GCC 




seq. 77 


rep78 


260 


GCG 


40 


seq. 78 


rep52 


260 


GCG 




seq. 79 


rep68 


260 


GCG 




seq. 80 


rep40 


260 


GCG 




seq. 81 


rep78 


263 


GCC 




seq. 82 


rep52 


263 


GCC 


45 


seq. 83 


rep68 


263 


GCC 




seq. 84 


rep40 


263 


GCC 




seq. 85 


rep78 


264 


GCG 




seq. 86 


rep52 


264 


GCG 




seq. 87 


rep68 


264 


GCG 


50 


seq. 88 


rep40 


264 


GCG 




seq. 89 


rep78 


334 


GCG 
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con Q/"\ 

seq .su 




334. 


GCG 




an O 1 

seq.w i 


I CfJUO 


334. 

W W*T 


GCG 

Vl^ 




seq.yz 




334. 

00*T 


GCG 




can 

acl(. 30 


rpn78 


335 

www 


GCT 


5 


can Q /I 




335 

WW w 


GCT 




can 


rpn68 

1 G fJ W W 


335 

WWW 


GCT 




can Qfi 


rpr>40 

1 wjJ^rw 


335 

www 


GCT 




con Q 7 

OClj.v? / 


rpn78 

1 C / w 


337 

WW / 


GCT 




can Q R 




337 

WW / 


GCT 


10 


can QQ 

seq 




337 

WW / 


GCT 




con 1 f^O 

seq. i uu 




337 

WW / 


GCT 




can 1 0 1 

seq . i u i 


rep / o 


34.1 

0*T 1 


GCC 




seq . i uz 




34.1 


GCC 




seq . i uo 


repoo 


34.1 


GCC 


1 o 


seq . i u*f 


rep*ru 


34.1 

w*+ 1 


GCC 




seq. i uo 


Top / O 


349 


GCC 




can 1 f\f\ 

seq. 1 uo 




34.9 


GCC 




can i II 

seq. 1 u / 


rpnflft 

(CjJUO 


34.9 


GCC 




can 1 Or Q 

seq. i uo 


f cptv 


34.9 


GCC 




seq* i uy 


rep / o 


34.7 


GCA 




seq. 1 i u 


ran C O 


34.7 

0*r / 


GCA 

\J v*^r^ 




seq. I i i 


repoo 


34.7 


GCA 




seq.i l z. 


r ep***u 


34.7 
w*r / 


GCA 




seq. i i o 


rep / o 


www 


AAT 




seq. i i *r 


ranC O 

repoz. 


o.cn 


AAT 




seq. i i o 


repoo 


3^n 


AAT 

nn i 




can 1 1 ft 

seq .mo 


rep*ru 


3RH 

O ww 


AAT 
^^^^ i 




seq. i i / 


ran "7 R 

rep / O 


3RD 
o uu 


GCT 




seq. i i o 


repoz, 


3*?n 

www 


GCT 




seq. r i y 


repoo 


3^n 

O wU 


GCT 




r««n 1 Ort 

seq. l zu 


rep^u 


oou 


GCT 




seq. 1 4. i 


rep /o 


oo*f 


CnCC 

VJ WW 




seq. i 


repoz 


*3R4. 
oo*f 


V3V-(V. 




seq. 1 Zo 


repoo 


OC/1 
wO*f 


err 




seq. i Z4 


rep^u 


OCA 


V3 w w 




seq.i zo 


rep / o 


ooo 


V9 O w 




r»on 1 OA 

seq. i zo 


ran K O 

rep oz 


ooo 


V3 w w 




seq. 127 


rep68 


363 


GCC 




seq. 128 


rep40 


363 


GCC 


40 


seq. 129 


rep78 


364 


GCT 




seq. 130 


rep 5 2 


364 


GCT 




seq. 131 


rep68 


364 


GCT 




seq. 132 


rep40 


364 


GCT 




seq. 133 


rep78 


367 


GCC 


45 


seq.i 34 


rep52 


367 


GCC 




seq. 1 35 


rep68 


367 


GCC 




seq. 1 36 


rep40 


367 


GCC 




seq. 137 


rep78 


370 


GCC 




seq. 138 


rep52 


370 


GCC 


50 


seq. 139 


rep68 


370 


GCC 




seq. 140 


rep40 


370 


GCC 
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10 



15 



20 



30 



35 



40 



45 



50 



Con 
oclj. 


1 A1 


ron7R 
rep / o 


376 


GCG 

Www 


con 
ocCj, 


1 A9 


ranC O 

rep o^. 


376 


GCG 


ocC{. 


1 A3 


rep do 


376 


u Ww 




1 /I /I 


rcp*rv 


376 


GCG 


DClj . 


1 AR 


rep / o 


381 


GCG 

www 


con 
ocL] . 


1 Afi 


ran Cj 9 


381 


GCG 
www 


con 

ocq. 


1 A 7 


ran Q 

repoo 


381 


GCG 

W WW 


con 

seq . 


1 Afl 


i ep*+v^ 


381 


GCG 

www 


con 

seq. 


1 AQ 


ran7R 
rep / o 


389 


GCG 

W WW 


con 


1 f»n 


rep o«c 


389 


GCG 

W WW 


con 

btJLj . 


19 1 


repoo 


389 


GCG 

W WW 


con 
acq . 


1 R9 


J CjJ'+v 


382 


GCG 

www 


con 

seq. 


1 R3 


ron7R 

rep / o. 


38Q 


GCG 

www 


con 

seq. 


I Of 


repo^ 


OOj 


GCG 

WWW 


seq. 


1 RS 
1 OO 


repoo 


3RQ 


GCG 
www 


con 

seq . 


1 OO 


ranAn 

rep*rw 


009 


GCG 
www 


con 

seq . 


1 R7 
1 o / 


rep / o 


AO 7 


GCC 
www 


con 

seq . 


1 oo 


repo^. 


Afi7 


GCC 
www 


con 

seq. 


1 0<7 


repoo 


A07 


GCC 
www 


con 

seq. 


1 OU 


rep*rw 


Afl7 


GCC 
www 


seq. 


1 O I 


rep / o 


A1 1 


GCA 


seq. 


169 
i OZ 


repoz. 


A1 1 


GCA 


jpt *m 

Seq. 


1 DO 


repoo 


4.1 1 


GCA 

W WSA 


con 

seq. 


I 6A 
1 0*f 


rep*Hj 


A1 1 


GCA 


seq. 


I DO 


rep / o 


A1 A 


GCT 

WW I 


seq. 


1 OO 


repoz 


A1 A 


WW I 


seq. 


1 O / 


repoo 




WV-r I 


seq. 


1 OO 


rep*tu 


A1 A 


WW I 


seq. 


1 AO 


rep /o 


Aon 


GCT 

WW I 


coo * 

seq. 


I /u 


repoz 


Aon 


WW I 


/> Art # 

Seq. 


I / I 


repoo 


A9H 


WW I 


seq. 


I 79 


■Tin yt ^1 

rep*t-u 


Aon 


GCT 

WW I 


seq. 


I 70 
I /o 


*»an~7 Q 

rep /o 


AO 1 


GCC 

WWW 


seq. 




repoz 


A0 1 


www 


seq. i 


f 7C 
1 /O 


repoo 




Www 


seq. 


1 "7 £5 
1 /O 


rep^u 




www 


seq." 


77 


rep 78 


yf OO 

4^^i 


www 


seq.1 


178 


rep52 


422 


GCC 


seq. 1 


179 


rep68 


422 


GCC 


seq.1 


180 


rep40 


422 


GCC 


seq.1 


81 


rep78 


424 


GCG 


seq.1 


82 


rep52 


424 


GCG 


seq.1 


83 


rep68 


424 


GCG 


seq.1 


84 


rep40 


424 


GCG 


seq.1 


85 


rep78 


428 


GCT 


seq.1 


86 


rep52 


428 


GCT 


seq.1 


87 


rep68 


428 


GCT 


seq.1 


88 


rep40 


428 


GCT 


seq.1 


89 


-rep78 


429 


GCC 


seq.1 


90 


rep52 


429 


GCC 


seq.1 


91 


rep68 


429 


GCC 
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sea 192 


reo40 


429 


GCC 




sea 193 


reo78 

1 W Km* m mmw 


438 


GCG 




sea 1 94 


reD52 

■ mm wmf mm* mm* 


438 


GCG 




sea 195 

www * I w W 


reD68 

I mm* m-w mm¥ mmw 


438 


GCG 


5 


sea 1 96 


reo40 


438 


GCG 




sea 1 97 


ren78 


440 


GCG 




sea 1 98 


T6d52 

V \-w mm* mm 


440 


GCG 




sea 1 99 


reo68 

■ \mf mdW mmw mmf 


440 


GCG 




can 900 


T6D40 


440 


GCG 


10 


<spa 201 


reD78 


451 


GCC 




can 202 


ran 5 2 

1 W f-f %mf --mm 


451 


GCC 




epn 903 


reo68 


451 


GCC 




can 904 


reD40 


451 


GCC 




sea 205 


reD78 


460 

• mr mr 


GCG 


15 


sea 206 

uykji mm W w 


reD52 


460 


GCG 




sea 207 


reo68 


460 


GCG 




sea. 208 


reo40 


460 


GCG 




sea 209 


reo78 


462 


GCC 




epn. 91 O 

tmm \ w 


reo52 

1 Wl^ \mf mm 


462 


GCC 


20 


sea 21 1 


ren68 


462 


GCC 




can 919 

oCLj.il 1 


reo40 


462 

■ \mf mm 


GCC 




can 913 


reo78 


462 


ATA 




can 914. 


reo52 


462 


ATA 




cpn 915 


ren68 


462 

■ mm 


ATA 


25 


can 916 


reo40 

1 w Lrf ^» 


462 


ATA 




sea 21 7 


reo78 


484 


GCC 




sea 21 8 


reo52 

| \mf +mf +m0 mm 


484 


GCC 




sea 21 9 


reo68 

4 W fmT \mW Kmf 


484 


GCC 




sea 220 


reD40 


484 


GCC 


30 


can 991 


reD78 


488 


GCG 




can 999 


1 ^^J*J mm 


488 


GCG 




can 993 


rpnfiR 


488 


GCG 




con 994. 


rpndO 

I C^*tv 


488 

"T w 


GCG 




can 99R 


rpn78 


495 


GCC 




con 99fi 


ronR9 


495 

^1 w w 


GCC 




con 997 


rpnfiR 


495 


GCC 




con 99R 


rpn40 


495 

W W 


GCC 




seq.229 


rep78 


497 


GCC 




seq.230 


rep52 


497 


GCC 


40 


seq.231 


rep68 


497 


GCC 




seq.232 


rep40 


497 


GCC 




seq.233 


rep78 


497 


CGA 




seq.234 


rep 5 2 


497 


CGA 




seq.235 


rep68 


497 


CGA 


45 


seq.236 


rep40 


497 


CGA 




seq.237 


rep78 


497 


CTC 




seq.238 


rep 5 2 


497 


CTC 




seq.239 


rep68 


497 


CTC 




seq.240 


rep40 


497 


CTC 


50 


seq.241 


rep78 


497 


TAC 




seq.242 


rep52 


497 


TAC 
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can 943 


reo68 


497 




con 94/1 
aoLj ■ *r 


reo40 


497 




can 945 


reo78 


498 




eon 946 


reo52 


498 


5 

w 


con 94.7 


T6D68 


498 




can 94.8 


reo40 


498 




cart 94. Q 


reo78 


499 




can 95fi 


reo52 


499 




can 9 51 


reo68 


499 


10 


can 959 


reD40 


499 




can 953 


reo78 


503 




can 954> 


reo52 


503 




con 955 


reD68 


503 




ran 95fi 


reo40 


503 


1 R 


can 957 


reo78 


510 




con 9 RR 


reo52 


510 




can 9 5Q 


ran68 


510 




can Ofin 




510 




can 9ft 1 




51 1 


9ft 


can 9R9 




511 




can 9fi*3 




511 


» 


can 9fi4. 


rpn.4-0 


51 1 




can 9fi 5 

acq. >J 


tpd78 


512 




can 9An 
btJq . iOO 




512 


9R 


can 9fi7 


rpn68 


512 




can 9fiR 

seq.zoo 


rpn4.0 


512 




can OfiQ 


rpn78 


516 




can 1~JC\ 
Seq . Z, / w 


rpn52 


516 




can 971 

seq .z / i 




516 




can 07° 

seq .z / z 


1 cp*rw 


516 




can 07*5 

seq . z / o 




517 




can Q7A 

seq . z. / *r 


pan R 9 


517 




can 07R 


rpn fi 8 


517 




seq. z / o 




517 




can 077 

seq. z / / 




517 




T"7Q 

seq.z /o 


ran 5 9 


517 




can 979 


reo68 


517 




seq. 280 


rep40 


517 




seq. 281 


rep78 


518 


40 


seq. 282 


rep52 


518 




seq. 283 


rep68 


518 




seq. 284 


rep40 


518 




seq. 285 


rep78 


519 




seq. 286 


rep 5 2 


519 


45 


seq. 287 


rep 68 


519 




seq. 288 


rep40 


519 




seq. 289 


rep78 


598 




seq. 290 


rep 5 2 


598 




seq. 291 


rep78 


598 


50 


seq. 29 2 


rep52 


598 




seq. 293 


rep78 


598 



TAC 

TAC 

GCT 

GCT 

GCT 

GCT 

GCC 

GCC 

GCC 

GCC 

GCG 

GCG 

GCG 

GCG 

GCA 

GCA 

GCA 

GCA 
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The above nucleic acid molecules are provided in plasmids, which 
are introduced into ceHs to produce the encoded proteins. The analysis 
revea | ec j the amino acid positions that affect Rep proteins activities. 
Changes of amino acids at any of the hit positions result in altered protein 

15 activity. Hit positions are numbered and referenced starting from amino 
acid 1 (nucleotide 321 in AAV-2 genome), also codon 1 of the protein 
Rep78 coding sequence under control of p5 promoter of AAV-2: 4, 20, 
22, 29, 32, 38, 39, 54, 59, 124, 125, 127, 132, 140, 161, 163, 193, 
196, 197, 221, 228, 231, 234, 258, 260, 263, 264, 334, 335, 337, 

20 342, 347, 350, 354, 363, 364, 367, 370, 376, 381, 389, 407, 411, 
414, 420, 421, 422, 424, 428, 438, 440, 451, 460, 462, 484, 488, 
495, 497, 498, 499, 503, 511, 512, 516, 517, 518, 542, 548, 598, 
600 and 601. The encoded Rep78, Rep68, Rep 52 and Rep 40 proteins 
and rAAV encoding the mutant proteins are provided. The corresponding 

25 nucleic acid molecules, Rep proteins, rAAV and cells containing the 
nucleic acid molecules or rAAV in which the native proteins are from 
other AAV serotypes, including, but are not limited to, AAV-1, AAV-3, 
AAV-3B, AAV-4, AAV-5 and AAV-6. 



30 175, 237, 250, 334, 429 and 519. 

Also provided are nucleic acid molecules, the rAAV, and the 
encoded proteins in which the native amino acid at each hit position is 
replaced with another amino acid, or is deleted, or contains additional 



Other hit positions identified include: 10, 64, 74, 86, 88, 101, 
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amino acids at or adjacent to or near the hit positions. In particular the 
following nucleic acid molecules and rAAV that encode proteins 
containing the following amino acid replacements or combinations 
thereof: T by N at Hit position 350; T by I at Hit position 462; P by R at 
5 Hit position 497; P by L at Hit position 497; P by Y at Hit position 497; T 
by N at Hit position 517; L by S at hit position 542; R by S at hit position 
548; G by D at Hit position 598; G by S at Hit position 598; V by P at Hit 
position 600; in order to increase Rep proteins activities in terms on AAV 
or rAAV productivity. The corresponding nucleic acid molecules, 
10 recombinant Rep proteins from the other serotypes and the resulting 
rAAV are also provided (see Figs. 5 and the above Table for the 
corresponding position in AAV-1, AAV-3, AAV-3B, AAV-4, AAV-5 and 
AAV-6) . 

Mutant adeno-associate virus (AAV) Rep proteins and viruses 
15 encoding such proteins that include mutations at one or more of residues 
64, 74, 88, 175, 237, 250 and 429, where residue 1 corresponds to 
residue 1 of the Rep78 protein encoded by nucleotides 321-323 of the 
AAV-2 genome, and where the amino acids are replaced as follows: L by 
A at position 64; P by A at position 74; Y by A at position 88; Y by A at 

» 

20 position 175; T by A at position 237; T by A at position 250; D by A at 
position 429 are provided. Nucleic acid molecules encoding these 
viruses and the mutant proteins are also provided. 

Also provided are nucleic acid molecules produced from any of the 
above-noted nucleic acid molecules by any directed evolution method, 

25 including, but are not limited to, re-synthesis, mutagenesis, recombination 
and gene shuffling and any way by combining any combination of the 

molecules, i.e., one, two by one, two by two, n by n, where n is the 

number of molecules to be combined ( i.e., combining all together). The 
resulting recombinant AAV and encoded proteins are also provided. 
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Also provided are nucleic acid molecule in which additional amino 
acids surrounding each hit, such as one, two, three . . . ten or more, 
amino acids are systematically replaced, such that the resulting Rep 
protein(s) has increased or decreased activity. Increased activity as 
5 assessed by increased recombinant virus production in suitable cells is of 
particular interest for production of recombinant viruses for use, for 
example, in gene therapy. 

Also provided are combinations of the above noted mutants in 
which several of the noted amino acids are changed and optionally 
10 additional amino acids surrounding each hit, such as one, two, three . . . 
ten or more, are replaced. 

The nucleic acid molecules of SEQ ID Nos. 563-725 and the 
encoded proteins (SEQ ID Nos. 1-562 and 726-728) are also provided. 
Recombinant AAV and cells containing the encoding nucleic acids are 
15 provided, as are the AAV produced upon replication of the AAV in the 
cells. 

Methods of in vivo or in vitro production of AAV or rAAV using 
any of the above nucleic acid molecules or cells for intracellular 
expression of rep proteins or the rep gene mutants are provided. In vitro 

20 production is effected using cell free systems, expression or replication 
and/ or virus assembly. In vivo production is effected in mammalian cells 
that also contain any requisite cis acting elements required for packaging. 

Also provided are nucleic acid molecules and rAAV (any serotype) 
in which position 630 (or the corresponding position in another serotype; 

25 see Figs. 5 and the table above) has been changed. Changes at this 

position and the region around it lead to changes in the activity or in the 
quantities of the Rep or Cap proteins and/or the amount of AAV or rAAV 
produced in cells transduced with AAV encoding such mutants. Such 
mutations include tgc to gcg change (SEQ ID No. 721). Mutations at any 
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position surrounding the codon position 630 that increase or decrease the 
Rep or Cap proteins quantities or activities are also provided. Methods 
using the rAAV (any serotype) that contain nucleic acid molecules with a 
mutation at position 630 or within 1, 2, 3 ... .10 or more bases thereof 
5 for the intracellular expression rep proteins or the rep gene mutants 

covered by claims 1 0 to 13, for the production of AAV or rAAV (either in 
vitro, in vivo or ex vivo) are provided. In vitro methods include cell free 
systems, expression or replication and/or virus assembly. 

Also provided are rAAV (and other serotypes with corresponding 

10 changes) and nucleic acid molecules encoding an amino acid replacement 
by N at Hit position 350 of AAV- 1 , AAV-3, AAV-3B, AAV-4 and AAV-6 
or at Hit position 346 of AAV-5; by I at Hit position 462 of AAV-1, 
AAV-3, AAV-3B, AAV-4 and AAV-6 or at Hit position 458 of AAV-5; by 
either R, L or Y at Hit position 497 of AAV-1, AAV-3, AAV-3B, AAV-4 

15 and AAV-6 or at Hit position 493 of AAV-5; by N at Hit position 517 of 
AAV-1, AAV-3, AAV-3B, AAV-4 and AAV-6 or at Hit position 535 of 
AAV-5; by S at hit position 543 of AAV-1 and AAV-6 or at hit position 
542 of AAV-3, AAV-3B and AAV-4 or at hit position 561 of AAV-5; by S 
at hit position 549 of AAV-1 and AAV-6 or at hit position 548 of AAV-3, 

20 AAV-3 B and AAV-4 or at hit position 567 of AAV-5; by either D or S at 
Hit position 599 of AAV-1, AAV-4 and AAV-6 or at Hit position 600 of 
AAV-3 and AAV-3B; by P at Hit position 602 of AAV-1, AAV-4 and AAV- 
6 or at hit position 603 of AAV-3 and AAV-3B or at hit position 589 of 
AAV-5 in order to increase Rep proteins activities as assessed by AAV or 

25 rAAV productivity. Methods using such AAV for expression of the 
encoded proteins and production of AAV are also provided. 
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Since modifications will be apparent to those of skill in this art, it is 
intended that this invention be limited only by the scope of the appended 
claims. 
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WHAT IS CLAIMED IS: 

1 . A process for the production of a peptide, polypeptide, or 
protein having a predetermined property, comprising: 

(a) producing a population of sets of nucleic acid molecules that 
5 encode modified forms of a target protein; 

(b) introducing each set of nucleic acid molecules into host cells 
and expressing the encoded protein, wherein the host cells are present in 
an addressable array; and 

(c) individually screening the sets of encoded proteins to identify 
10 one or more proteins that have activity that differs from the target 

protein, wherein each such protein is designated a hit. 

2. The process of claim 1, wherein each set of nucleic acid 
molecules is individually designed and synthesized. 

3. The process of claim 2, wherein each set is deposited at a 
15 locus in an addressable array. 

4. The process of claim 1 , wherein each polynucleotide in a set 
encodes a protein that differs by at least one amino acid from the target 
protein. 

5. The process of claim 1, wherein the array comprises a solid 
20 support with loci for containing or retaining cells; and each locus contains 

one set of cells. 

6. The process of claim 1, wherein the array comprises a solid 
support with wells for containing or retaining cells; and each well contains 
one set of cells. 

25 7. The process of claim 1, wherein the nucleic acid molecules 

comprise viral vectors; and the cells are eukaryotic cells that are 
transduced with the vectors. 

8. The process of claim 1, wherein the nucleic acid molecules 
comprise plasmids and the cells are bacterial cells. 
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9. The method of claim 1 , further comprising: 

(d) modifying the nucleic acid molecules that encode the hits, to 
produce a set of nucleic acid molecules that encode modified hits; 

(e) introducing the each set of nucleic acids that encode the 
5 modified hits into cells; and 

(f) individually screening the sets of cells that contain the nucleic 
acid molecules that encode the modified hits to identify one or more cells 
that encodes a protein that has activity that differs from the target protein 
and has properties that differ from the original hits, wherein each such 

10 protein is designated a lead. 

10. The process of claim 9, wherein each set nucleic acid 
molecules in step (d) is individually designed and synthesized. 

1 1 . The method of claim 1 , wherein the nucleic acid molecules in 
step (a) are produced by a method selected from among nucleic acid 

15 shuffling, recombination, site-directed or random mutagenesis and de 
novo synthesis. 

1 2. The method of claim 9, wherein the nucleic acid molecules in 
step (d) are produced by a method selected from among nucleic acid 
shuffling, recombination, site-directed or random mutagenesis, and de 

20 novo synthesis. 

13. The method of claim 2, wherein the nucleic acid molecules in 
step (a) are produced by systematically changing each codon in the target 
protein to a pre-selected codon. 

14. The method of claim 13, wherein the codon is selected from 
25 a codon encoding Ala (A), Ser <S), Pro (P) or Gly (G). 

15. The method of claim 13, wherein the codon is selected from 
a codon encoding Arg (R), Asn (N), Asp (D), Cys (C), Gin (Q), Glu (E), His 
(H), He {I), Leu (L), Lys (K), Met (M), Phe (F), Thr (T), Trp (W), Tyr (Y) or 
Val (V). 
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16. The method of claim 9, wherein the nucleic acids of step (d) 
are produced by systematically replacing each codon that is a hit, with a 
codon encoding the remaining amino acids, to produce nucleic acid 
molecules each differing by at least one codon and encoding modified hits 

5 to identify leads. 

17. The method of claim 9, further comprising: 
recombining the nucleic acid molecules encoding the leads; 
introducing those nucleic acid molecules into cells; and 

■ 

screening the cells to identify nucleic acid molecules that encode 

10 optimized leads. 

18. The method of claim 17, wherein the recombining is two, 
three or more up to all of the nucleic acids encoding the leads. 

19. The method of claim 17, wherein the recombining is effected 
by a method selected from among nucleic acid shuffling, recombination, 

15 site-directed or random mutagenesis and de novo synthesis. 

20. The method of claim 1 , wherein the modifications are 
effected in a selected domain of the target protein. 

* 

21. The method of claim 1, wherein the modifications are 
effected alqng the full length of the target protein. 

20 22. The method of claim 1 , wherein the change in activity is at 

least about 10%, 20%, 30%, 40% or 50%. 

23. The method of claim 1, wherein the change in activity is at 
least about 75%, 100%, 200%, 500% or 1000%. 

24. The method of claim 7, wherein at step (b) the titer of the 
25 viral vectors in each set of cells is assessed. 

25. The method of claim 24, wherein titering is effected by real 
time virus titering, comprising: 

(i) incubating the nucleic acid molecules or a vector 
(biological agent) comprising the nucleic acid molecules at an initial 
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concentration C, which is the unknown titer, with the host cells at a 
constant known concentration, D; 

(ii) measuring at successive times, an output signal, i; 

(iii) determining the time tfi, wherein: 
5 t/2 corresponds to i =£; 

£ min < yff < £ max ; 

jff min and /? max correspond to values of i at the 
inflection point of the curve i = f(t), for the minimal and maximal values, . 
respectively, of the concentrations of a reference biological agent for 
10 which the curve t£ = f(c) is predetermined; and 

(iv) determining the initial concentration C. 

26. The method of claim 24, wherein titering is effected by 
Tagged Replication and Expression Enhancement, comprising: 

(i) incubating with host cells a reporter- virus vector 
15 with a titering virus of unknown titer, wherein the titering virus increases 

or decreases the output signal from the reporter virus; and 

(ii) measuring the output signal of the reporter virus 
and determining the titer of the reporter virus; 

(ii) determining the titer of the interfering virus by 
20 comparing the titer of the reporter virus in the presence and absence of 
the interfering virus. 

27. The process of claim 9, wherein the nucleic acid molecules 
comprise viral vectors; and the cells are eukaryotic cells that are 
transduced with the vectors. 

25 28. The method of claim 27, wherein at step (f) the titer of the 

viral vectors in each set of cells is determined. 

29. The method of claim 28, wherein the target protein is a 
protein involved in viral replication. 
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30. The method of claim 1 , wherein the performance of the 
screened proteins is evaluated by a Hill analysis or by fitting the output 
signal to a curve representative of the interaction of the target protein and 
a test compound. 

5 31. The method of claim 30, wherein the Hill analysis, comprises: 

(a) preparing a sample of each nucleic acid molecule or a plasmid or 
vector that comprises each nucleic acid molecule (biological agent), 
wherein each sample is obtained by a serial dilution of the molecules or 
vector or plasmid at a concentration R1, 
10 (b) incubating each sample of the dilution obtained in (a) with the 

host cells (target cells) at a constant concentration R2, 

(c) determining a P product from the reaction R1 + R2, at a t 
moment, in each the sample; and 

(d) preparing a theoretical curve H from the experimental points R1 
15 and P, for each biological agent by iterative approximation of parameters 

of the reaction Rl + R2 -> P, at the t moment, in accordance with the 
equation: 

P = Pmax (TRU r / + ("R1) r ) r=1,...,n (2) 

in which: 

20 R1 represents the biological agent concentration in a sample 

from the scale; 

R2 is concentration of target cells (in vitro or in vivo) 
P (output) represents the product from the reaction R1 + 
R2 at a t moment; 
25 P max represents the reaction maximal capacity; 

k represents, at a constant R2 concentration, the biological 
system for responding to the biological agent (resistance cpnstant R2); 

r represents a dependent coefficient of R1 and corresponds 
to the Hill coefficient; and 
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rr represents the intrinsic power of the R1 biological agent to 
induce a response in the biological system (P production at the t 
moment), and 

{e) sorting the k and rr values obtained in (d) for each protein 
5 encoded by the nucleic acid molecules or plasmids or vectors and the 
cells, and then ranking according to the values thereof. 

32. The process of claim 1 that is automated. 

33. The process of claim 32 that is computer-controlled. 

34. A non-random method for generating proteins with a desired 
10 property, comprising: 

identifying a target protein; 

preparing sets of variant nucleic acid molecules that each encode a 
protein that differs by one amino acid from the target protein; 

screening and selecting the sets of variant nucleic molecules to 
15 identify those that encode proteins that have activity that differs by a 
predetermined amount from the activity of the target protein, thereby 
identifying proteins that are hits; 

identifying the residues in the hit proteins encoded by the variant 
nucleic acid molecules that differ from the target proteins; 
20 preparing further sets of variant nucleic acid molecules in which 

each codon in the nucleic acid molecule encoding each of the identified 
residues in each of the hits is replaced with codons encoding each of the 
remaining 1 8 amino acids to produce the further sets of variant nucleic 
acids, wherein in each set differs from each other set by one codon; and 
25 screening the further sets of nucleic acid molecules to identify 

those that encode proteins that have activity that differs from the activity 
of the hits, thereby identifying nucleic acid molecules that encode leads. 

35. The method of claim 34, wherein the replaced amino acid 
positions comprise a functional domain of the protein. 
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36. The method of claim 34, wherein the positions in the protein 

« 

in which amino acids are replaced comprise at least about 50% of the 
amino acids in the protein. 

37. The method of claim 34, wherein the positions in the protein 
5 in which amino acids are replaced comprise at least 90% of the amino 

acids in the protein. 

38. The method of claim 34, wherein the positions in the protein 
in which amino acids are replaced comprise at least 95% of the amino 
acids in the protein. 

10 39. The method of claim 34, wherein the positions in the protein 

in which amino acids are replaced comprise all of the amino acids in the 
protein. 

40. The method of claim 34, wherein each set of nucleic acid 
molecules is generated, processed and screened separately or in parallel. 
15 41. A method for producing a protein having modified properties, 

comprising: 

(a) preparing a population of nucleic acid molecules that encode 
rationally modified proteins; 

(b) inserting the population into expression vectors; 

20 (c) introducing each vector into host cells therefor, and expressing 

the modified proteins, 

(d) screening each modified protein, and selecting one or more that 
has (have) a modified property. 
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10 20 30 40 50 60 

1 MPGi^IVIJCVPSDIiDEHLPGXSDSFVSWAEKEWELPPDSDMDI^LISQAPLrVAEKZiQ 60 

2 MPGFYEIVIJCVPSDLDEKI^GXSDSFVNfcA/Aff 60 

3 MPGFYEIVIJCVPSDLDEHLPGXSNSFWWVAEKEWELPPDSDMDPNLISQ 60 

4 MPGCTEIVIJCVPSDLDEHLPGXSNSFWWAZK^ 6 0 

5 MPGFYEIVIJeVPSDLDEHLPGXSDSFVSfvrVASKEWELPPDSDrro^^^ 60 

6 MPGryEIVIXVPSDLDGHLPGXSDSFVNWASKEWELPPDSDMDLNLIEQAPLrVAEKL^ 60 

7 MATrYSVIVRVPFDVEEHLPGXSDSFVDWVTGQIWELPPESDLNLTLVSQPQLIVADRIR 60 
C M**fYE** : *VP*D***HLPGXS+SFV: WV****WEIiPP*SD** + *L*EQ**LrVA**** 

70 80 90 100 110 120 

1 RDFI,VQWRRVSKAPEALFFVQFEKG£SyFHLKILVETTGViCSMVLGRFLSQIRDKLVQTI 120 

2 RDFLVQWRRVSKAPEALFFVQFEKGESrFHLHILVETTGVJCSMVLGRFLSQIRDKLVQTI 120 

3 RE FLVEWRRVS KAPEALF FVQ FEKGETYFHLHVL I E T I GViCS MWGRYV SQIK2 KLVTR I 12 0 

4 RE FLVEWRRVS KAPEALF FVQ FE KGET YFHLHVL I ET I G V JCSMWGR YVS Q I KE KLVTR I 12 0 

5 RE FLVEWRRVS KAPE ALFFVQFEKGDS YFHLHI LVETVGVJCSMVVGRYVS Q I KE KLVTR I 12 0 

6 RDFLTEWRRVSKAPEALFFVQFEKGES5TFHMHVLVETTGVJCSMVLGRFLSQIREKLIQRI 12 0 

7 RVFL YE WNKFS KQ - ESKFFVQFEKGSE YFKLHTL VETSG I S SMVLGR YVS Q I RAQLVKW 119 
C R:FL++W***SK**E**FFVQFEKG+:YFH*H:L+ET:G**SMV:GR: : SQI : : *L* : :* 

130 140 150 160 170 180 

1 YRGXrPTLPNWFAVTKTRNGAGGGNKVVDECYIPNYLLPKrQPELQWAWTNMEEyiSACL 180 

2 YRGXEPTLPNWFAVTKTRNGAGGGNKVVDECYI PNYLLPKFQPELQWAWTNMEE YISACL 180 

3 YRGVSPQLPNWFAVTKTRNGAGGGNKVVDDCYI PNYLL P KTQPELQWAWTNMDQ YLS AC L 180 

4 YRGVJSPQLPNWFAVTKTRNGAGGGNKWDDCYI PNYLLPKTQPELQWAWTNMDQ YLSACL 180 

5 YRGVSPQLPNWFAVTKTRNGAGGGNKVVDDCYI PNYLLPKTQPELQWAWTNMDQ YI SACL 180 

6 YRGXEP Th PNW^AVTKTRNGAGGGNKWDEC YI PNYLL P KTQPELQWAWTNMEQ YLS AC L 180 

7 FQGJEPQINDWVAITKVKK- -GGANKVVBSGYIPAYLLPKVQPELQWAWTNLDEYKLAAL 177 
C * *G : E9 : * * * W*A*TK* * * * *GG*NKWD : * YIP* YLLPK*QPELQWAWTN* : : Y : * A* L 

190 200 210 220 230 240 

1 NIiAERKRLVAQHLTHVSQTQEQNKENLNPNSDAPVIRSKTSARYMELVGWLVDRGirSEK 24 0 

2 NLAERKRLVAHDLTHVSQTQEQNKENLNPNSDAPVIRSKTSARYMELVGWLVDRGirSEK 240 

3 NLAERKRL VAQHLTHVSQTQEQNKENQNPNSDA PVI RS KTS ARYMS L VGWL VDRG I TS EK 24 0 

4 NLAERKRLVAQHLTHVSQTQEQNKENQNPNSDAPVIRS KTSARYMEL VGWL VDRG I TSEK 24 0 

5 NIAERK^VAQHLTHVSQTQEQNKENQNPNSDAPVIRSKTSARYMELVGWLVDRGITSEK 24 0 

6 NLTERKRLVAQHL THVS QTQEQNKENQNPNSDAP V IRS KTSARYME L VGWL VDiCG I TS E K 24 0 

7 NLEERKRLVAQFLAESSQRS-QEAASQREFSADPVIKSKTSQKYMALVNWLVEHGIXSEK 236 
C NL+ERK21LVA*+L***SQ***Q****+***S**PVI*SKTS**YM*LV*WLV*+GIXSEK 

250 260 270 280 290 300 

1 Q W I QEDQAS YI S FNAAS2JSR S QXXAALDNAGKI MALTKS A PD YL VG PAP PAD I KTNRI YR 300 

2 QWIQEDQASYISFNAASNSRSQXJCAALDNAGKIMALTKSAPDYLVGPAPPADIKTNRIYR 3 00 

3 QW I QEDQAS YI S FNAAS^SRSQXJCAALDNAS KIMSLTKTAPDYLVGSNPPED I TKMR I YQ 3 0 0 

4 QWIQEDQASYISFNAAStfSRSQXXAALDNASKIMSLTKTAPDYLVGSNPPEDITKNRIYQ 3 00 

5 QWI QEDQAS YI S FNAAS2tfSRS QXICAALDNAS Kl MS LTKTAPDYLVGQNP PED I S SNR I YR 3 00 

6 QW I QEDQAS YI S FNAAStfSRS QXJCAALDNAGKIMS LTKTAPDYLVGQQP VED I SSNRIYK 300 

7 QWIQENQESYLSFNSTGWSRSQXJCAALDNATKIMSLTKSAVDYLVGSSVPEDISKNRIWQ 2 96 

C QWIQE*Q*SY*SFN***tfSRSQXJCAALDNA:KIM+LTK:A*DYLVG: :**+DI: :NRI*: 

310 320 330 340 350 360 

1 ILELNGYEPAYAGSVFLGWAQK^FGIO^IWLFGPArTGKri/IAEAXAHAVPFYGCVNWT 3 6 0 

2 I LEIiNGYDPAYAGSVFLGWAQKRFGKRNTIWLFGPATTGKrWI AEAXAHAVPFYGCVNWT 3 60 

3 I LEIjNGYDPQYAASVFLGWAQKKFGKRNTIWLFGPATTGKI^I AEAXAHAVPFYGCVNWT 36 0 

4 ILELNGYDPQYAASVFLGWAQKKFGKRNTIWLFGPATTGKrinLAEAXAHAVPFYGCrV^ 360 

5 ILEMNGYDPQYAASVFLGWAQKKFGKROTIWLFGPArTGKIT/IAEAXAHAVPFYGCVNWT 360 

6 ILEIJKTGYDPQYAASVFLGWATKKFGKRNTIWLFGPArTGKr^ 3 6 0 
-n I FEMNGYDP AYAG S I LYGWCQRS FNKRNTVWL YGPATTG KZ!NT AE AXAHTVP F YGCVNWT 35 6 

C I * E +NGY *P:YA:S***GW***:F* KRNT * WL * GPATTGKTOT: AE AXAH+ VP F YGCVNWT 
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370 380 390 400 410 420 

1 NENFPFWDCVDKMVlWWEEGJOifrAKVVESA 420 

2 NENFP FND C VDKMV I WWEEGJCMTAKVVESAKAI LGGS KVRVDQKCKS SAQXD PTPVI VTS 420 

3 NEWFPFm)CVDKMVIWWEEGiOCTAKW 420 

4 NENFP FMD C VDKMV I WWEEGiCMTAKVVESAKAI LGGS KV^ 420 

5 NENFPFMDCtaDKKVIWWEEGiOifrAKVVE^ 420 

6 NENFPFJTDCVDKMVIWWEEGiCOTAKVVESAK 420 

7 NEOTPFNDCVDKMLI&WEEGKM^ 416 
C NEWFPF1^CVTDKM*IF^EEGJCMT*KVVESAKAILGGSKVR 

430 440 450 460 470 480 

1 OTNMCAVXDGNSTTFEHQQPLQDRMFKFELrRRLEHDFGJCVrKQEVKEFFRWAQDHVTEV 4 8 0 

2 imDtfCAVXjDGNSTTFEHQQPLQDRMFKFELrRRLEHDFGJO/rKQEVKEFFRWAQDHVTEV 4 80 

3 tfTNMC AVXJDGNS TT FEHQQPLQDRMFKFELTRRIiDHDFGiCVrKQE VKD F FRWAS DHVTDV 4 80 

4 jmiMCAVXIX^STTFEHQQPLQDRMFKFELrRRLDHDFGi^ 4 80 

5 J/TNMCAVXDGNS TTFEHQQPLQDRMF KFE LTKRLEHD FGitVITCQE VKD FFRWASDHVTE V 4 8 0 

6 OTNMCAVXI)GNSTTFEHC?QPLQDRMFKFELrRRLDHDFGiCVrKQE VKDFFRWAKDHWEV 4 80 

7 NTNWCVVVPGNSTTFEHOQPLEDRMFKFELrKRLPPDFGJCI ZTCQE VKDFFAWAKVNQVPV 476 
C N7W4C * V*Z>GNSTTFEHQQPL *DRMFKFELZ*+RL : *DFGJC*rKQEVK+FF*WA :***+: V 

490 500 510 520 

1 AHEJTYVRKGGANKRPAPD.DAD KSE P KRA CPS VM PS TSDAEG 522 

2 AHEFYVRJCGGANKRPAPDDADKSE PKRA C PS VAD PSTSD AEG 52 2 

3 AHEJ jr YVRJCGGAKKRPASNDADVSEPKRQ CTSLAQPTTSDAEA 522 

4 AHRFYVRJCGGAKKRPASND AD VSE P KRQ CTSLAQPTTSDAEA 522 

5 THRFYVRiCGGARKRPAPNDADISEPKRA CPSVAQPSTSDAEA 522 

6 EH EFYVKfCGG AKKRP APSD AD ISEP KR V - RESVAQPSTSDAEA 522 

7 THEFKVPRELAGTKGAEKSLICRPIjGDVTNTSYKSLSKRARLSFVPETPRSSDVTVDPAPL 53 6 
C :HEJ?*V+***A: ***A ::***.♦**** : +:*:*:*: ***A* : 

530 540 550 560 570 580 

1 APVDFADRYQNKCSRHAGMLQMLFPCKTCERMNQNFNICFTHGTRDCSECFP- -GVSESQ 580 

2 AP VD F ADR YQNXC S RKAGMLQMLF PCKTCERMNQNFN I C r THGTRDC S E C F P - -GVSESQ 580 

3 P-ADYADRYQ^CSRHVGMNLMLFPCKTCERMNQISNVCFTHGQRDCGECFPGMSESQPV 581 

4 P - ADYADRYQNKCSRHVGMNI/MLFPCKTCERMNQISNVCFTKGQRDCGECFPGMSESQPV 581 

5 P - VDYADRYQNKCSRKVGM^ILMLFPCRQCERMNQNVDI CFTHGVMDCAECFP - VSESQPV 580 

6 S-INYADRYQNKCSRHVGMNLMLFPCRQCERMNQNSNICFTKGQKDCIiECFP- -VSESQP 579 



7 RPLNWNSRYDCKCDYHAQFDNISNKCDECEYLNRGKNGCICHNVTKCQICHG 588 

C : : : + : **RY**KC**K : ** : : ****C : :CE**N* : : * : C**H* : : *C. *C** . . : : + : : : 

590 600 610 620 

1 PWRKRTYRKLCAIHHLLGRAPEIACSACDIiVNVDLDDCVSEQ 62 3 

2 PWRKRTYRKLCAIHHLLGRAPEIACSACDLVNVDLDDCVSEQ 623 

3 S WKKKTYQKLC P IHKI LGRAPE IACS ACDLANVDLDDCVS EQ 624 

4 SWKKKTYQKLCP IHH I LGRAPE I ACS ACDLANVDLDDC VS EQ 624 

5 SWRKRTYQKLCPIHHIMGRAPEVACSACELANVDLDDCDMEQ 623 

6 VSWKKAYQKLCYIHHIMG-KVPDACTACDLVNVDLDDCIFEQ 621 

7 IPPWEKENLSDFGDFDDANKEQ 610 

C : * + + + + : ** + **+* : **D*DD* : : EQ 
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